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SECTION  1 


INTRODUCTION 


This  technical  report  summarizes  the  work  performed  on 
>JCEL  contract  number  M6247‘i-a7-C-3056 ,  "Determination  of  th 
Best  Gr ound- Penet r at ing  Radar  Source  Signal  Type  for  Accura 
location  of  Underground  Utilities."  The  report  is  di video, 
several  sections  describing  the  work  performed  on  the  pro.cr 
and  the  results  obtained. 


There  are  many  applications  today  for  ground-penetrating 
radar  ( G?R )  systems  vith  varying  requirements.  These 


applications  encompass  utility  pipe  detection,  subsurface 
highway  surveys,  archaeological  surveys,  mine  shaft  detection 
and  more.  ^■lany  G?R  systems  have  been  developed  to  meet  these 


1  H 


s.  The  designs  of  the  G?R  systems  are  as  varied  as 


applications . 

Each  application  requires  that  the  G?R  provide  informat: 


about  the  subsurface  target(s)  of  interest  and  that  the  sign: 
from  the  desired  target(s)  be  distinguishable  from  other 
signals,  for  exaxaple,  from  those  generated  by  the  surrounding 
media.  Thus,  the  potential  G?R  user  identifies  system 


r f ormance 

requirements 

such  as  depth  penet 

ration,  target 

solution , 

range  accuracy,  measurement  spee 

d,  power 

nsumpt ion , 

,  and  the  GPR 

physical  di.mensions 

.  These 

rf ormance 

r  equi  r  eme.nts 

translate  into  the 

radar  system  des 

parameters  of  bandwidth,  pulse  width,  resolution,  range 
accuracy,  maximum  range,  pulse  repetition  frequency  (PR?), 
required  radiated  power  levels,  permitted  radiated  power 
levels,  measurement  scan  soesds ,  eoc . 


All  radar  designs,  whether 


nc' 


Ley  are  use; 


aoo 1  cat 


3ns  ,  i.ncoroor at 


well-established  signal  orocessi.ng 


.ec.nni  ques 


•Uost  of  t.he  basic  radar  parameters  are  determi.n: 


y  using  Fourier  Transform  analysis.  Thus,  radar  system 


:  s  1  g: 


■  s  mu: 


aooroacn 


is  ign 


:uirements 


nsidering  time-dcmain  and  f requency-dcmain 


umen' 


e  aoorca: 


3  m 


ind 


•  Ti3. 


insect 


sy  s  ■ 


^0  .n 


h a ve  beer,  built 
techniques,  and 

Most  of  the  time-domain  G?R  systems  use  baseba.nd  pulsers 
vich  pulse  widths  and  power  levels  determined  by  the  design 
requirements.  However,  optimum  pulse  shapes  are  difficult  to 
enerate.  These  time-domain  systems  often  use  simple, 
roadband  sampling  receivers  which  produce  low-frequency  cutp 
waveforms  that  are  easily  recorded  and  displayed.  Many 
different  signal -processing  techniques  have  been  applied  to 
cime-icmain  G?R  output  data  to  enhance  target  detection  and 
1  dent  if i cat  ion . 

Fr equency-dcmai.n  G?R  systems  typically  use  RF  frequency- 
modulation  techniques.  These  systems  are  usually  larger  and 
m.ore  expensive  than  t.heir  time-domain  counterparts.  However, 
the  f r equency-dc.mai.n  systems  offer  more  flexibility  in  RF 
bandwidth  selection.  Furthermore,  receivers  and  transmitters 
can  be  designed  that  yield  better  overall  system  clutter  and 
noise  performance  than  their  time-dcm.ain  counterparts.  Cutpu 
data  from  the  frequency-domain  radars  must  usually  undergo 
complex  signal  processing  before  it  can  be  easily  interpreted 
by  the  operator. 

The  purpose  of  this  GPR  signal  source  evaluation  progra.m 
was  to  identify  which  type  of  GPR  system  provided  the  best 
performance,  based  on  the  criteria  of  depth  penetration, 
resolution,  and  sig.nal -to-clutter  performance.  GAR  tested  tw 
pulsed  GPR  systems,  as  well  as  a  F?^-CW  based  GPR  system,  and 
stepped  FM-based  system  to  determi.ne  which  technology  yielded 
t.ce  best  results.  Each  of  the  systems  was  tested  on  a  scale- 
mcdel  test  tank  containing  homogeneous  soil  with  known  buried 
targets  and  probes.  Two  different  soil  types  were  modeled  wi 
the  tank  to  simulate  a  medium-loss  environ.ment  and  a  low-loss 
environment. 


using  time-domain  techniques,  f r equency-dcmai 
combinations  of  both. 


o 

d. 


U)  tl. 


SECTION  2 


SUMMARY  OF  WORK  PERFORMED 


Early  work  on.  the  program  involved  the  specification  cf 
the  candidate  scale-model  GPR  source  signal  types,  the  desigr 
of  a  suitable  soil  scale-model  tank  for  testing  the  GPR 
systems,  the  analytical  modeling  of  the  GPR  pipe  detection 
problem,  and  the  fabrication  of  antennas  to  be  used  fcr  test; 
one  candidate  GPR  systems.  The  remainder  cf  the  prcgra.m  vcrl- 
included  soil-tank  measurements  vith  each  of  the  candidate  T; 


sysoems,  the  examination  of  the  corresponding  data  collec 


from  these  measurements,  and  the  evalual 
obtained  from  each  of  the  candidate  systems 


the  perfcrmanc 
A  r eccmmendat i c 


:ased  on  this  evaluation  vas  made  which  specified  the  cotimr. 


GPR  source  signal  for  accurate  utility  location.  Finally,  a 
full-scale  version  of  the  GPR  system  producing  the  best  scale 
model  performance  was  qualitatively  tested  on  an  outdoor  pipe 
test  field  at  the  GAR  facility. 


Four  types  of  GPR  systems  u-ere  selected  for  testing. 


systems  were  time-domain  short-pulse  radars;  one  was  a  swep- 
frequency  FM-CV  system,  and  one  was  a  stepped-FM  system. 
four  systems  represent  radar  technology  that  had  been  appli; 
to  GPR  problems  in  the  past.  One  of  the  short-pulse  radar 
systems  is  built  by  Gulf  Applied  .Research  (GAR)  and  utilize: 
cne-nanosecond  mcnocycle  pulser  and  an  equivalent-time  samp: 


’he  other  short-pulse  radar  is  built  by  Gecphysica 


Survey  Systems,  Inc.  (GSSI)  and  uses  a  balanced,  one-nanose; 
unipolar  pulser  and  a.n  equivalent-time  sampling  receiver. 


two  f r eque.ncy-dc.mai .n  syste.ms  were  simulated  using  a  Hewlett- 


Packard  35  ’  0  .network  analyzer  .  The  .network  analyzer  was 


operated  i.n  its  t  i.me -domai.n  mode  with  synthesized  one- 
nancseco.nd  unipola.'’’  pulses.  Thus,  the  performance  of  all 
systems  could  be  app  r  cpr  i  at  e  ly  compared  by  exa.mini.ng  their 
t  ime-dcmai.n ,  or  si.mulated  ti.me-dcraain  responses.  Function 
escriptions  of  the  candidate  GPR  systems  are  included  in 
e  c  1 1  o  n  6  . 


’he  pipe  target  detection  problem  was  modeled  analytical 
ier  to  evaluate  the  expected  performance  of  GPR  system.s  . 


apprcxim.ate  a.nalytical  m.odel  based  on  the  radar 
was  used.  The  model  accounted 


uat  1 ' 


o  r  a.n  ^ 


n  a  oar  am.  e  t  e  r 


3 


soil  parameters 


and  the  pipe  target  parameters.  The  analysis 
igh  signal  losses  were  expected,  for  deeply 
n  moderately  lossy  soil.  The  predicted  signal 


less  was  higher  than  expected.  GAR  concluded  that  the 
appr oximat i ons  used  in  the  analytical  model  may  have  introduce 
seme  error  into  the  calculations.  Actual  results  cctai.ned 
while  testing  the  GPR  systems  on  the  soil-model  tank  indicated 
that  the  analytical  model  indeed  produced  pessi.mistic  results. 


We  concluded  that  accurate  measure.ments  of  the  G?F.  antenna 
pattern  in  the  soil,  ana  of  the  radar  cross  sections  of  various 
pipe  targets,  would  facilitate  the  improvement  of  the 
analytical  model.  Some  of  the  results  from  the  GPR  tests  i.o 
this  program  could  be  used  to  i.mprove  the  analytical  model. 
However,  an  extensive  measurement  program  to  establish  an 
accurate  analytical  model  was  beyond  the  scope  of  this  progra-ti. 

GA.R  elected  to  build  a  scale-model  soil  environment 
containing  buried  pipe  a.nd  plate  targets  for  the  purpose  of 
testing  the  oa.ndidate  GPR  systems.  This  was  dons  because  a 
full-scale  testing  area  with  buried  targets  would  be  difficult 
to  construct  and  maintain.  The  design  of  the  scale-model  soil 
tank  is  discussed  in  Section  5.  Some  of  the  important  concepts 
and  factors  influencing  the  design  of  the  scals-mcdel  tank  are 
outlined  below. 


so 


*  0 
me 


.me 

3.1 


S  3 

de 


Essentially,  scaling  allows  the  reduction  of  all 
nsions  of  the  full-scale  environment  by  the  corrsspondi.n 
e  factor.  This  i.ncludss  all  target  dimensions,  burial 
hs ,  and  the  pulse  widths  of  the  GPR  transmitter  signals, 
efore,  the  dime.nsions  of  tne  test  facility  and  the  GPR 
nnas  were  reduced  by  the  scale  factor.  The  electromagne 


nnar  ac  t  e  r  1  s  1 1  c  c: 


;he  diel 


material  used 


'  r.  R 


*  tanx  was  scaied  nigher 
1  .n  the  full-scale  e  n  v  i  r  ■ 
same  total  soil  loss  ove: 


;nan 

..ment  .  This  was  done  to 
tne  scale-model 


.ne  loss  characteristic  oi 

ainta: 


tnat 
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GAR  initially  cor 
ta.nx  using  a  li 
e  V  e  r  a  1  a  d  v  an 


s  i  de  r 
:uid  ; 
lages 

r>  ** 

^  ^  w  * 

'j  1  i  1 


tier,  or  a  soi*  sca*e- 
A  liquid  dielectric 

^  W 


;ver  a  sci*  aieiectric. 

:  liqu:.d  could  be  accura 
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ddp  ...  ^ 


1  n  s  e  r  t  j 


G  . 


me ve c  ; r 


4. 


the  liquid,  and  the  electromagnetic  effects  cf  the  tank' 
b cun dan ies  could  be  r emoved  from  the  G?R  test  data.  Despite 
these  advantages,  a  liquid  soil-model  tank  would  have  been 
costly  to  build  and  maintain.  Also,  there  was  a  significa.nt 
technical  risk  that  the  liquid  dielectric  properties  could  net 
be  mai.ntained  with  satisfactory  homogeneity  in  a  large  volume 
tank,  such  as  the  cne  that  would  have  been  required  for  the  G? 
tests  . 


cost  cin 


3 


)rcblems  associated 


maintaining  a  liquid  sca-e-moue 
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prompted  GAR  to  use  hcmogenecus  soils  as  scale-model  dielectri 
materials  instead.  G.AR  located  commercial  sources  of  two  loss 
clay  soils  and  used  these  soils  as  the  scale-model  dielectrics 


These  clays  exhibited  high  enough  attenuation  factors 


;o  a-  - ; 


;  .n  e  m 


le  used  as  scaled  dielectrics. 


GAR  filled  the  soil-model  tank  with  one  of  the  clays  and 
buried  various  plastic  pipes,  .metal  pipes,  and  metal  plates  in 
the  clay  at  different  depths.  After  measurements  were  made 
wioh  the  first  clay  soil,  the  tank  was  emptied  and  re-filled 
with  the  other  clay.  The  same  targets  were  buried  in  the  sa.me 
locations  as  they  were  with  the  original  clay,  and  measurement 
were  then  made  with  the  new  soil.  Using  two  different 
materials  in  the  test  tank  allowed  the  prediction  of  GPR 
performance  in  two  full-scale  environments.  This  approach  als 
facilitates  the  extrapolation  of  the  results  obtained  to  other 
similar  soil  types. 

Antennas  were  specially  designed  to  be  used  for  the  tests 
of  the  candidate  G.^.R  syste.ms.  .Accurate  comparisons  betwee.n  t.o 
G.^.R  syste.ms  were  possible  only  if  the  same  antennas  were  used 
for  eac.n  of  them.  GAR  designed  and  built  a  pair  of  TEM  horn 
a.ntennas  and  a  pair  of  shielded  triangular-sheet  dipole 
a.ntennas  to  be  used' with  the  radars.  The  performance  cf  each 
pair  of  anten.nas  was  evaluated  and  the  best  pair  was  used  for 
the  GPR  ocmparisc.n  tests.  Section  3.3.T  establishes  the  desig 
requirements  of  GPR  ante.nnas  and  Section  6.4-  describes  the  two 
types  of  a.nten.nas  constructed  by  GAR  for  the  tests. 

T.he  TEM  antennas  had  .more  desirable  pulse  performa.nce 
characteristics  tnan  the  dipole  antennas.  Therefore, 
considerable  care  was  taken  in  their  fabrication.  I.n  addition 
t.'ie  TE.M  horns  were  designed  a.nd  built  very  carefully  i.n  a.n 
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attempt  to  minimize  the  effects  of  clutter  signals  on  the 
received  target  signals.  However,  we  discovered  during  the 
tests  of  the  two  types  of  antennas  that  the  dipole  antennas 
provided  noticeably  better  signal-tc-clutter  performance  than 
the  TEM  horns.  Since  the  signal-to-clutter  ratio  was  the 
dcmi.nant  limitation  in  target  detectability,  the  triangular- 
sheet  dipoles  were  used  for  the  subsequent  GPR  comparison 
tests  . 

The  four  candidate  G?R  sources  were  tested  with  the  two 
clay  model  environments .  The  tr iangular -sheet  dipole  antennas 
were  used  wizh  all  of  the  systems  except  the  GSSI  radar.  The 
antennas  could  not  be  used,  with  the  GSSI  radar  because  it  had  a 
single  integrated  receiver,  transmitter,  and  antenna  asse.mbly. 
Target  signal  return  amplitudes  were  measured  from  the  pipe 
targets  and  the  performances  of  the  radars  were  compared.  The 
measurement  results  are  presented  in  Section  8. 

Section  3  also  presents  the  results  of  measurements  using 
a  full-scale  version  of  the  scale-model  GPR  system  that  gave 
the  best  performance  in  the  scale-model  tests.  The  system  was 
tested  on  an  outdoor  pipe  test  field  containing  buried  plastic 
and  metal  pipes.  Oualit'ative  data  from  the  full-scale 
measurements  demonstrated  the  general  performance 
characteristics  of  that  GPR  system. 

A  recommendation  for  a  high-performance  GPR  based  on  an 
evaluation  of  the  measurement  results  is  given  in  Section. 9. 
After  reviewing  the  results,  GAR  felt  that  a  single,  fixed  RE 
bandwidth  GPR  is  not  flexible  enough  to  meet  all  of  the 
requirements  of  resolution  and  deep  penetration.  Several  RE 
candwidths  are  required.  The  frequency-domain  systems  provide 
more  performance  flexibility  than  the  time-domain  pulsed  GPR 
systems  at  the  expense_of  higher  cost  and  complexity.  If  the 
receiver  dynamic  range  and  system-ge.nerated-clutter  performance 
e  used  as  evaluation  criteria,  a  pulsed,  ti.me-dcmain  GPR 
cannot  typically  match  the  capabilities  of  a  frequency-domain 
T?R  that  utilizes  a  synthesized  signal  source.  Reducing 
clutter  and  maximizing  dynamic  range  are  critical  in  order  to 
cistinguish  desired  target  signals  from  clutter  returns. 


A 


SECTION  3 


DESCRIPTION  OF  THE  GROUND  PENETRATION  PROBLEM 

Ground-penetrating  radars  are  designed  to  operate  in  a 
considerably  different  environment  than  more  conventional 
radars.  The  GPR  must  operate  with  a  medium  that  is 
electrcmagnetical ly  inhomogeneous  and  lossy.  Furthermore,  the 
electromagnetic  properties  of  the  soil  containing  one  targecs 
are  usually  unknown  or  are  known  only  near  the  surface.  These 
properties  may  vary  on  a  daily  basis  with  varying  moisture 
content . 

An  ideal  ground-penetrat ing  radar  system  should  be  able  t 
detect  targets  buried  at  any  desired  depth,  and  resolve 
multiple  targets  if  they  are  closely  spaced.  However,  the 
relative  dielectric  constant  and  attenuation  properties  of  the 
soil  generate  severe  limitations  on  one’s  ability  to  construe" 
a  G?R  system  that  provides  such  ideal  performance.  The  GPR 
design  goals  usually  represent  a  compromise  between  such 
desired  ideal  performance  and  the  best  performance  achievable 
under  given  conditions.  For  example,  a  high  resolution  GPR 
typically  is  not  capable  of  great  depth  penetration.  In 
addition,  most  GPR  systems  have  maximum  range  capabilities  of  . 
few  meters  or  less.  Although  the  radar  designer  is  typically 
not  able  to  achieve  ideal  performance,  he  can  in  practice 
approach  optimum  performance,  that  is,  the  best  performance 
obtainable  with  available  technology  under  existing  conditions. 


3.1  SOIL  DIELECTRIC  PROPERTIES 

The  di e 1 ectr ic  'proper t ies  of  the  soil  affect  the 
performance  of  a  GPR  more  than  any  other  factor.  These  soil 
propercies  vary  with  geographical  location,  as  well  as  with 
depth  from  the  surface.  The  attenuation,  conductivity, 
relative  dielectric  constant,  loss  tangent,  and  permeability 
are  interrelated  factors  that  determine  the  way  electromagnet: 
energy  propagates  through  the  soil.  The  effects  that  these 
soil  characteristics  have  on  the  performance  requirements  cf  a 
GPP.  are  explai.ned  below. 
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3.1.1  ELECTROMAGNETIC  PROPAGATION  IN  SOIL 

Th_  propagation  characteristics  of  an  electromagnetic  wave 
in  any  medium  are  determined  largely  by  the  propagation 
constant  of  the  medium.  The  propagation  constant  includes  the 
effects  of  attenuation,  conductivity,  and  the  relative 
dielectric  constant.  The  following  description  of 
electromagnetic  wave  propagation  in  soil  illustrates  the 
effects  of  these  factors. 

To  simplify  the  description  of  electromagnetic  wave 
propagation  in  soil,  or  amy  other  medium,  it  is  helpful  to  make 
several  assumptions.  Let  us  assume  that  the  soil  medium  is 
infinite  in  spatial  extent,  linear,  homogeneous,  isotropic,  and 
electromagnet ically  source-free.  Assume  also  that  the 
electromagnetic  field  energy  propagating  in  the  ground  has  a 
transverse  electromagnetic  field  structure  (the  TEM  mode),  has 
a  sinusoidal  variation  with  time,  and  operates  as  a  plane  wave. 

Under  these  assumptions.  Maxwell’s  equations  may  be 


expressed  as: 

V2  E 

- 

7^ 

E  = 

0 

(3.1) 
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72 

H  = 

0 

(3.2) 

V  . 

H 

= 

0 

(3.3) 

V  . 

E 

= 

0 

(3.4) 

7  = 

(j' 

u;  n 

(  <7  + 

(3.5) 

Here,  S  and  H  are  vector,  sinusoidal  representations  of  the 
electric  and  magnetic  field  components  of  the  electromagnet i o 
wave.  The  scalar  quantity  7  is  the  complex  propagation 
constant,  determined  by  characteristics  of  the  medium.  It  is 
described  in  detail  later.  The  term  ^J.  is  the  permeability  of 
the  medium,  a  is  the  conductivity,  and  €  is  the  dielectric 
constant . 
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Assume  that  a  rectangular  coordinate  system  is  establishe 
having  x,  y,  and  z  axes,  with  the  positive  z  direction  being 
the  direction  of  propagation  of  the  electromagnetic  wave  into 
the  ground.  This  assumption  yields  solutions  of  Maxwell’s 
equations  of  the  form: 


E 


X 


e 


(3.6 
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H.  =  H,  e 


72 


( 3.7 
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Where  E^ ,  ,  and  are  scalar  components  of  the  vectors  z 

and  H  in  the  x  and  y  directions.  The  +  sign  on  the  E  and  H 
field  components  indicates  that  the  wave  is  traveling  in  the 
positive  z  direction.  Solutions  also  exist  for  waves  traveli.n 
in  the  negative  z  direction,  but  they  are  neglected  here  for 
simplicity.  Equations  (3.6)  through  (3.9)  describe  the  form  c 
the  electromagnetic  plane-wave  components  in  the  soil  with  the 
parameters  established  above. 


The  complex  variable  7  in  the  expo 
through  (3.9)  is  of  great  interest  in  p 
performan.cs  of  a  GPR .  This  term  is  def 
and  contains  all  of  the  electromagnetic 
medium  through  whiciq  the  wave  propagate 
=  27rf,  where  f  is  tTie  frequency  of  the 
electromagnetic  energy.  Also,  j  =  (-1  ) 
imaginary  number . 


ne.nt  of  Equations  (3.5 
redicting  the 
ined  in  Equation  (3.5) 
parameters  of  the 
s.  In  Equation  (3.5). 
propagating  sinusoidal 
and  denotes  an 


Often , 
parts  as  in 


7  is  expressed  in  terms 
Equation  (3.10)  below: 


of  its  real  and  imaginary 


7  =  a  + 


(3.1G 
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Solving  Equation  (3.5)  for  its  real  and  imaginary  parts 
produces  two  roots;  however,  only  the  solution  yielding  a 
positive  imaginary  part  is  valid  for  real  media.  The  quantity 
a  is  the  attenuation  factor  and  ,3  is  the  phase  propagation 
factor . 


Once  the  quantities  /j,  ,  a,  and  6  are  determined  for  a 
given  medium,  the  attenuation  and  phase  factors  in  Squuf'on 
(  3  .  3  )  are  known.  Solving  for  a  and  3  yields: 


a 


UJ  fl  e 


■'Z  1  +  {a / u  e 


(3.11) 


3 


=  U  yflU 


^  1  +  ( cr/  w  £  )^ 


(3.12) 


Often  the  electromagnetic  soil  attenuation  characteristic  is 
expressed  in  d3/meter.  This  method  of  expressing  the  soil 
attenuation  addresses  power  loss  and  is  computed  as  follows; 


Attenuation  (dB/meter)  =  20  log,0  e“  . 


(3.13) 


Other  quantities  used  to  characterize  dielectric  media  and 
soil  are  derived  from  q,  ij, ,  a,  and  e.  These  quantities  are 
the  relative  dielectric  constant,  the  relative  permeability, 
the  loss  tangent,  and  the  intrinsic  impedance.  The  definitions 
of  these  quantities  are  given  below: 

Relative  Dielectric^ Const,  =  f/e^  (3.14) 

Relative  Permeability  (3.15) 

Loss  Tangent  5  =  a / e  (3.16) 


Intrinsic  Impedance  q  =  +  ju  e  ))'''^ 


(3.17) 


Here,  and  fj  ^  are  the  permittivity  and  permeability  of  free 
space,  respectively.  For  most  soils,  ^^=1  . 

The  intrinsic  impedance  relates  the  amplitudes  of  the  E 
and  H  field  components  by  the  equation: 

U  I  =  I  E  1  /  [H  I  (3.1S) 

In  this  report,  the  most  commonly  used  terms  to  describe  soil 
characteristics  will  be  a,  the  attenuation  in  dB/m,  and 


3.1.2  DETERMINING  SOIL  PARAMETERS 

One  of  the  most  challenging  problems  confronting  any  GPR 
designer  is  to  design  the  radar  to  operate  in  a  soil 
environment  that:  Is  quite  diverse.  The  designer  must  address 
the  most  demanding  environment  and  tailor  his  design  to  meet 
those  requirements.  The  worst  soils  often  exhibit  high 
relative  dielectric  constants  coupled  with  high  less 
characteristics.  These  lossy  soils  impose  stringent 
requirements  on  maintaining  high  signal -to-clutter  and  signal- 
to-noise  ratios  in  the  radar. 

Soil  characteristics  vary  greatly  with  geographical 
location.  Figure  3.1  is  a  map  of  the  United  States  indicating 
the  different  soil  types  that  may  be  found  in  different 
geographical  areas\  The  map  is  generalized  and  does  not  show 
the  range  of  soil  types  that  occur  even  within  the  areas 
indicated . 


Table  3.1  is  a  list  of  common  materials,  including  soils, 
with  their  approximate  relative  dielectric  constants.  The 
relative  dielectric  constants  range  from  1  for  air  to  81  for 
water.  The  relative  dielectric  constant  determines  the  speed 
of  electromagnetic  energy  in  that  medium.  That  speed  is 
inversely  proportional  to  the  square  root  of  the  dielectric 
constant  for  the  medium.  This  relationship  is  given  by: 


Speed  in  medium  =  c  / 


(3.19) 


wrier e  c  is 


the  speed  of  light  in  a  vacuum. 


TABLE  3-1.  TYPICAL  RELATIVE  DIELECTRIC  CONSTANTS  FOR  SELECTED 
MATERIALS . 


Material 

Air 

1 

Pure  Water 

81 

Seawater 

SI 

Freshwater  ice 

4 

Seawater  ice 

6 

Snow  ( firm) 

1  .  4 

Sand  (dry) 

5 

Sand  (saturated) 

30 

Clay  (saturated) 

1  0 

Granite  (dry) 

5 

Granite  (wet) 

7 

Limestone  (dry) 

7 

Limestone  (wet) 

8 

Shale  (wet) 

7 

Sandstone  (wet) 

6 

Soil 

sandy  dry 

2-4 

sandy  wet 

20-25 

loamy  dry 

2-6 

loamy  wet 

1  5-20 

clay  dry 

2-6 

clay  we t 

1  0-20 

Permafrost 

6-13 

Strong  Conor 

ete  -  dry 

5-9 

soaked 

20  h  r  s 

10-15 

Cracked  concrete  -  dry 

4-5 

soaked 

20  hrs 

1  3-20 

Asphal  t 

12-16 

Asphal 


For  example,  assume  that  an  electromagnetic  wave  is 
propagating  through  a  soil  with  a  relative  dielectric  constant 
cd  16.  Since  the  speed  of  light  in  a  vacuum  is  approximately 
3x13®  meters/second,  the  speed  of  the  wave  in  the  soil  is 
3x1  0®/ (  1  6  ,  or  7.5x10"^  meters/second. 

Since  the  speed  of  electromagnetic  energy  in  a  medium  is 
i.nversely  proportional  to  the  square  root  of  the  relative 
dielectric  constant  of  that  medium,  the  wavelength  of  the 
energy  is  also  inversely  proportional  to  this  quantity.  This 
derives  from  the  relationship  between  frequency  (f),  wavelengnn 
(A1,  and  speed( v  )  of  electromagnetic  energy: 

f\  =  u  (3.20) 

The  variation  in  speed  and  wavelength  with  respect  to  the 
relative  dielectric  constant  of  the  medium  i.mplies  that  range 
accuracy  and  resolution  are  medium-dependent.  The  impact  of 
this  observation  on  the  operation  of  a  G?R  is  discussed  in  the 
next  section. 

In  order  to  evaluate  GPR  per formance • in  a  given  soil,  one 
must  know  all  of  the  important  soil  parameters  in  addition  to 
the  (iielectric  constant.  Unfortunately,  all  of  these 
parameters,  including  the  dielectric  constant,  are  frequency- 
dependent.  To  measure  the  parameters  accurately  with  respect 
to  frequency,  a  soil  sample  must  be  taken  and  measured. 

Typical  test  equipment  for  making  soil  parameter 
measurements  at  radio  frequencies  ( RF )  includes  a  test  cell  and 
a  microwave  network  analyzer.  The  test  cell  is  usually  a 
section  of  RF  transmission  line  that  is  filled  with  the  soil 
sample.  The  propagation  characteristics  of  the  RF  energy  in 
the  test  cell  filled  with  the  sample  are  measured  with  the 
network  analyzer.  Since  the  propagation  char acter ist ics  of  the 
unfilled  test  cell  are  known,  the  soil  parameters  can  be 
determi.ned  by  comparing  the  results  of  the  measurements  with 
the  filled  and  unfilled  test  cells. 

Appendix  A  contains  plots  of  the  measured  parameters  of 
various  soil  samples.  All  of  the  measurements  were  performed 
by  Dr.  Glenn  Smith  of  the  Georgia  Institute  of  Technology, 
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School  of  Electrical  Engineering.  Dr.  Smith  used  a  Kevlett- 
Packard  8438  network  analyser  and  a  calibrated,  open-circuir 
test  cell. 

The  first  group  of  measurements  were  made  with  a  sample 
clay  soil  taken  near  the  Gulf  Applied  Research  facility  in 
Marietta,  Georgia.  The  moisture  content  of  the  clay  was  var 
and  a  set  of  measurement  curves  was  generated.  The  curves 
indicate  the  variations  in  relative  dielectric  constant  and 
attenuation  with  respect  to  moisture  content. 

As  the  clay  moisture  content  increases,  so  does  the 
relative  dielectric  constant  and  the  attenuation  of  the  ener 
as  it  travels  through  a  given  depth  of  soil.  One  should  not 
that,  in  the  frequency  range  of  the  measurements,  the 
attenuation  for  a  given  moisture  content  increases  with 
frequency.  This  is  true  for  typical  soils.  Also,  note  that 
the  attenuation-per-meter  curves  are  for  one-way  travel  thro 
t.ne  soil.  Thus,  if  a  target  is  located  one  meter  below  the 
surface,  this  one  meter  represents  the  one-way  travel  distan 
However,  the  radar  energy  actually  travels  two  meters  as  it 
travels  from  the  transmitter  antenna  to  the  target  and  back 
the  receiver  antenna. 

The  next  group  of  measurements  was  made  with  a  sample  c 
sand.  At  low  to  moderate  moisture  contents,  the  sand  has  a 
relative  dielectric  constant  and  low  attenuation.  The 
characteristics  of  the  sand  are  considerably  different  than 
those  of  the  clay  discussed  above. 

The  last  group  of  measurements  in  Appendix  A  was  made  f 
several  samples  of  soil  taken  in  Japan.  The  locations  where 
the  soil  samples  were  taken  are  unknown.  No  attempt  was  mad 
to  vary  the  moisture  content  of  the  soil  samples. 
.Nevertheless,  the  'vrariation  in  the  soil  parameters  on  a  sam.p 
to- sample  basis  is  remarkable. 


3.1.3  EFFECTS  OF  SOIL  PARAMETERS  ON  GPR  SYSTEM  PERFORMANCE 

The  properties  of  the  soil  greatly  affect  the  performan 
of  a  GPR.  The  relative  dielectric  constant  of  the  soil 
determines  the  maximum  dept.h  resolution  possible  for  a  give.n 
ba.ndwidth.  Gepth  resolution  is  defined  as  the  capability  of 
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distinguishing  two  targets  that  are  closely  spaced  in  range 
(time).  The  resolution  of  any  radar  system  is  related  to  its 
R?  bandwidth,  or  equivalently,  to  its  RF  pulse  width.  As  the 
relative  dielectric  constant  of  the  soil  increases,  the 
resolution  of  a  GPR  with  a  given  R?  bandwidth  improves. 

The  resolution  capabilities  of  a  GPR  are  easiest  to 
visualize  if  a  pulse-based,  time-domain  system  is  considered. 
Nonetheless,  a  frequency-domain  GPR,  such  as  one  that  uses 
continuous  frequency  modulation,  has  resolution  capabilities 
based  o.n  its  RF  bandwidth  as  well. 

Resolution  for  a  pulse  radar  is  usually  defined  as  half  of 
the  pulse  length  in  the  medium.  A  GPR  that  generates  a  one- 
nanosecond  pulse  in  air  has  a  pulse  length  of  0.3  meter.  This 
gives  the  GPR  a  resolution  of  0.3/2  =  0.15  meter  in  air.  In  a 
soil  with  a  dielectric  constant  of  16,  the  resolution  improves 
to  0.15/(16)''^  =  0.0375  meter. 

Not  only  will  the  relative  dielectric  constant  affect  the 
GPR  resolution,  this  constant  must  be  known  exactly  in  order  to 
measure  the  target  depth  accurately.  Although  the  time 
difference  between  the  surface  return  and  a  subsurface  target 
return  can  be  accurately  measured  with  the  GPR,  a  knowledge  of 
the  speed  of  the  RF  energy  is  required  to  determine  the 
corresponding  depth  of  the  target.  Without  a  knowledge  of  the 
soil  relative  dielectric  constant,  the  speed  of  the  RF  energy 
in  the  soil  is  unknown. 

Fortunately,  in  many  GPR  applications,  precise  depth 
accuracy  is  not  required.  An  experienced  operator  can  often 
estimate  the  relative  dielectric  constant  of  the  soil  well 
enough  for  the  purposes  of  most  field  operations.  In  practice, 
it  is  usually  possibl  i  .to  obtain  soil  samples  for  the  purpose 
of  measuring  the  dielectric  constant  in  regions  near  the 
surface.  The  dielectric  constant  of  the  surface  sample  can 
sometimes  be  used  as  a  reasonable  approximation  of  the 
dielectric  constant  of  deeper  soil. 

The  soil  attenuation  characteristics  can  affect  the  GPR 
resolution  and  depth  penetration.  The  attenuation  curves 
usually  show  increasing  Icss-per-meter  with  increasing 
frequency.  This  means  that  soil  acts  like  a  low-pass  filter  for 
the  RF  energy  passing  through  it.  Removing  some  of  the  high 
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frequency  energy  from  an  RF  pulse  has  the  effect  of  slowing 
rise  time  of  the  pulse  and  widening  it  as  well.  The  degree 
which  this  occurs  depends  cn  the  severity  of  the ' at tenuac i on 
curves  for  a  particular  soil.  As  a  GPR  pulse  widens  on  its 
passage  through  lossy  soil,  the  effective  resolution  of  the 
decreases.  Thus,  one  finds  that  most  GPR  systems  designed  f 
deeper  penetration  (and  corresponding  poorer  resolution)  use 
wider  pulses  than  those  designed  for  shallow  penetration  and 
higher  resolution. 

The  depth  penetration  capability  of  a  given  GPR  is 
primarily  dependent  on  the  attenuation  characteristics  of  th 
soil.  As  RF  energy  passes  through  the  soil,  it  encounters 
targets  that  reflect  a  portion  of  the  energy.  The  energy  is 
attenuated  as  it  travels  through  the  soil  toward  the  target 
well  as  back  to  the  surface  from  the  target.  The  amount  of 
attenuation  due  to  the  soil  parameters  is  based  on  the 
attenuation  curves  fcr  that  particular  soil.  Unfortunately, 
the  attenuation  is  not  linear  with  depth  and  it  is  usually 
measured  in  decibels-per-meter  (d3/m). 

For  example,  assume  that  an  RF  pulse  passes  through  soi 
with  a  23  d3/m  attenuation  characteristic  and  that  this 
attenuation  is  uniform  with  respect  to  frequency.  Thus,  e.ne 
traveling  through  one  meter  of  soil  is  attenuated  23  d3 .  Th 
means  that  the  pulse  power  level  after  one  meter  is  3.31  ti.m 
the  power  level  as  the  pulse  enters  t.he  soil.  If  the  pulse 
travels  two  meters,  it  is  attenuated  A0  d3  and  the 
corresponding  power  level  at  that  poi.nt  is  3.303'!  ti.mes  the 
original  power  level. 


3.2  THE  TARGET  IDENTIFICATION  PROBLEM 


The  main 


purpose  ot 

-nt-.  , 


a  GPR  is  to  detect  subsurface  burie 
oblects,  or  targets.  These  targets  can  be  utility  pipes, 
tun.nels,  cavities,  buried  cables,  etc.  Signal  returns  from 
GPR  are  typically  processed  and  displayed  to  identify  the 
returns  from  desired  targets.  However,  the  signals  received 
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The  RF  energy  transmitted  from  the  antenna  system  can 
encounter  many  interfaces  and  objects.  When  the  propagating  ■ 
energy  impinges  on  an  interface  or  object,  a  portion  of  the 
energy  is  reflected  back  toward  the  source.  The  rest  of  the 
energy  either  scatters  in  other  directions  or  passes  through 
the  object  or  interface.  The  ajnount  of  energy  reflected  back 
toward  the  source  is  determined  by  the  size  and  shape  of  the 
interface  as  well  as  by  the  difference  between  the  relative 
dielectric  constants  of  the  materials  at  the  i.nterface. 
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target  returns.  Clutter  also  refers  to  any  undesired  signals 
that  appear  in  the  receiver  that  are  coherent  with  the 
transmitted  signal.  This  can  include  connector  reflections  i.n 
the  R'^  circuitry,  multiple  bounces  in  the  antenna,  and  noise 
ccupli.ng  from  the  tra.nsmitter  to  the  receiver  over  the  power 
supp ly  1  -nes . 

Clutter  has  the  tendency  to  mask,  or  interfere  with,  the 
returns  f-om  desired  targets.  With  most  broadband  G?R  systems, 
the  internally  generated  or  ''system”  clutter  (not  related  to 
false  targets)  limits  the  detectabi 1 i ty  of  a  desired  target. 

One  refers  to  such  a  system  as  being  "clutter-limited, ”  as 
opposed  to  a  "noise-limited”  system,  in  which  target 
detectability  is  limited  by  the  thermal  noise  level. 

If  the  target  signal  level  falls  below  the  clutter  _evei, 
It  becomes  more  difficult,  if  not  impossible,  to  reliably 
detect  the  target.  Fortunately,  certain  signal  processing 
techniques  can  reduce  the  effects  of  the  clutter.  Such 
techniques  for  G?R  systems  can  be  quite  diverse.  Most  of  them 
are  developed  from  techniques  used  with  conventional  radar 
systems.  Examples  of  these  techniques  are  clutter  subtraction, 
time-domain  synt.het  ic 'aperture  processing^,  various  types  of 
digital  filter  processing^,  and  probabilistic  target 


:r iminat i cn“ .  *hese  techniques  vary  in  implementation, 
jlexity.  and  speed.  The  selection  of  the  specific  technique 
;se  in  a  given  case  is  dependent  upon  the  application  of 


21 .  the  task  of  identifying  a  desired  target  becomes 
inguis.ning  its  signal  return  from  the  returns  of 
argets,  system  clutter,  and  system  noise. 


Furthermore,  the  need  to  satisfy  this  requirement  is 
independent  of  the  type  of  RF  modulation  used  by  the  GPR .  In 
essence,  this  is  a  key  performance  criterion  of  any  GPR  system 


3.3  DESIGNING  THE  GPR  SYSTEM 

This  section  presents  some  of  the  design  considerations 
that  Influence  the  development  of  a  GPR  system.  The  design 
goals  are  usually  based  on  the  required  penetration  depth, 
target  resolution,  RF  power  radiation  limits,  the  GPR  system 
size,  and  the  system  cost.  These  requirements  must  always  be 
referenced  to  some  specific  soil  parameters.  This  section 
focuses  on  the  relationship  between  GPR  resolution  and 
bandwidth,  as  well  as  on  RF  power  requirements,  general  antenr. 
design  considerations,  and  the  elimination  of  system-generated 
clutter . 


3,3.1  GPR  BANDWIDTH  AND  RESOLUTION 

The  penetration  depth  and  resolution  specifications 
determine  the  RF  bandwidth  required  for  a  GPR.  The  bandwidth 
may  be  developed  using  a  direct  time-domain  pulse;  or  by  using 
a  frequency-domain  system,  such  as  FM-CV  or  stepped-FM.  The 
frequency-domain  systems  provide  the  most  flexibility  in 
tailoring  the  RF  bandwidth,  since  generating  optimum  time- 
domain  pulses  directly  is  difficult.  With  frequency-domain 
systems,  the  frequency  components  can  be  individually  phase-  c 
amplitude-modulated  to  achieve  beneficial  system 
characteristics.  This  is  not  as  easy  to  accomplish  with  pulse 
r  adar  s . 

Due  to  the  soi i' attenuation  characteristics,  good  depth 
pe.netration  requires  the  use  of  low-frequency  RF  energy.  In 
contrast,  high  resolution  requires  high-frequency  energy. 
Figure  3.2  plots  two  rectangular  pulses.  The  rectangular 
pulses  each  have  normalized  amplitudes  of  one  volt.  Figure 
3.3a  is  the  frequency  spectrum  of  the  two-nanosecond  pulse  frc 
Figure  3.2a  and  Figure  3.3b  is  the  spectrum  of  the  higher- 
resolution  one-nanosecond  pulse.  Note  that  the  bandwidths  of 
the  pulses  are  inversely  proportional  to  their  pulse  widths. 
The  rectangular  pulse  shapes  were  selected  purely  as  an  examp  1 
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and  do  not  necessarily  represent  optimum  pulse  shapes  for  a  G?R 
design . 


Since  the  soil  acts  as  a  lowpass  filter  to  electromagnetic 
energy,  a  compromise  will  exist  in  the  GPR  design  between 
resolution  and  depth  penetration.  In  fact,  with  increasing 
depth  penetration,  a  GPR  with  a  given  RF  bandwidth  will  have 
better  resolution  at  shallow  depths  than  it  will  at  deeper 
depths.  For  this  reason,  many  GPR  systems  utilize  different  RF 
bandwidths  and  antennas  depending  on  their  application. 


3.3.2  GPR  TRANSMITTER  POWER  CONSIDERATIONS 

Once  the  GPR  RF  frequency  range  has  been  selected,  the 
transmitter  power  requirement  must  be  identified.  The  power 
requirement  is  a  function  of  the  desired  depth  penetration,  the 
soil  attenuation,  the  target  radar  cross  section,  the 
efficiency  of  the  antennas,  the  receiver  noise  level,  and  the 
applicable  regulations  regarding  RF  energy  emissions.  The  G? 
designer  must  estimate  the  expected  signal  level  from  a  burie 
target  given  a  reference  power  level.  If  the  estimate 
indicates  that  the  signal  level  will  fall  below  the  minimum 
detectable  receiver  signal,  then  the  power  transmitted  must  be 
increased.  The  limits  on  the  maximum  radiated  power  are 
usually  regulated  by  governmental  agencies  and  are  based  on 
safety  and  RF  interference  requirements. 

Unfortunately,  merely  increasing  the  GPR  transmitter  power 
does  not  necessarily  mean  that  signals  from  deeper  subsurface 
targets  will  be  identifiable  at  the  receiver.  The  system¬ 
generated  clutter  usually  proves  to  be  the  limiting  factor. 

The  clutter  will  increase  linearly  with  increasing  transmitter 
power,  as  will  the  desired  target  returns.  Therefore,  if  a 
target  signal  is  masked  by  an  internally-generated  clutter 
signal  when  a  low-power  transmitter  is  used,  increasing  its 
power  will  increase  the  clutter  signal  as  well,  so  that  the 
target  signal  will  still  be  masked. 

The  internally-generated  clutter  often  appears  in  fixed 
ti.me  positions  within  the  receiver  output  signal  range  window. 
Thus,  this  type  of  clutter  does  not  always  interfere  with  the 
target  return  signal.  If  the  receiver  window  is  positioned 
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optimally  in  time,  or  the  internal  clutter  is  moved  away  from 
the  receiver  time  window,  then  the  GPR  will  usually  benefit 
from  increased  transmitter  power.  The  optimal  positioning  of 
the  internal  clutter  is  not  always  possible,  however. 

The  reader  should  note  that  soils  have  one-way  power 
attenuation  characteristics  that  are  on  the  order  of  tens  of 
decibels  per  meter.  To  increase  the  transmitter  output  power 
10  dB ,  the  output  signal  voltage  must  increase  by  a  factor  of 
3.16.  To  increase  it  by  20  d3 ,  the  output  voltage  must  be  10 
times  higher.  Large  increases  in  transmitter  power  can  place 
heavy  demands  on  the  transmitter  design,  especially  for  short- 
pulse  systems.  This  is  because  higher-power  transmitters 
usually  require  electronic  components  that  can  withstand  large 
power  supply  voltages  and  currents.  Such  components  typically 
cannot  produce  pulses  as  narrow  as  the  ones  produced  by  lower- 
power  devices. 


3.3.3  GPR  A3rrEN?fA  CHARACTERISTICS 

One  of  the  key  elements  in  a  GPR  design  is  the  seleotion 
of  a  proper  antenna.  A  GPR  can  use  a  single  antenna  for  both 
transmitting  and  receiving  (monostatic  configuration),  or  one 
antenna  for  transmitting  and  a  separate  antenna  for  receiving 
(bistatic  configuration).  In  any  case,  the  antenna(s)  used 
must  provide  a  good  RF  match  to  the  transmitter  and  receiver 
aind  should  be  as  efficient  as  possible.  In  addition,  the 
antenna  size  and  radiation  pattern  may  be  important 
considerations  for  particular  GPR  applications. 

Typical  GPR  applications  require  relatively  broadband  RF 
systems.  The  antenna  must  have  the  capability  of  radiating  as 
much  of  the  transmitter  output  energy  as  possible.  Thus,  the 
RF  bandwidth  of  the  GPR  must  match  the  radiation  bandwidth  of 
the  antenna.  Cr ,  equivalently,  the  antenna  must  be  designed  c 
match  the  RF  bandwidth  of  the  transmitter . 

Broadband  antennas  act  much  like  highpass  or  bandpass 
filters.  They  cannot  radiate  energy  below  their  low-frequency 
cutoff.  Therefore,  if  the  RF  frequency  band  of  the  GPR  is  sue 
that  it  extends  below  the  low-frequency  cutoff  of  the  antenna, 
that  portion  of  the  energy  below  the  cutoff  will  not  be 


radiated.  A  similar  situation  exists  at  the  antenna  upper 
cutoff  frequency. 

One  desires  that  the  transmitting  antenna  absorb  the 
energy  that  it  cannot  radiate.  Otherwise,  the  antenna  will 
present  a  high  reflection  coefficient  for  the  frequencies 
outside  its  efficient-radiation  band.  Any  reflected  energy 
from  the  antenna  input  will  propagate  back  toward  the  source. 
This  can  be  disastrous,  especially  for  a  monostatic  system, 
because  the  reflected  energy  will  contribute  to  the  internal 
system  clutter  and  can  mask  weak  targets  if  it  should  appear  ir. 
the  receiver  window. 

An  additional  requirement  for  the  antenna  is  that  it  be 
non -dispersive .  This  implies  that  the  group  delay  for  the 
frequencies  in  the  antenna  passband  is  constant.  If  a  pulse 
system  uses  an  antenna  with  a  constant  group  delay  and  a 
bandwidth  that  is  wide  enough  for  the  RF  frequency  range,  the 
radiated  and  received  pulses  will  not  ring  or  stretch 
appreciably.  This  will  maintain  a  minimum  pulse  resolution 
cell.  These  benefits  are  also  desirable  for  CW  systems. 

However,  one  can  measure  the  antenna  characteristics  on  a 
discrete-frequency  basis  and  correct  for  certain  antenna 
distortions,  whether  in  ajnplitude  or  phase.  Unfortunately, 
these  corrections  are  not  always  easily  performed.  The  best 
approach  to  minimizing  antenna-related  distortion  problems  is 
with  careful  antenna  design,  thereby  eliminating  as  much 
waveform  pre-processing  or  post-processing  as  possible. 

The  antenna  propagation  and  radiation  characteristics  are 
related  to  the  length  and  shape  of  the  antenna.  For  a  given 
antenna  with  a  known  radiation  bandwidth,  doubling  its  length 
will  reduce  the  low-frequency  cutoff  by  a  factor  of  two.  Thus, 
antennas  for  GPR  systems  that  are  designed  f-'r  deep  penetrai.ion 
are  usually  larger  than  those  used  with  high-reso 1 ut ion  GPR 
systems.  This  is  because  soils  typically  attenuate  high  RF 
frequencies  more  than  low  RF  frequencies.  Therefore,  deep- 
penetrating  GPR  systems  have  transmitter  bandwidths 
concentrated  at  lower  frequencies.  Such  concentration  at  lower 
frequencies  creates  larger  pulse  widths  in  time-dcmain  systems. 
In  contrast,  the  high-resolution  GPR  systems  require  a 
relatively  large  RF  bandwidths,  which  is  equivalent  to  a  narrow 
pulse. 
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The  antenna  size  also  affects  its  radiation  pattern.  Many 
antennas  have  radiation  patterns  that  can  be  approximated  as 
though  originating  from  apertures.  Larger  antennas  approximate 
larger  apertures.  For  a  given  frequency  of  electromagnetic 
energy  being  radiated,  the  radiation  pattern  from  a  small 
aperture  will  be  broader  than  that  from  a  large  aperture. 

.\lso,  for  a  given  aperture  size,  the  radiation  pattern  will  be 
broad  for  low  frequencies  and  narrow  for  high  frequencies. 

This  implies  that  the  antenna  radiation  pattern  will  not  be 
constant  over  the  typical  bandwidths  of  GPR  systems.  This 
phenomenon  will  affect  system  resolution  at  target  angles  off- 
boresight  from  the  antenna. 


3.3.4  REDUCING  INTERNAL  CLUTTER  IN  A  GPR  SYSTEM 

Internal  clutter  in  a  GPR  system  is  minimized  by  careful 
design  and  construction  practices.  A  critical  factor  affecting 
the  clutter  levels  is  the  design  of  the  transmitter  and 
receiver  transmission  systems.  All  R?  components  should  be 
matched  as  well  as  possible  to  the  transmission  lines  (usually 
50  ohm  coaxial  cables)  to  eliminate  reflection  points.  This 
includes  antennas,  transmi tters ,  samplers,  RF  amplifiers, 
mixers,  etc.  Any  reflections  should  be  absorbed  or  moved  in 
time  with  RF  delay  lines  to  keep  them  out  of  the  receiver  ran 
window.  Any  discontinuities  in  the  transmission  lines,  such 
connectors,  will  contribute  to  the  reflections. 

Clutter  may  also  enter  the  transmitter  and  receiver 
systems  through  the  power  supply  lines  for  the  system 
components.  For  example,  a  transmitter  pulser  can  place  high- 
frequency  noise  on  its  power  supplies  and  ground  plane.  If 
this  power  supply  noise  is  not  filtered  or  removed  by  other 
techniques,  it  can  enter  the  receiver  section  through  the 
supply  lines.  This  type  of  noise  can  easily  add  significant 
clutter  to  the  receiver  output. 
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SECTION  4 


MODELING  THE  PIPE  DETECTION  PROBLEM 


The  initial  approach  to  determining  the  optimum  GPR  signal 
source  required  detailed  investigation  of  the  problem  of 
subsurface  pipe  detection.  GAR  performed  an  approximate 
analysis  to  estimate  the  maximum  depth  capability  of  a  GPR  in  a 
known  environment.  The  analytical  model  selected  was  flexible 
e.nough  to  include  effects  of  soil  attenuation,  antenna 
parameters,  surface  reflection  loss,  wave  diffraction  at  the 
surface  with  variable  antenna  height,  and  the  target  radar 
cross  section.  The  model  did  not  include  the  effects  of 
clutter,  whether  generated  internally  or  externally.  GAR 
adjusted  the  model  to  simplify  the  computations.  The 
dominating  factors  controlling  the  signal  propagation  and 
attenuation  were  retained  in  the  analysis.  These  factors 
included  the  antenna  gain-bandwidth  product,  the  soil 
attenuation  characteristics,  the  transmission  coefficient  at 
the  air/soil  interface,  and  the  radar  cross  section  of  the 
pipe  . 


4.1  SUBSURFACE  TARGET  DETECTION  MODEL 

The  model  selected  was  developed  by  the  National  Bureau  of 
Standards  ( N3S )  as  part  of  a  research  project  for  the  U.S. 
Department  of  the  Interior,  Bureau  of  Mines’.  This  section 
describes  that  model  and  the  approximations  used  by  GAR  to 
simplify  the  pipe  target  analysis.  The  basis  for  the  model  is 
the  radar  range  equation.  The  model  estimates  the  received 
signal  power  as  a  function  of  GPR  system  parameters,  soil 
parameters,  and  target  parameters.  The  model  accounts  for  only 
o.re  RF  frequency  at  a  ti.me;  thus,  the  model  should  be  evaluated 
fsr  all  of  the  discrete  frequencies  in  the  RF  frequency  band  cf 
interest.  The  radar  equation  is  given  below: 


PG. 

S  =  - 
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where 


S 


received  signal  power. 


transmitter  output  power, 

gains  of  the  transmit  and  receive  antennas, 
respectively 

target  radar  cross  section, 

wavelength  of  electromagnetic  energy  in  air, 

the  attenuation  factor  of  the  complex  propagation 
constant , 


I 


two-way  power  transmission  at  the  air/soil 
interface,  accounting  for  dielectric  differences, 


D 


dispersion  effect  due  to  the  medium  resulting  in 
pulse  stretch  and  subsequent  amplitude  reduction, 


R  =  range  from  the  radar  antenna  to  the  target, 

Ri  =  the  part  of  the  range  in  the  soil  only. 


C 


d 


factor  accounting  for  the  two-way  refraction  ao 
the  air/soil  interface. 


In  Equation  (4.1),  the  left-most  fraction  is  the  power 
density  at  the  target,  and  the  second  is  the  reflected  power 
density  at  the  recei'ver,  based  on  the  target  radar  cross 
section.  These  power  densities  assume  a  free  space 
environment.  The  rig.ht-most  fraction  is  the  effective  aperture 
of  t.he  receiving  antenna.  The  remaining  terms  account  for  the 
losses  associated  with  the  soil  and  the  air/soil  interface. 

Several  simplifying  assumptions  were  made  for  the  GAR 
analysis  using  the  model  in  Equation  (4.1).  First,  the  antenna 
and  its  phase  center  were  assumed  to  be  on  the  surface  so  that 
R-  =  R-  Second,  the  factors  and  D  were  neglected.  Third,  a 
monostatic  G?R  configuration  was  assumed,  so  G,  =  G, .  Fourth, 
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the  soil  was  assumed  to  be  homogeneous.  The  first  two 
simplifications  were  instituted  because  detailed  information 
concerning  the  antenna  radiation  pattern  for  each  of  the 
frequencies  of  interest  was  required  if  the  factors  and  D, 
as  well  as  the  antenna  phase  center  were  included  in  the 

analysis.  This  information  was  not  available  and  it  was  felt 

that  these  parameters  did  not  have  the  significant  impact  on 
the  analysis  that  the  remaining  parameters  had. 

The  radar  cross  section  of  a  target  is  defined  as  the 
ratio  of  the  power  density  of  the  scattered  wave  at  the 
receiver  to  the  power  density  of  the  incident  wave  at  the 
carget ,  multiplied  by  4tR^,  where  R  is  the  distance  from  the 
source  to  the  target.  For  this  analysis,  a  pipe  target  was 
assumed  to  be  an  infinitely  long  cylinder  and  the  incident 
electromagnetic  field  was  assumed  to  be  polarized  parallel  to 
the  pipe  axis.  The  N3S  radar  cross  section  model  of  an 

infinite  cylinder  assumed  that  the  distance  from  the  source  t 

the  cylinder  was  very  large  with  respect  to  the  wavelength  cf 
the  electromagnetic  energy,  so  that  there  was  plane-wave 
incidence  at  the  cylinder. 

Equation  (4.2)  gives  the  radar  cross  section  for  the 
cylinder  at  plane  wave  incidence: 


cTpyj  =  8R  I  P  I  ^  /  k 


where,  P  is  given  by: 


P 


-  E  (-1)"  (ka) 

n  =  0 


cos  no 
(ka) 


(4.2) 


(4.3) 


Here,  e^=l  for  n  =  0  and  e^  =  2  otherwise,  J„  are  Bessel  functions 
are  Hankel  functions,  ,  and  a  is  the  cylind 

r  adi us  . 


For  regions  where  ka  <  0.05,  the  solution  for  Equation 
(4.3)  may  be  approximated  by  a  Rayleigh  function.  For  region 
where  ka  >  6,  Equation  (4.3)  may  be  approximated  by  a  gecmetr 
optics  solution.  For  the  intermediate  region  of  ka,  the  e.xac 
solution  for  Equation  (4.3)  is  used  for  maximum  accuracy. 
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The  value  of  ka  for  most  GPR  applications  falls  into  the 
intermediate  region.  However,  to  simplify  the  analysis 
considerably,  at  the  risk  of  some  tolerable  inaccuracy,  rhs  GAF, 
approach  utilizes  the  geometric  optics  approximation  for  the 
radar  cross  section  of  the  pipe  using  a  monostatic  radar.  This 
radar  cross  section  is  given  by: 

<7-  Rrra  (4.4) 

j 


4.2  EXAMPLE  CALCULATION  USLNG  THE  NBS  MODEL 

The  following  example  will  illustrate  the  use  of  the  N 
model  for  predicting  the  radar  return  from  a  buried  pipe, 
simplifying  assumptions  outlined  in  Section  4.1  will  apply  for 
this  example.  Although  the  accuracy  of  the  model  is 
compromised  by  the  assumptions,  the  reader  should  keep  in  mind 
that  the  object  of  the  assumptions  was  to  allow  a  quick 
assessment  of  potential  GPR  performance  in  a  given  environment. 
Note  that  the  radar  ra.nge  equation  model  that  is  used  assumes  a 
homogeneous  soil  and  neglects  the  effects  of  clutter.  Also, 
the  assumptions  are  made  that  the  GPR  antenna  system  is 
monostatic  and  the  antenna  (and  its  phase  center)  is  in  contact 
with  the  ground  surface. 

The  key  element  in  the  radar  range  equation  model  is  the 
specification  of  the  antenna  parameters.  GAR  selected  a  TEri 
horn  antenna  for  this  application.  This  particular  antenna  is 
used  by  GAR  for  its  RODAR*-  ground-penetrating  radar  system. 

The  antenna  was  originally  developed  by  R.  Wohlers  at  Calspan, 
and  its  radiation  characteristics  in  air  are  well  documented®. 
The  antenna  has  many  of  the  desirable  properties  that  a  pulse 
GPR  system  requires..  It  has  a  low  voltage  standing  wave  ratio 
( VSWR ) ;  it  is  extremely  broadband;  it  has  good  group  delay 
properties;  and  it  has  low  pulse  stretch. 

The  antenna  gain,  i.n  air.  was  measured  for  one  of  the 
Calspan  TSM  antennas  and  is  given  approximately  by®: 

G  =  4t  (LVA‘)  sin  9^  ,  (4.5) 

where  L  is  the  antenna  length,  9^  is  equal  to  the  angle  betweer 
the  two  antenna  plates  divided  by  two,  and  A  is  the  wavelength 
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of  the  RF  energy.  The  antenna  under  consideration  had  a 
length,  L,  of  one  meter  and  an  angle  between  the  plates  of  11.6 
degrees.  The  effective  aperture  of  the  antenna  was  measured  in 
air^  and  a  plot  of  the  measured  effective  aperture  vs. 
frequency  is  shown  in  Figure  A . i . 

The  performance  of  the  antenna  in  contact  with  the  soil  is 
different  than  its  performance  in  air.  Specifically,  its  gain 
and  effective  aperture  increase  somewhat  when  the  antenna 
contacts  the  ground  due  to  the  larger  dielectric  constant  of 
Che  soil.  Also,  the  soil  interface  is  in  the  near-field  region 
of  the  antenna.  Thus,  it  is  difficult  to  evaluate  the 
transmitted  power  density  at  the  soil  interface.  For  this 
reason,  the  parameters  of  the  antenna  measured  in  air,  rather 
than  those  measured  in  the  ground-contacting  mode,  are  used  in 
the  radar  range  equation  model. 

The  soil  parameters  selected  for  this  example  were 
obtained  from  a  red  clay  soil  found  near  the  GAR  facility  in 
>Iarietta,  Gecrgia.  It  was  assumed  that  the  clay  had  a  moisture 
content  of  20  percent  by  dry  weight.  The  soil  data  are 
presented  in  Appendix  A. 

The  RF  frequency  ratige  for  the  example  was  chosen  to  be 
from  100  megahertz  (MHz)  to  300  MHz.  This  bandwidth  is 
equivalent  to  about  a  five-nanosecond  pulse  width  in  the  time- 
domain.  The  radar  range  equation  was  evaluated  at  100  MHz,  200 
MHz,  and  300  MHz. 

The  subsurface  target  was  selected  to  be  a  12-inch 
diameter  metal  pipe  buried  at  a  depth  of  t.hree  meters  (measured 
from  the  surface  to  the  top  of  the  pipe).  The  pipe  was 
considered  to  be  infinitely  long.  Its  radar  cross  section  was 
assumed  to  be  described  by  the  NBS  model  for  a  cylinder  in  the 
geometric  optics  regidn. 

Equation  ( i  .  6  )  below  is  the  modified  N3S  radar  range 
equation  that  was  used  for  the  GAR  evaluation.  The  equation 
expresses  the  ratio  of  received  power  to  transmitted  power  as  a 
function  of  the  soil  parameters,  a.ntenna  parameters,  and  target 
parameters . 

S/P  =  - ^ - ^  G.  I  e  .  (4.6) 

47rR^  4,tR2 
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horn  antenna  effective  aperture 


where,  G,)/47r  is  the  effective  receiving  aperture  of 

the  antenna,  and  the  factors  and  D  in  Equation  (4.1 )  have 
been  neglected. 

Table  4.1  summarizes  the  soil,  antenna,  and  pipe 
parameters  used  in  Equation  (4.6)  for  the  frequencies  of  100 
MHz,  200  MHz,  and  300  MHz.  The  results  obtained  from  Equation 
(4.6)  for  the  S/P  ratios  reflect  tremendous  signal  attenuation, 
extending  from  -142  dB  at  100  MHz  to  -162  d3  at  300  MHz. 

We  feel,  however,  that  these  calculations  are  somewhat 
pessimistic  because  the  performance  of  the  antenna  near  the 
soil  interface  is  unknown,  but  is  likely  better  than  that 
assumed  .i.n  the  calculations.  Thus,  a  more  exact  analytical 
approximation  of  the  received  target  signal  power  would 
probably  yield  a  more  optimistic  answer  for  the  received  power 
level.  However,  the  results  obtained  based  on  our 
approximations  lend  insight  into  the  performance  requirements 
that  must  be  imposed  on  a  successful  GPR  system. 
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TABLE  4.1 


RESULTS  OF  XUMERICAL  ESTI.MATION  OF  PIPE  TARGET 
DETECTION  FOR  FREQUENCIES  OF  100  MHZ, 

200  MHZ,  AND  300  MHZ. 


FREQUENCY  (MHz) 


Parameter 

1  30 

230 

3  33 

a  ( ) 

1  .44 

1  .  44 

A,  (m^) 

0.01 

0.01 

3.01 

G 

0.014 

0 . 058 

3.133 

a  (Nepers) 

1.15 

1  .  50 

1  .73 

I 

3 . 426 

0 . 458 

3 . 458 

(radians/s) 

6 . 28x1  0® 

1  .  26x1  0’ 

1 .88x1 3 

S/P 

6.82x1  0*’5 

4,55x1  0"'^ 

6 . 16x1  3'" 

S/P  (d3) 

-142 

-153 

C\J 

1 
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SECTION  5 


DESIGN  OF  THE  SCALE-MODEL  TEST  TANK 

A  soil  test  tank  was  required  to  test  the  various  GPR 
signal  source  types.  The  test  tank  was  designed  to  contain 
various  pipe  targets,  metal  plate  targets,  and  electr omagnet ic 
field  probes  in  order  to  test  candidate  GPR  systems.  The  GPR 
system  requirements  provided  by  NCSL  included  the  spec i f icat i on 
zo  detect  and  resolve  plastic  and  metal  pipes  buried  at  depths 
of  up  to  23  feet.  Care  was  taken  to  ensure  that  the  dielectric 
medium  in  the  test  tank  was  homogeneous  to  obtain  repeatable 
measurements  and  reduce  the  effects  of  random  external  clutter. 

GAR  decided  to  evaluate  the  GPR  sources  using  a  higher- 
frequency  scale  model  than  the  operating  frequency  of  a  full- 
scale  GPR.  Such  a  scale  model  is  much  easier  to  use  while 
making  laboratory  measurements  than  is  a  full-scale  system 
because  of  its  reduced  size.  However,  the  scale  model  provides 
measurement  results  that  apply  directly  to  the  full-scale 
system . 


5.1  DESCRIPTION  OF  THE  GPR  SCALE-MODEL  CHARACTERISTICS 

To  create  a  high-frequency  scale  model  for  a  GPR,  the 
radar  environment,  as  well  as  the  radar  system  parameters,  muse 
be  adjusted.  All  spatial  dimensions  are  reduced  by  the  scaling 
factor,  resulting  in  a  much  smaller  testing  environment  and  a 
radar  system  that  is  more  manageable.  The  target  dimensions 
and  burial  depths  are  also  reduced  by  the  scale  factor. 

Several  types  of  ■scale  models  are  outlined  in  Antennas  in 
Matter  by  King  and  Smith^.  The  one  selected  by  GAR  permits  the 
use  of  the  same  relative  dielectric  constants  in  the  model  that 
are  found  in  the  full-scale  environment.  This  means  that  air 
can  be  used  as  a  dielectric  in  the  model.  Otherwise,  if  the 
dielectric  constants  were  scaled,  air  could  not  be  used  in  the 
model  environment,  since  its  dielectric  constant  would  have  to 
be  scaled.  The  dielectric  parameter  that  is  scaled  i.n  the  GAR 
model  is  t.he  soil  conductivity,  and  hence,  its  attenuation. 
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For  example,  assume  that  a  full-scale  GPR  system  is  to  be 
modeled  using  a  scale  factor  of  5.  Let  the  RF  bandwidth  of  th 
full-scale  GPR  extend  from  100  MHz  to  300  MHz.  Furthermore, 
assume  that  the  full-scale  soil  parameters  include  a  relative 
dielectric  constant  of  16  and  an  attenuation  characteristic  of 
10  dB/m.  The  scale-model  RF  bandwidth  would  then  extend  from 
500  MHz  to  1500  MHz.  The  model-soil  dielectric  constant  would 
remain  at  16,  but  its  attenuation  characteristic  would  increas 
to  50  dB  per  meter . 

GAR  chose  the  scale-model  soil  tank  dimensions  to  be  six 
feei  wide  by  eight  feet  long  by  four  feec  high.  This 
corresponds  to  full  scale  dimensions  of  30  feet  by  40  feet  by 
20  feet.  The  targets  were  appropriately  placed  in  the  model 
soil  to  avoid  boundary  effects.  They  were  also  placed  in  such 
a  manner  as  to  avoid  any  interference  with  each  other  . 

The  scale  factor  selected  for  the  model  was  five.  This 
scale  factor  allowed  GAR  to  use  its  one-nanosecond  monocycle 
short-pulse  radar  to  model  a  full-scale  GPR  system  with  a  five 
nanosecond  pulse.  A  five-nanosecond  pulse  is  of  the  duration 
that  might  be  used  for  a  GPR  system  designed  for  deep 
penetration  and  medium  resolution. 

GAR  investigated  two  approaches  to  building  a  scale-model 
environment.  One  approac.h  used  a  soil  that  had  the  required 
parameters  of  the  model  dielectric.  The  second  approach  used 
liquid  dielectric.  The  liquid  dielectric  used  an  emulsion  of 
oil,  water,  and  salt  to  achieve  the  required  properties  of  the 
model  dielectric.  Cr .  Glenn  Smith  of  the  Georgia  Institute  of 
Technology  has  previously  investigated  and  used  such  liquid 
dielectrics  as  soil  models. 


S2  THE  LIQUID-DIELECTRIC  SOIL  MODEL 

The  liquid-dielectric  model  had  three  key  advantages  over 
a  soil  dielectric.  First,  the  liquid  tank  has  the  potential  c 
being  more  homogeneous  than  a  soil  dielectric.  Second,  the 
targets  and  their  positions  can  be  easily  changed.  Third,  the 
effects  of  the  liquid  tank  boundaries  can  be  easily  eliminated 
in  any  target  measurements  after  a  "baseline"  measurement  is 
made  of  a  target-free  tank.  In  general,  eliminati.ng  the 


35 


boundary  effects  in  a  soil  dielec.  tank  is  not  possible 
because  the  removal  and  insertion  of  the  targets  would  disturb 
the  soil,  which  could  never  be  replaced  in  exactly  the  same 
configuration.  Furthermore,  a  controlled  clutter  environment 
could  easily  be  placed  in  the  liquid-dielectric  tank  to  examine 
the  effects  of  external  clutter  on  the  GPR  systems. 

The  emulsion  created  for  use  in  the  liquid  tank  could  be 
made  to  model  selected  soil  dielectrics.  Dr.  Smith  uses  a 
high-grade,  light  mineral  oil,  water,  and  salt  as  the  emulsion 
ingredients.  He  has  investigated  many  different  combinations 
of  emulsion  ingredients  and  their  ratios.  The  dielectric 
properties  of  these  emulsions  were  measured  and  used  to  create 
sets  of  curves  that  allow  the  specification  of  a  particular 
emulsion  mix  ratio  to  achieve  a  given  liquid  dielectric 
constant  and  loss  characteristic. 

The  primary  drawback  in  using  a  liquid  tank  with  a  large 
vclume  was  the  need  to  mix  and  maintain  a  homogeneous  emulsion 
of  the  liquid  dielectric.  tank  with  dimensions  of  six  feet 
by  eight  feet  by  four  feet  holds  over  1400  gallons  of  liquid. 

To  correctly  mix  and  maintain  the  emulsion,  a  shearing  gear 
pump  would  be  required  that  must  operate  at  a  high  flow  rate. 
The  emulsion  must  be  circulated  very  quickly  during  mixing  to 
ensure  a  homogeneous  mix.  Furthermore,  the  emulsion  must  be 
re-circulated  and  mixed  prior  to  each  measurement  to  ensure  its 
homogeneity  and  the  accuracy  of  its  dielectric  properties. 

At  the  time  the  GAR  model-soil  parameters  were  specified. 
Dr.  Smith  had  not  built  a  liquid-dielectric  tank  of  such  large 
dimensions.  The  cost  and  logistics  of  fabricating  a  large 
tank,  purchasing  the  necessary  pumps,  and  constructing  a 
working  model  tank  appeared  to  be  quite  formidable.  In 
addition,  there  was  a  significant  technical  risk  in  making  a 
large  tank  work  properly  as  a  model.  GAR  felt  that  the 
technical  risks  and  t.he  costs  associated  with  t.he  liquid  tank 
precluded  its  use  in  the  measurement  p'^ogram.  GAR  elected  to 
construct  a  soil  tank  instead. 
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5.3  THE  SOIL-DIELECTRIC  MODEL 

A  model  tank  containing  soil  is  much  easier  and  less 
costly  to  construct  and  maintain  than  a  corresponding  liquid 
tank.  However,  any  change  in  target  placement  once  objects  are 
buried  is  much  more  difficult  with  the  soil  tank.  In  addition, 
the  effects  of  any  inhomogeneities  and  of  the  tank  boundaries 
are  difficult  to  characterize  and  their  effects  on  the  target 
measurement  are  very  difficult  to  remove. 

In  order  to  build  the  soil  tank,  GAR  had  to  locate  a 
source  of  homogeneous  soil  wich  electromagnetic  properties 
suitable  for  use  as  a  model  soil.  Thus,  a  local  brick 
manufacturer  was  contacted  ci.nd  arrangements  were  made  to  test 
the  properties  of  the  clay  soils  that  they  used  for  the  bricks. 
The  brick  manufacturer  used  two  types  of  clays  in  the  bricks, 
which  were  identified  as  a  clay  shale  and  a  clay  schist. 
Appendix  B  contains  the  measured  characteristics  of*  these 
clays.  The  GPR  source  evaluation  and  .measure.ment  program 
required  the  use  of  at  least  two  soil  dielectrics.  This 
allowed  evaluation  of  the  performances  of  the  GPR  systems  with 
different  soil  types..  Therefore,  the  two  clays  from  the  brick 
manufacturer  would  satisfy  this  requirement. 

The  properties  of  the  clay  shale  approximated  those  of 
actual  clay  with  20  percent  moisture  as  given  in  Appendix  A. 

The  shale  had  an  attenuation  characteristic  of  50  dB/m  at  1 
GHz.  This  attenuation  scaled  to  10  d3/m  at  200  MHz.  The 
schist  had  an  attenuation  characteristic  of  53  d3/m  at  1  GHz, 
which  scaled  to  6.5  d3/m  at  200  MHz. 

GAR  decided  to  use  the  clay  shale  as  the  first  soil  in  the 
soil-model  tank.  Sixteen  targets  were  buried  in  a  six-foot  by 
eight-foot  by  four-foot  high  soil  test  tank.  Thirteen  pipes 
with  diameters  ranging  from  0.88  inch  to  6.63  inches  were  used. 
The  pipe  diameters  increased  with  burial  depth.  Each  pipe 
target  was  two-feet  seven-inches  long.  For  each  metal  pipe 
there  was  a  corresponding  PVC  plastic  pipe  of  similar  diameter. 
The  one  exception  was  the  6.63  inch  diameter  metal  pipe  that 
was  five-feet  six-inches  long  and  buried  at  a  three-foot  depth. 
Mo  plastic  pipe  of  a  comparable  size  was  buried  in  the  model 
tank . 
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Figures  5.1  and  5.2  are  sketches  of  a  top  view  and  a  side 
view,  respectively,  of  the  soil-model  tank.  The  sketches  show 
the  locations  of  the  pipe  targets.  In  general,  one  side  of  the 
tank  was  reserved  for  metal  pipes  and  the  other  was  reserved 
for  plastic  pipes.  The  exception  to  this  was  the  6.63  inch 
diameter  pipe  which  extended  across  the  entire  tank.  In 
addition  to  the  13  pipes,  three  metal  plate  targets  of  various 
sizes  and  burial  depths  were  also  included  in  the  tank.  Each 
of  these  plates  had  a  half-loop  H-field  probe  mounted  on  it. 

The  probes  had  a  0.25  inch  radius  and  used  the  plates  as  an 
image  plane.  The  probes  were  included  as  sensors  to  detscc  any 
energy  reaching  them.  The  reader  should  note  that  energy 
reaching  the  loop  probes  was  only  subjected  to  one-way 
attenuation  through  the  soil. 

The  targets  in  the  soil  tank  were  arranged  to  minimize 
potential  electromagnetic  interference  between  each  other  and 
the  tank  boundaries.  The  pipe  diameters  and  burial  depths  were 
largest  at  the  center  of  the  tank.  The  shallow,  smaller  pipes 
were  located  at  the  tank  edges. 

During  construction ,  care  was  taken  to  maintain  soil 
homogeneity  in  the  tank.  This  was  accomplished  by  tamping  the 
clay  every  six  inches  as  the  tank  was  filled.  The  tank  was 
covered  with  a  thin  sheet  of  polyethylene  plastic  to  minimize 
the  evaporation  of  moisture  from  the  model  soil.  If  the 
moisture  had  evaporated  from  the  soil  tank,  the  electromagnetic 
soil  properties  would  have  been  altered. 

After  the  measurements  were  performed  with  the  clay  shale, 
the  shale  was  removed  and  replaced  with  the  clay  schist. 

Target  and  probe  locations  were  the  same  as  they  were  with  the 
clay-shale  tank. 


38 


Skel.c^h  ol'  l.lie  L»>i»  view  of  l.he  GAR  bo  11  model  tank 
showing  burled  target  locations. 


tank  allowing  burled  target  locatlonu. 


SECTION  6 


DESCRIPTION  OF  THE  GPR  SYSTEMS  TESTED 

I 

I 

GAR  tested  four  GPR  systems  representing  time-domain  and 
[  frequency-domain  technologies.  Two  of  the  systems  used  short- 

pulse,  time-domain  techniques;  and  the  other  two  systems  used 
frequency-modulation  techniques.  The  basic  operations  of  the 
two  time-domain  radar  systems  were  quite  similar,  and  the 
operations  of  the  two  frequency-domain  radars  were  also  very 
similar.  However,  the  time-domain  radar  systems  functioned 
quite  differently  in  concept  from  the  frequency-domain  radars. 

One  of  the  time-domain  radar  systems  tested  was  built  by 
GAR  for  the  RODAR"”^  highway-surveying  short-pulse  radar  system, 
and  the  other  was  a  Geophysical  Survey  Systems,  Inc.  (GSSI) 
short-pulse  radar.  Both  FM  systems  were  based  on  the  H?  85' 3 
microwave  network  analyzer.  The  radar  systems  used  for  the 
tests  are  described  below. 

Comparisons  between  the  FM  systems  and  the  GAR  short-pulse 
radar  were  made  using  antennas  built  by  GAR  specifically  for 
this  testing  program.  However,  due  to  the  use  of  an  integrated 
transmitter,  receiver,  and  antenna  assembly  in  the  GSSI  system, 
this  unit  could  not  be  used  with  the  GAR-built  antennas. 
Therefore,  it  was  not  possible  to  make  direct  comparisons 
between  the  GSSI  radar  and  the  others. 


6.1  THE  GAR  SHORT-PULSE  RADAR  SYSTEM 

The  GAR  short-pulse  radar  system  used  for  the  testing 
program  included  a  transmitter  pulser ,  equi val ent - t ime  sampling 
receiver,  antenna,  oscilloscope  display,  and  chart  recorder 
display.  A  block  diagram  of  the  GAR  short-pulse  radar  is  shcwr 
in  Figure  6.1  and  the  system  specifications  are  given  in  Table 
6.1. 

The  radar  uses  a  15  MHz  timing  reference,  from  which  all 
other  timing  signals  are  derived.  These  ti.ming  signals  are 
developed  on  the  radar  timing  board.  The  15  MHz  reference  is 
divided  down  to  a  5  MHz  pulse  repetition  frequency  ( ?RF )  and  a 


41 


l>n  1  .so  radar 


TABLE  6.1.  GAR  SHORT-PULSE  RADAR  SYSTEM  SPECIFICATIONS. 


Transmitte 

Pulse 

r« 

shape 

Sine  monocycle 

Pulse 

repetition  frequency  ( PRF ) 

5  MHz 

Pulse 

wi  dth 

1  nanosecond 

Pulse 

amplitude 

23  V  peak-to-peak 

Pulse 

peak  power 

8  watts 

Pulse 

average  power 

43  mw  (+16  dBm ) 

Sampler 

Type 

Single  diode  gate 
53  ohm  terminated 

Sample  strobe  width 

433  picoseconds 

Receiver 

Type 

Equivalent-t ime 

sampling 

Frequency  division  ratio 

533,030 

Scan  repetition  frequency 

50  Hz 

Real-time  receiver  window  width 

36  nanoseconds 

Minimum  discernible  signal 

2  mV,  terminated  i 

50  ohms 

Dynamic  range 

60  dB 

43 


50  Hz  scan  frequency.  The  timing  board  generates  a  50  Hz  scan 
ramp  that  is  used  in  the  receiver,  along  with  the  5  MHz  PR? 
signal,  to  generate  its  sweeping  sampler  trigger  signal. 

The  GAR  short-pulse  radar  uses  a  transmitter  that 
generates  a  one-nanosecond  pulse  shaped  as  a  single  cycle  of  a 
sine  wave.  The  transmitter  operates  at  a  PRF  of  5  .MHz  and 
produces  a  power  output  of  8  watts  peak  and  ^0  milliwatts 
average.  When  specified  in  decibels  referenced  to  a  milliwatt 
( d3m ;  ,  the  power  output  is  -r39  dBm  peak  and  +16  dBm  average. 

The  pulse  has  approximately  a  20  volt  amplitude  into  a  50  chm 
load.  The  transmitter  output  waveform  and  frequency  spectrum 
are  shown  in  Figure  6.2. 

The  receiver  is  an  equivalent-time  sampling  type  that 
generates  a  low-frequency  replica  of  the  received  RF  waveform. 
An  equi val ent - t ime  sampler  does  not  sample  a  signal  at  the  true 
Nyquist  frequency,  where  the  Nyquist  frequency  is  defined  as 
the  minimum  sampling  frequency  required  to  sample  a  waveform 
without  aliasing.  This  frequency  is  equal  to  twice  the  highest 
frequency  present  in  the  waveform.  Clearly,  since  the  GAR 
short-pulse  radar  generates  a  one-nanosecond  pulse,  the  true 
Nyquist  frequency  is  in  the  GHz  region. 

The  sampling  receiver  implements  a  500,000:1  frequency 
division  ratio.  The  receiver  output  waveform  occurs  at  a  50  Kz 
scan  rate,  or  a  20  millisecond  scan  period.  However,  the 
actual  sampled  waveform  duration  is  only  18  milliseconds,  due 
to  a  2  millisecond  reset  time  for  the  receiver.  The  real-time 
window  is  set  to  56  nanoseconds.  Thus,  a  one-nanosecond  pulse 
is  converted  to  a  500  microsecond  pulse  by  the  equivalent-time 
receiver , 

The  equi valent -t ime  sampling  technique  avoids  the  problems 
associated  with  an  extremely  hig.h  real-time  sampling  frequency 
by  utilizing  frequency  division.  The  only  requirement  is  that 
the  sampled  waveform  be  repetitive.  For  a  short-pulse  radar, 
this  is  equivalent  to  requiring  that  the  received  radar  signal 
be  stationary  in  time  over  a  receiver  scan  period. 

The  sampler  takes  one  sample  during  each  period  of  the  PRF, 
delaying  it  slightly  each  ti.me  relative  to  t.he  previous  one. 
Bach  successive  sample  modulates  the  charge  on  a  sampling 
storage  capacitor.  The  GAR  short-pulse  waveform  is  sampled 


44- 


Figure  6.2 


(A)  MONOCYCLE  TRANSMITTER  OUTPUT  WAVEFORM. 
VERTICAL  SCALE  =3.55  VOLTS/DIV. 
HORIZONTAL  SCALE  =  500  PS/DIV. 


(B)  FREQUENCY  SPECTRUM  OF  THE  ABOVE  PULSE. 
HORIZONTAL  SCALE  =  180  MHZ/DIV. 


GAR  short-pulse  radar  transmitter  (a)  output 
waveform  and  (b)  its  frequency  spectrum. 
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once  every  period  of  the  5  MHz  P?.F  for  IS  aii  1 1  i seconds .  Since, 
90,300  periods  of  the  PRF  exist  in  the  18  millisecond  sampled- 
waveform  window,  90,300  samples  are  required  to  replicate  a  36 
nanosecond  real-time  window.  The  equivalent  sample  spacing  of 
the  receiver  is  approximately: 

36  ns/90.000  =  0.4  picosecond,  (6.1) 

which  easily  meets  the  Nyquist  criterion. 

The  receiver  RF  sampler  is  a  50  ohm  terminated,  single- 
diode-gate  sampler  using  a  Schottky  sampling  diode.  The  sa.mple 
strobe  is  nominally  a  10  volt,  400  picosecond  pulse.  The 
.maximum  input  level  that  the  sampler  can  handle  is  4  volts 
peak-to-peak .  There  is  approximately  a  6  dB  conversion  loss 
through  t.he  sampler.  The  minimum  discernible  signal  at  the 
sampler  'output  is  2  millivolts.  Therefore,  the  sampler  dynamic 
range  is  63  d3 . 


6.2  THE  GSSI  SHORT-PULSE  RADAR  SYSTEM 

The  GSSI  short-pulse  radar  used  for  the  tests  with  the 
soil-model  test  tank  operates  in  a  similar  fashion  as  the  GAR 
short-pulse  radar.  However,  the  specifications,  timing 
functions,  and  operating  modes  for  the  GSSI  radar  system  are 
noc  as  well  known  by  GAR  personnel  as  those  for  the  GaR  short- 
pulse  radar.  During  tests,  the  radar  was  operated  wit.h  mi.ni.mu.m 
use  of  its  signal  processing  functions  in  order  to  make 
obiective  comparisons  with  the  other  GPR  systems. 

The  GSSI  radar  uses  a  one-nanosecond  transmitter  pulse  and 
an  equivalent- t ime  receiver.  The  transmitter  pulser  generates 
a  balanced,  one-nanoseccnd ,  unipolar  impulse.  It  operates  at  a 
system  ?RF  of  approximately  53  kHz.  Its  peak  power  output  is 
unxnown ,  but  is  comparable  to  the  GAR  monocycle  transmitter 
power  output.  Its  average  power  output,  however,  is  less  than 
the  GAR  transmitter  because  it  operates  at  a  PRF  that  is  133 
times  lower . 

The  GSSI  receiver  is  a  balanced  diode-bridge  sampler  that 
operates  in  a.n  equi  val  ent  - 1  ime  mode.  It  has  an  adjustable 
real-time  wi.ndow  and  a.n  adjustable  output  waveform  scan  rate. 
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The  radar  was  adjusted  for  the  tests  to  operate  with  a  scan 
frequency  of  51.2  Hz  (a  19.5  millisecond  period).  The 
equivalent-time  frequency  division  ratio  was  adjusted  to 
520,000:1  to  match  the  GAR  system.  The  GSSI  receiver  has  a 
reset  time  of  about  4.00  microseconds.  Therefore,  the  real-c 
window  is  39  nanoseconds  long.  Based  on  a  50  kHz  PR?,  977 
samples  are  used  to  generate  the  replica  waveform.  This 
translates  into  a  40  picosecond  equivalent-time  sample  spaci 
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Figure  6.3  is  a  block  diagram  of  the  GSSI  short-pulse  G 
radar  timing  and  control  signals  are  ge.nerated  in  the 
rol  unit.  Trigger  signals  are  sent  to  the  transmitter  a 
ling  receiver  through  a  cable.  The  sampler  output  wavef 
ent  by  cable  to  the  control  unit  for  further  processing, 
receiver  output  waveform  can  be  displayed  on  an 
lloscope  and  a  strip  chart  recorder. 


The  transmitter  and  sampler  are  housed  with  their 
respective  antennas  in  an  integrated  transducer  enclosure, 
radar  operates  in  a  bistatic  mode  using  tr iangular-sheet  bow 
tie  antennas.  The  exact  construction  of  the  receiver, 
transmitter,  and  antennas  is  not  known  because  it  is  not 
possible  to  gain  access  to  any  of  the  components  in  the 
transducer . 


6.3  THE  FREQUENCY-DOMAIN  GPR  SYSTEMS  USING  THE  HP  8510  NETWORK 
ANALYZER 


This  section  gives  an  overview  of  the  Hewlett-Packard  S' 
network  analyzer  system  and  describes  its  application  to  che 
GPR  signal  source  testing  program.  An  explanation  of  all  cf 
tne  functions  and  capabilities  of  the  HP  8510  is  beyond  the 
scone  of  this  docume_nt  .  However,  Appendix  C  contains  an 
excerpt  from  the  HP  8510  manual  and  provides  a  good 
introduction  amd  syste.m  overview.  Appendix  C  also  contains 
simplified  block  diagram  of  the  network  analyzer  system.  Th 
present  section  only  addresses  those  functions  of  the  HP  85' 
that  are  directly  applicable  to  the  GPR  measureme.nt  program.. 

The  HP  8510  was  used  in  two  modes  for  the  GPR 
measurements.  The  first  was  a  swep t - f r equency  mode  that  use 
synthesized  sweeper  i.n  its  linear  frequency  sweeping  m.cde  t.i 
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generate  the  increasing-f requency  sweep  of  RF  signal 
frequencies  required  for  the  measurements.  This  mode 
represented  a  simulation  of  a  swept  FJ-l-CW  GPR  system.  The 
second  mode  of  operation  for  the  HP  8510  simulated  a  sxepped-FF 
GPR  by  using  a  stepped-f requency  signal  source. 

Section  6.3.1  below  explains  some  of  the  common 
capabilities  of  the  HP  8510  that  were  used  in  both  the  FM-CV 
and  stepped-F?-l  system  tests.  Section  6. 3.  1.3  describes  the 
differences  between  the  swept  and  stepped  operating  modes  of 
xhe  HP  8510. 


6.3.1  HP  8510  NETWORK  ANALYZER  CAPABILITIES 

The  Hewlett-Packard  8510  network  analyzer  is  a  precision 
instrument  used  for  making  RF  and  microwave  measurements  of 
two-port  networks.  The  8510  measurement  system  includes  the 
following  four  subsystems:  a  frequency  source,  a  test  set,  a 
signal  detector  and  analog-to-digital  converter,  and  a 
microprocessor  ajnd  display.  The  HP  8510  system  conf igur ax  ion 
used  by  GAR  for  this  program  included  an  HP  834-lA  synthesized 
sweeper,  an  HP  851 3A  transmission  test  set,  an  HP  851 02A  IF 
detector,  and  an  HP  851 01 A  display  processor.  This  HP  8510 
configuration  allows  measurements  to  be  made  over  the  frequency 
range  from  4-5  MHz  to  20  GHz. 


6, 3. 1,1  SCATTERING  PARAMETER  MEASUREMENTS  WITH  THE  HP  8510 

The  HP  8510  is  designed  to  make  vector  (signal  amplitude 
aind  phase)  scattering  parameter  ( S-par  ameter )  measurements  of  a 
two-port  electrical  network.  Measurements  of  devices  at  RF  an„ 
microwave  frequencies  primarily  use  S-parameters  because  S- 
parameters  are  well-suited  to  describing  the  performance  of 
xransmissicn  line-based  circuits  and  systems.  The  S-par  a.mexer  s 
are  quite  useful  for  characterizing  networks  at  RF  and 
microwave  frequencies  because  they  are  directly  related  to  xhe 
amount  of  power  reflected  at  the  test  network  ports  and  the 
amount  of  power  transmitted  through  the  device.  The  derivation 
and  advanced  applications  of  S-parameters  are  described  in 
texts  treating  the  subject  of  RF  and  microwave  applications-. 
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The  S-parameters  are  often  normalized  to  impedances  of  52 
ohms.  This  implies  that  the  netvork  has  an  input  impedance  and 
output  impedance  of  50  ohms.  Thus,  any  connecting  transmission 
lines  and  test  cables  should  have  50  ohm  characteristic 
impedances  as  well.  GPR  test  equipment  configurations 
including  the  H?  8510  used  50  ohm  coaxial  cables  to  connect  the 
antennas  to  the  test  set  ports.  Also,  the  antennas  have 
nominal  input  impedances  of  50  ohms. 

The  two-port  network  measurements  using  S-parameters  are 
sub-divided  into  those  related  to  transmission  through  the 
network  (forward  and  reverse),  and  those  related  to  reflections 
at  the  two  network  ports.  Transmission  measurements  include 
insertion  loss  or  gain,  transmission  coefficient,  electrical 
delay,  and  group  delay.  Reflection  measurements  include  return 
loss,  voltage  standing  wave  ratio  (VSWR),  reflection 
coefficient,  and  impedance. 

The  measurement  quantities  of  primary  interest  in  this 
testing  program  are  the  reflection  coefficient,  S, ,  ,  and  the 
transmission  coefficient,  S21  .  If  50  ohm  impedance 
normalization  is  used,  S,,  is  the  ratio  of  the  reflected 
voltage  to  the  incident  voltage  at  port  one  of  the  two-port 
network  under  test.  The  measurement  of  S,,  requires  that  port 
two  of  the  network  be  properly  terminated  while  port  one  is 
being  driven  by  a  signal  source  with  a  50  ohm  driving 
impedance.  Sj-  is  the  ratio  of  the  voltage  measured  at  port 
two  of  the  network  to  the  incident  voltage  at  port  one,  with 
port  one  being  driven  by  the  signal  source.  The  50  ohm 
termi.nation  and  driving  impedance  requirements  apply  to  the  S,. 
measurement  as  well. 

If  an  antenna  is  connected  to  port  one  of  the  HP  851 3A 
test  set,  the  parameter  S,,  represents  a  measurement  of  the 
energy  reflected  by  the  antenna  itself  and  any  targets  that 
reflect  radiated  energy  back  to  the  antenna.  This 
conf igurat ion  represents  both  a  one-port  network  and  a 
monostatic  radar  system. 

If  an  antenna  is  also  connected  to  port  two  of  the  test 
set,  the  parameter  S21  represents  the  energy  received  at  port 
two  through  the  two-antenna  network.  The  energy  received  at 
port  two  includes  reflections  from  targets  irradiated  by  the 
port  one  ante.nna.  This  conf iguration  represe.nts  a  two-port 
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network  and  a  bistatic  radar  system.  The  measurements 
conducted  with  the  soil-model  tank  were  performed  with  a 
bistatic  arrangement,  measuring  the  Sj,  parsimeter  . 

When  S,i  and  Sj,  measurements  are  made  with  the  network 
analyzer,  the  RF  frequency  range  of  interest  must  be  specified. 
Whether  in  a  ramp  mode  or  a  stepped  mode,  the  network  analyzer 
adjusts  the  frequency  of  its  signal  souroe  and  makes  S- 
parameter  measurements  at  a  number  of  individual  frequencies 
throughout  the  specified  bandwidth.  The  number  of  freque.ncies 
is  sej-ectable,  with  the  available  choices  being  5"  ,  '3'  ,  2  3":  . 

and  iiSl  .  The  effects  of  choosing  a  particular  number  of 
measurement  frequencies  are  discussed  later. 

Once  the  frequency  range,  the  number  of  frequency  points, 
and  the  ramp  or  stepped  mode  are  selected,  the  network  analyze, 
requires  calibration.  Calibrating  the  network  analyzer  reduce: 
measurement  errors.  The  vector  calibration  identifies  the 
measurement  reference  planes  at  the  two  network  ports  and 
corrects  for  any  losses  or  phase  aberrations  up  to  the 
reference  planes.  Thus,  the  calibration  procedure  can  be  used 
to  establish  the  reference  planes  for  a  Sji  measurement  at  the 
ends  of  the  antenna  cables.  The  HP  8510  then  corrects  for  any 
cable  loss  and  phase  effects. 

The  network  analyzer  calibration  procedure  requires 
measurements  of  known  terminations  at  the  reference  planes 
before  the  test  network  is  inserted  for  measurements.  The 
standard  terminations  include  a  shielded  open  circuit,  a  shore 
circuit,  and  a  50  ohm  load  (a  50  ohm  sliding  load  termi.naeicn 
is  required  if  precise,  high-frequency  measurements  are  to  be 
made).  Once  the  standard  terminations  are  measured  and  the 
network  analyzer  computes  its  vector  correction  factors,  all 
subsequent  measurements  automatically  apply  the  correction 
factors  to  the  data.^  The  HP  8510  uses  an  elaborate  S-para.meo9' 
error  model  to  compute  the  correction  factors,  the  details  cf 
which  are  beyond  the  scope  of'  this  document.  However,  the  H? 
3510  manual  contains  a  complete  d-^scr ipt ion  of  the  vector 
error-correction  model. 
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6. 3. 1.2  HP  8510  TIME-DOMAIN  MEASUREMENT  MODE 


The  KP  8510  can  display  the  results  of  its  measurements  in 
the  frequency-domain  or  the  time-domain.  However,  all  of  its 
measurements  are  made  in  the  frequency-domain.  The  HP  8510  has 
a  time-domain  bandpass  mode  of  operation  that  synthesizes  an 
impulse  function  with  a  pulse  width  inversely  proportional  to 
the  RF  bandwidth  specified  for  the  measurement.  A  type  of 
digital  Fourier  transform  known  as  the  Chirp-Z  transform  is 
used  to  transform  measur emenzs  of  S,,  and  Sj.  from  a  frequency- 
domain  representation  to  a  time-domain  representation  . 

Appendix  D  contains  excerpts  from  the  HP  8510  manual  explaining 
the  time-domain  synthesis  mode  of  operation. 


6. 3. 1.2.1  TIME-DOMAIN  PULSE  SYNTHESIS  WITH  WEIGHTING  FUNCTIONS 

Once  an  RF  bandwidth  is  selected  for  measurements,  and  the 
time-domain  synthesis  mode  is  selected,  frequency-domain 
weighting  functions  are  applied  to  adjust  the  sidelobe  levels 
of  the  synthesized  impulse.  Three  selectable  weighting 
functions  are  provided  which  give  maocimum  sidelobe  levels  cf  - 
13  dB ,  -43  d3 ,  and  -90  d3  below  the  pulse  main  lobe.  These 
weighting  functions  are  identified  as  "minimum,"  "normal,"  and 
"maximum,"  respectively.  For  a  given  bandwidth,  the  cost  of 
increased  sidelobe  suppression  is  an  increased  pulse  width  for 
the  synthesized  pulse. 

Figures  6.4  through  6.6  are  magnitude  plots  of  one- 
nanosecond  synthesized  pulses  generated  by  the  network  analyzer 
using  the  three  frequency-domain  weighting  functions,  or 
windows.  The  pulse  i.n  Figure  6.4  uses  the  minimum  window  with 
a  45  MHz  to  1.2  GHz  bandwidth;  the  pulse  in  Figure  6.5  uses  the 
normal  window  with  a  45.  MHz  to  2.0  GHz  bandwidth:  and  the  pulse 
in  Figure  6.6  uses  the'  maximum  window  wit.h  a  45  MHz  to  2.8  GHz 
bandwidth.  All  of  the  plots  use  a  one-nanosecond  per  division 
horizontal  scale  and  t.he  pulse  magnitude  is  normalized  to  o.ne 
unit,  with  a  linear  vertical  scale.  Note  that  more  RF 
bandwidth  is  required  to  maintain  a  given  pulse  width  with 
suppressed  sidelobes. 

All  network  analyzer  measurements  performed  with  the  soil- 
model  test  tank  used  a  one-nanosecond  pulse  having  a  *i5  MHz  to 
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Figure  6. 
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HP85iO  plot  of  the  magnitude  of  a  Ins  pulse 
synthesized  using  the  minimum  weighting  function 
and  a  45MHz-1.2GHz  bandwidth. 
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Figure  6.5;  ::?3C10  plot  of  ^-he  magnitude  of  a  Ins  pulse 

synthesized  using  the  normal  weighting  function 
and  a  45MHz-2.0GHz  bandwidth. 
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2.(3  GHz  bandwidth,  and  a  normal  window.  The  one-nanoseoond 
pulse  width  was  selected  to  permit  a  more  accurate  comparison 
of  the  network  analyzer  system  performance  with  the  performa.nce 
obtained  with  xhe  pulsed  GPR  systems.  In  addition,  GAR  felt 
that  the  bandwidth  required  for  the  maximum  window  was  too 
great,  and  the  bow-tie  antennas  used  for  the  tests  oould  not 
support  the  extended  bandwidth.  Also,  the  soil  attenuation 
properties  would  oause  greater  attenuation  of  the  higher 
frequencies,  thus,  widening  the  pulse. 


6. 3. 1.2. 2  EFFECT  OF  THE  NUMBEK  OF  FREQUENCY  MEASUREMENT  POINTS 
SELECTED 

As  indicated  earlier,  the  HP  8510  network  analyzer 
performs  its  measurements  at  a  finite  number  of  points  over  the 
bandwidth  specified  by  the  user.  The  number  of  points  is 
selectable  from  51,  101,  201,  and  401.  When  the  network 

analyzer  is  operated  in  the  time-domain  bandpass  mode,  it 
generates  the  same  number  of  measurement  points  as  specified  by 
the  user  during  the  HP  8510  calibration  in  the  frequency- 
domain  . 

The  number  of  measurement  points  and  the  bandwidth 
selected  determine  the  unambiguous  range  in  the  time-domain 
mode.  The  unambiguous  range  is  the  maximum  range  time  for 
which  the  time-domain  waveform  will  not  repeat.  This  is 
analogous  to  the  reciprooal  of  the  PRF  for  a  time-domain  pulse 
radar  system.  The  unambiguous  range  for  the  network  analyzer 
is  equal  to  the  reoiprocal  of  the  frequenoy  spacing.  For 
example,  if  201  points  are  seleoted,  and  a  frequenoy  range  of 
45  MHz  to  2.0  GHz  is  specified,  the  unambiguous  range  is  (201- 
1)/(2,0  GHz-45  MHz)  =  102.3  nanoseconds,  which  corresponds  to 
slightly  more  than  25  feet  in  soil  with  a  dielectric  constant 
of  1  6  . 

The  selection  of  the  number  of  measurement  points  also 
determines  the  range  measurement  accuracy  of  the  network 
analyzer.  For  the  GPR  measurements,  the  range  window  was 
selected  to  be  36  nanoseconds  long  and  the  number  of 
measurement  points  was  selected  to  be  201.  Using  these 
quantities  as  an  example,  the  range  measurement  accuracy  is 
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given  by  36  ns/(20l-1)  =  180  picoseconds,  which  corresponds  to 
about  0.5  inch  in  soil  with  a  dielectric  constant  of  16. 

Another  important  factor  affected  by  the  selection  of  the 
number  of  measurement  points  is  the  time  required  to  complete 
measurement  over  tne  specified  bandwidth.  Increasing  the 
number  of  measurement  points  requires  an  increase  in  the 
measurement  time.  Therefore,  doubling  the  number  of  points 
from  201  to  4-01  will  require  approximately  twice  the  time  to 
complete  the  measurement. 


6.3. 1.3  DIFFERENCES  BETWEEN  THE  HP  8510  SWEPT -FREQUENCY  AND 
STEPPED-FREQUENCY  MEASUREMENT  MODES 

The  H?  8510  network  analyzer  can  perform  its  measurements 
over  a  specified  frequency  band  by  stepping  its  frequency 
source  in  discrete  frequency  steps,  or  by  sweeping  its  source 
continuously.  Stepping  the  source  provides  the  highest 
measurement  accuracy  and  dynamic  range.  In  the  stepped  mode, 
the  HP  834.1  A  synthesizer  is  stepped  and  phase-locked  to  each 
measurement  frequency.  This  reduces  -ne  phase  noise  effects 
and  the  frequency  inaccuracies  that  appear  in  the  swept  mode. 

When  the  swept  mode  is  used,  the  synthesizer  is  locked  to 
the  start  frequency  and  swept  linearly  to  the  stop  frequency. 
However,  the  synthesizer  is  not  locked  to  any  other  frequency 
besides  the  start  frequency.  The  network  analyzer  in  effect 
assumes  that  the  synthesizer  will  be  at  certain  frequency 
points  at  certain  times  during  its  sweep.  It  has  no  way  to 
determine  the  exact  frequency  of  the  synthesizer.  Thus, 
inaccuracies  in  determining  the  exact  frequency  of  the 
synthesizer  will  exist  as  it  is  swept  across  the  specified 
bandwidth.  The  primary  advantage  to  the  swept-f r equency 
measurement  moae  is  its  speed.  In  the  swept  mode,  the  network 
analyzer  produces  a  measurement  on  the  order  of  i  00  times 
faster  than  in  the  stepped  mode. 
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6.3.2  SUMMARY  OF  NETWORK  ANALYZER  CONFIGURATIONS  AND 

SPECIFICATIONS  FOR  SOIL-MODEL  TANK  MEASUREMENTS 

All  measurements  performed  with  the  network  analyzer  on 
the  soil-model  test  tank  used  the  HP  8510  time-domain  bandpass 
mode.  The  time-domain  range  window  for  the  shale  olay 
measurements  was  set  to  36  nanoseconds  in  duration.  It  was 
changed  to  4-0  nanoseconds  for  the  clay-schist  measurements 
because  the  sohist  had  a  higher  relative  dielectric  constant 
than  the  shale.  The  frequency  range  of  45  MHz  to  2.0  GHz  with 
201  points  was  selected.  A  normal  weighting  function  was  used 
to  synthesize  a  one-nanosecond  pulse  with  maximum  sidelobe 
levels  of  -43  dB . 

The  sweep  time  was  set  for  100  milliseconds.  In  the  swept 
measurement  mode,  the  45  MHz  to  2  GHz  sweep  takes  place  in  100 
milliseconds.  However,  in  the  stepped  measurement  mode,  a 
measurement  takes  approximately  100  times  longer,  or  10 
seconds . 

The  measurements  utilized  a  bistatic  antenna 
conf iguration .  Hence,  the  S21  scattering  parameter  was  of 
primary  interest.  Recall  that  S21  measures  the  insertion  loss 
of  a  two  port  network,  which  is  essentially  what  a  bistatic  GPR 
does  when  it  radiates  energy  from  one  antenna  and  receives 
energy  at  another. 

The  HP  8510  RF  output  power  was  set  to  +10  dBm  (0.1  watt) 
average  for  all  frequencies  over  the  measurement  range.  At 
this  power  level,  tne  network  analyzer  had  approximately  a  90 
dB  dynamic  range.  However,  the  full  dynamic  range  capabilities 
of  the  network  analyzer  were  not  utilized.  This  was  because 
the  surface  return  signal  from  the  S21  measurements  with  the 
GPR  antennas  was  about  40  dB  below  the  calibration  signal 
level.  Thus,  the  usable  dynamic  range  was  reduced  by  40  dB. 

The  dynamic  range  mig.ht  be  i.mproved  at  the  expense  of 
calibration  accuracy  by  i.nserting  a  low-noise  amplifier  at  port 
two  (the  receiver  port)  of  the  HP  8513A  transmission  test  set. 
Unfortunately,  calibrating  the  analyzer  wit.h  the  amplifier 
would  require  adding  an  attenuator  in  front  of  the  amplifier 
during  calibration  to  avoid  overloading  the  port  two  i.nput  . 

The  attenuator  would  then  be  removed  during  actual 
measurements.  However,  the  vector  calibration  would  then  be 
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inaccurate.  It  was  decided  not  to  use  an  amplifier  at  the 
receiver  port  to  maintain  the  accuracy  of  the  calibration. 

The  network  analyzer  has  the  capability  of  performing 
trace  averaging.  Averaging  traces  reduces  the  noise  effects, 
and  hence  increases  the  measurement  dynamic  range.  The 
measurement  time  in  the  swept  mode  increases  linearly  with  th 
number  of  averages,  but  in  the  stepped  mode,  at  least  530 
averages  are  required  to  lengthen  the  measur e.ment  time.  It  w 
decided  not  to  utilize  the  averaging  capabilities  of  the  H? 

35 '0  because  the  pulsed  G?R  systems  tested  did  not  have 
averaging  as  a  built-in  feature.  Utilizing  averaging  for  onl 
the  frequency-domain  G?R  systems  would  not  allow  an  ob;;;ective 
comparison  between  all  of  the  tested  radars. 


6.4  ANTENNAS  USED  FOR  THE  GPR  SOIL-.MODEL  TESTS 

Two  types  of  antennas  were  designed  and  built  for  testing 
the  GPR  systems.  One  type  was  a  TSM  horn  si.mila,’*  to  the  ones 
used  with  the  GAR  RCDAR*^  highway-survey  GPR.  The  other  was  a 
tr iangul ar -sheet  (bow-tie)  dipole.  Each  of  these  antennas 
offered  advantages  over  the  other.  Tests  were  performed  wich 
the  two  types  of  antennas  to  decide  which  type  gave  the  best 
overall  performance. 


6.4.1  TE.M  HORN  ANTENNA 

The  TEM  horn  has  several  performance  capabilities  that 
make  it  useful  for  GPR  applications.  It  has  broad  bandwidth, 
low  dispersion,  and  a  low  reflection  coefficient.  Such  facccr 
make  it  ideal  for  short-pulse  radars,  which  require  low  pulse 
stretch  and  hig.h  ba.ndwidth.  Some  of  t.he  considerations 
impacting  the  design  of  a  TEM  antenna  for  GPR  applications  are 
outli.ned  in  Section  5.3.3. 

GAR  spent  considerable  ti.me  constructing  a  pair  of  TE.M 
horns  for  use  in  t.he  tests.  We  used  a  design  similar  to  the 
.RCDAR*^  a.ntennas,  except  t.hat  a  wider  aperture  and  a  different 
shape  for  the  antenna  foil  plates  were  employed.  The 
dimensions  of  the  TEM  antennas  were  selected  for  optimum 
operation  with  the  cne-nanosecond  time-domain  pulses  that  wer-- 
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used  with  the  test  GPR  systems.  Furthermore,  we  decided  to 
an  alternative  fabrication  technique  for  the  test  antennas  i 
order  to  see  if  enhanced  antenna  clutter  and  reflection 
coex'f i c i e.n t  performance  could  be  achieved. 

Each  antenna  was  fabricated  with  two  39-inch  long  brass 
foil  plates  that  had  varying  width  cross-sections  along  the! 
lengths.  The  foil  plates  were  adjusted  into  a  wedge  shape  w 
the  feed  point  at  the  narrow  end.  The  foil  plates  were  glue 
to  two  0.75-inch  thick  foa.m  boards  and  spaced  with  the  prcpe 
wedge  angle  by  additional  foam  side  panels.  The  plates  face 
each  Qt.ner  on  the  inside  of  the  wedge.  The  i.ntention  was  to 
have  only  air  between  the  foil  plates  to  minimize  adverse 
performance  anomalies  that  could  result  with  ot.her  materials 
betwee.n  the  plates. 

At  the  feed  toint,  a  coaxial  connector  made  the  transit! 
to  the  parallel  plate  structure  of  t.he  antenna.  The  tips  of 
the  antenna  plates  were  attached  to  resistive  cards  that 
terminated  the  surface  currents  flowing  on  the  plates.  The 
a.nte.nna  feed  point  and  tip  were  tuned  to  reduce  time-doraai.n 
reflections.  This  tuning  procedure  was  a  time-consuming  tas 
requiring  considerable  patience. 

The  antennas  were  tested  using  both  short-pulse  time- 
domain  ref 1 ectometry  and  HP  8510  frequency-domain  tests  to 
verify  their  performance.  We  found  that  the  antennas  had 
reflection  waveform  amplitudes  that  were  about  -36  d3  down  f 
tne  incident  pulse  amplitude.  Their  effective  radiation 
bandwidth  extended  from  250  MHz  to  2.0  GHz.  However,  time 
limitations  did  not  allow  formal  and  extensive  antenna 
parameter  and  pattern  tests  to  be  performed. 


6.4.2  TRIA.VGULAS- SHEET  DIPOLE  ANTENNA 

The  ccplanar,  tr  ia.ngular -sheet  dipole  a.ntenna  is  a  simpl 
iesign  that  has  been  used  extensively  in  GPR  applications, 
e.xhibits  a  .moderately  broadband  response,  and  achieves  it  wi 
a  simple,  compact  structure.  It  does  not,  however,  e.xhibit 
wide  bandwipth,  low  dispersion,  or  low  reflection  coefficisn 
of  the  TEM  horn. 


These  antennas  are  deployed  with  the  two  triangular  sheet.: 
lying  on  the  surface  of  the  ground.  If  the  antenna  is 
electrostatically  shielded  on  the  back  side,  most  of  the 
radiated  signal  energy  from  the  antenna  couples  into  the 
ground.  The  close  coupling  of  the  dipole  to  the  ground  gives 
it  a  tremendous  advantage  over  the  TE?^  horn. 

The  length  of  the  dipole  arms  is  directly  related  to  the 
radiation  bandwidth  of  the  antenna.  Long  triangular  sheets  ca: 
radiate  lower  frequency  RF  energy  than  short  triangular  sheets 
Typically,  each  ante.nna  arm  .must  be  at  least  1/4.  wavele.ngth 
long  at  the  lowest  desired  RF  frequency  in  order  to  radiate 
properly.  Fortunately,  when  the  bow-tie  antenna  is  in  direct 
contact  with  the  ground,  the  dielectric  constant  of  the  soil 
acts  to  decrease  the  wavelength  of  the  RF  energy,  and  thus 
appears  to  make  the  antenna  arms  effectively  longer.  This 
implies  that  for  a  given  RF  bandwidth,  the  GPR  bow-tie  antenna 
dimensions  can  be  reduced  from  ‘the  di.mensions  of  an  antenna 
designed  to  radiate  in  air. 

GAR  designed  one  pair  of  bow-tie  antennas  for  the  GPR 
tests.  The  antennas  were  designed  to  have  a  1  ow-f r eque.ncy 
cutoff  of  approximately  230  MHz  when  operated  over  soil  with  a 
relative  dielectric  constant  of  16.  The  design  parameters  of 
the  antenna  were  extracted  from  a  report  published  by  Brown  an 
Woodward'.  The  bow-tie  was  built  from  brass  shim  stock  and 
i.ncluded  a  shield  on  the  back  to  minimize  clutter  returns  and 
maximize  the  e.nergy  coupled  i.nto  the  ground. 

A  sketch  of  the  bow-tie  antenna  is  shown  in  Figure  6.". 
coaxial  cable  feed  is  coupled  to  the  two  antenna  leaves  throug 
a  50-ohm  to  233-chm  balun  transformer.  Each  of  the  antenna 
leaves  has  two  230-ohra  resistors  connecti.ng  the  antenna  foils 
to  the  antenna  shield.  These  resistors  help  to  absorb  the 
e.nergy  tnat  is  not  "radiated  or  reflected  by  the  antenna. 


6.4.3  QUALITATIVE  COMPARISONS  BETWEEN  THE  TEM  HORN  ANTENNAS 

AND  THE  TRIANGULAR -SHEET  BOW-TIE  ANTENNAS 

Qualitative  tests  were  used  to  compare  the  performances 
t.ne  TEM  horns  a.nd  the  bow-tie  anten.nas  in  order  to  determine 
wnic.n  type  of  a.nten.na  to  use  for  the  GPR  source  tests.  The 
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antennas  were  tested  in  the  ground-penetration  mode  on  the 
soil-model  tank  using  the  responses  from  the  buried  targets  an 
loop  probes  as  the  criteria  of  selection.  No  formal  antenna 
pattern,  VSWR,  or  pulse  tests  were  performed.  The  antennas 
were  tested  in  the  bistatic  and  monostatic  modes  with  both  the 
GAR  time-domain  G?R  and  the  HP  85.0  network  analyzer -based 
systems.  The  TEM  horn  antennas  were  tested  in  a  gr ound-ccntac 
mode  and  an  air-coupled  mode.  Most  of  the  expected  performanc 
characteristics  of  the  antennas  were  verified  during  the  tests 

The  TEM  horns  exhibited  good  bandwidth  response  and  a  lew 
reflection  coeff icie.nt  in  both  the  ground-contact i.ng  and  air- 
coupled  modes.  The  target  returns  and  probe  signals  showed 
little  pulse  stretch.  The  response  of  the  bow-tie  antennas 
showed  more  pulse  stretch  sind  less  bandwidth.  The  bow-tie 
antennas  had  a  higher  reflection  coefficient  than  did  the  TEM 
horn  antennas.  The  bow-tie  antennas  also  coupled  more  energy 
into  the  ground  than  did  the  TE?^  horns.  This  was  evidenced  by 
the  higher  target  return  signals. 

An  unexpected  result  led  to  the  selection  of  the  bow-tie 
antennas  for  the  subsequent  GPR  tests.  Specifically,  the  bow- 
tie  antennas  had  a  signal-to-clutter  ratio  that  was  at  least  1 
d3  larger  than  that  of  the  horns.  This  phenomenon  was  probabl 
due  to  the  better  antenna- to-ground  energy  coupling  of  the  bow 
tie  antennas.  In  addition,  the  TEM  horns  radiate  along  their 
entire  39-inch  structure  and  thus  are  susceptible  to  above¬ 
ground  clutter. 
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SECTION  7 


DESCRIPTION  OF  THE  MEASUREMENTS  USLNG  THE  CANDIDATE  GPR  SYSTEMS 

The  GPR  measurement  program  was  divided  into  three  phases. 
Phases  1  and  2  involved  scale-model  measurements  with  the  soil 
tank  using  two  different  claj  dielectrics.  Phase  3  consisted 
of  full-scale  GPR  measurements  using  the  HP  851(3  in  a  stepped- 
frequencj  mode  on  an  outdoor  -cest  field.  The  stepped-f requencj 
system  was  chosen  exclusively  for  the  Phase  3  measurements 
because  it  showed  the  best  performance  of  the  four  scale-model 
GPR  systems  tested. 


7.1  DESCRIPTION  OF  THE  PHASE  1  AND  PHASE  2  SCALE-MODEL  GPR 
MEASUREMENTS  USING  THE  SOIL-MODEL  TANK 

Phase  1  involved  soil-model  tank  measurements  using  clay 
shale  as  the  model  dielectric.  The  clay  shale  provided  the 
most  demanding  of  the  two  scale-model  testing  environments. 

This  was  because  the  shale  had  considerably  higher  loss 
characteristics  than  the  clay  schist  used  in  the  Phase  2 
measurements.  Phase  1  and  Phase  2  measurements  were  conducted 
in  a  similar  fashion.  Figures  7.1  through  7.3  are  photographs 
of  some  of  the  GPR  equipment  configurations  used  for  the  soil 
tank  measursmer: 0  3  . 

When  the  target  measurements  were  made,  the  bistatic 
antenna  systems  were  placed  directly  over  eac.h  of  the  pipes. 

The  antennas  were  oriented  in  such  a  manner  that  the  electric 
field  polarization  vector  was  parallel  with  the  pipe  axial 
dimension  to  ensure  maximum  signal  return.  Recall  that  the 
GSSI  short-pulse  radar  required  the  use  of  its  own  antenna 
system,  so  the  spec ial ly-constructed  test  antennas  could  not  be 
used  with  it.  This  should  be  taken  into  account  when  comparing 
the  results  of  the  GSSI  radar  measurements  with  those  from  the 
other  GPR  systems. 

The  output  data  of  each  of  the  candidate  GPR  system.s  were 
presented  in  a  time-domain  format  so  that  the  results  could 
easily  be  compared.  Photographs  were  taken  of  the  GPR  time- 
domain  displays.  In  the  case  of  the  HP  a5i3-based  frequency- 
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domain  systems,  plots  were  also  made  of  the  time-domain 
displays  for  better  visual  clarity. 

Care  was  taken  to  avoid  any  uncalibrated  gain  adjustments 
with  the  GPR  systems  in  order  to  maintain  measurement  accuracy 
and  repeatability.  The  target  signals  from  both  of  the  short- 
pulse  radars  were  measured  at  the  receiver  outputs.  No  gain 
adjustments  were  possible  in  the  receiver  signal  chain  before 
the  measurement  points  with  the  short-pulse  radars.  The 
network  analyzer  calibration  procedure  performed  before  each 
measurement  removed  the  possibility  of  any  uncalibrated  gain 
adjustment  errors  in  its  display. 

The  performances  of  the  GPR  systems  were  not  compared 
based  on  absolute  received  signal  power  levels.  Instead,  the 
primary  performance  criteria  used  for  the  GPR  comparisons  were 
based  on  target  detectability  and  the  signal-to-clutter  levels 
of  the  target  returns. 

Strip-chart  recordings  were  produced  for  the  two  short- 
pulse  radars  as  the  antenna  systems  were  pulled  across  the  soil 
tank.  The  strip-chart  recordings  were  intended  primarily  as  a 
qualitative  indication  of  the  relative  performances  of  the  two 
short-pulse  radars.  They  were  most  useful  when  searching  for 
target  return  signals  having  signal-to-clutter  ratios  less  than 
6  dB. 

A  strip-chart  recording  provides  additional  information  to 
the  GPR  operator  by  displaying  the  radar  scans  on  paper  in  a 
sequential  fashion.  The  radar  output  waveform  is  digitized  and 
converted  to  a  grey-level  format  and  printed  on  paper  by  a 
stylus.  The  paper  scrolls  forward,  and  as  each  radar  scan  is 
produced,  the  stylus  scans  across  the  scrolling  paper  in  a 
direction  perpendicular  to  the  scrolling  direction,  thereby 
producing  a  line  with  modulated  grey-level  intensity.  As  more 
and  more  scans  are  recorded  on  the  paper,  the  strip-chart 
recording  presents  a  display  that  allows  the  observer  to 
visually  "integrate"  many  scans  to  see  the  changing  target 
returns . 

Unfortunately,  the  HP  8510  receiver  waveform  data  was  not 
readily  available  for  output  to  a  strip-chart  recorder.  To 
generate  a  strip-chart  recording  of  the  HP  8510  data,  the 
digitized  display  output  would  have  to  be  exported  to  a 
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computer,  stored,  and  then  output  to  the  strip-chart  recorder. 
Time  did  not  permit  implementation  of  a  strip-chart  interface 
for  the  HP  8510. 

Simple  signal-to-clutter  measurements  were  performed  with 
each  of  these  candidate  GPR  systems.  The  purpose  of  these 
measurements  was  to  evaluate  the  internal  clutter  level 
generated  in  the  GPR  systems.  However,  this  clutter  level  is  a 
function  of  the  antennas  and  their  deployment  because  any 
impedance  mismatches  at  the  antenna  cause  reflections  that 
appear  in  the  receiver  range  window.  It  was  desirable  to 
include  this  form  of  clutter  in  the  signal-to-clutter 
measurement . 

It  was  very  difficult  to  structure  an  accurate  test  for 
the  signal-to-clutter  measurement.  The  bow-tie  antennas  are 
very  sensitive  to  radiation  impedance  loading.  They  must  be 
placed  on  or  near  the  surface  of  the  ground  to  avoid  producing 
a  response  with  excessive  ringing.  The  excess  ringing 
contributes  to  the  clutter  present  in  the  receiver  output 
waveform.  Ideally,  the  antenna  must  be  placed  over  a 
dielectric  medium  that  absorbs  all  of  the  energy  radiated  from 
the  antenna;  thus,  no  signals  from  external  targets  and  signals 
will  appear  in  the  receiver  waveform. 

GAR  engineers  were  unable  to  locate  such  an  ideal  medium. 
The  closest  available  approximation  to  this  ideal  absorbing 
medium  that  could  be  used  for  the  signal-to-clutter  tests  was 
the  clay  shale  in  the  Phase  1  test  tank  (the  clay  schist  in  the 
Phase  2  tank  did  not  have  high  enough  attenuation).  The  shale 
exhibited  a  very  high  attenuation  characteristic.  Although  the 
shale  tank  contained  pipe  targets,  the  GPR  systems  were 
incapable  of  producing  detectable  signals  from  any  of  the 
targets  deeper  than  27  inches.  Therefore,  the  antenna  system 
was  placed  over  the  center  of  the  tank,  on  the  surface,  where 
the  pipe  targets  were  a  maximum  distance  from  the  antennas. 
Additionally,  the  antenna  system  was  oriented  so  that  it  was 
cross-polarized  with  respect  to  the  pipe  targets,  thereby 
minimizing  potential  target  return  signals.  With  this 
configuration,  we  expected  that  any  target  signal  returns  would 
fall  below  the  GPR  receiver  internal  clutter  level. 

The  signal-to-clutter  measurement  required  the 
establishment  of  a  reference  signal  level.  The  reference 
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signal  chosen  was  the  return  from  a  19  inch  by  24  inch  metal 
plate  target  located  two  inches  below  the  antenna  system.  Air 
was  used  as  the  dielectric  between  the  antenna  system  and  the 
plate.  The  plate  chosen  as  the  reference  target  had  large 
enough  dimensions  that  it  intercepted  almost  all  of  the  signal 
radiated  from  the  antenna.  Thus,  a  larger  plate  would  have 
produced  an  insignificant  difference  in  the  plate  return 
amplitude . 

The  amplitude  of  the  signal  return  from  the  plate  was 
measured  and  used  as  the  reference  for  each  of  the  GPR  systems. 
The  signal  return  from  the  plate  was  within  the  dynamic  range 
of  the  GAR  radar  system  sampler  and  the  HP  8510  receiver  so 
that  there  was  no  risk  of  signal  compression.  However,  the 
exact  compression  point  of  the  GSSI  sampler  was  unknown. 

The  antenna  systems  were  then  placed  over  the  center  of 
the  soil  tank  in  the  cross-polarized  orientation.  The  receiver 
waveforms  beyond  the  surface  return  were  assumed  to  be 
dominated  by  internal  clutter  signals.  Therefore,  a 
measurement  was  made  of  the  signal  level  at  a  location 
approximately  15  nanoseconds  from  the  beginning  of  the  receiver 
range  window.  The  largest  signal  present  in  a  four  nanosecond 
region  centered  at  the  15  nanosecond  point  was  used  as  the 
representative  clutter  level  for  the  signal-to-clutter 
measurement . 

A  signal-to-clutter  ratio  was  computed  for  each  GPR  system 
and  expressed  in  decibels.  This  ratio  was  computed  by  dividing 
the  measured  amplitude  of  the  return  from  the  metal  plate  by 
the  measured  clutter  signal  level. 

The  difficulty  in  making  an  accurate  and  objective  signal- 
to-clutter  measurement  that  included  the  effects  of  the  antenna 
system  was  recognized.  Therefore,  the  signal-to-clutter 
measurement  described  above  was  used  as  a  figure  of  merit  and 
was  not  intended  to  be  used  as  a  binding  measure  of  the  clutter 
performance  of  the  candidate  GPR  systems.  Future  investigation 
may  well  reveal  a  better  method  of  measuring  the  signal-to- 
clutter  ratio  that  would  include  the  antenna  systems  and  avoid 
the  use  of  a  lossy  dielectric  as  a  matched  load  for  the 
antennas.  Nonetheless,  the  signal-to-clutter  ratios  determined 
under  this  program  (as  described  above)  were  felt  to  be 
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adequate  in  comparing  the  relative  merits  of  the  candidate  GPR 
systems . 


7.2  FULL-SCALE  GPR  FIELD  MEASUREMENTS 

Once  the  soil-model  tank  measurements  of  Phases  1  and  2 
^ere  completed,  the  corresponding  performance  of  each  of  the 
GPR  systems  was  evaluated.  GAR  engineers  concluded  that  the 
stepped-FM  system  using  the  HP  8510  gave  the  best  signal-to- 
clutter  and  dynamic  range  performance.  The  basis  for  the 
selection  of  the  stepped-FM  system  is  explained  in  more  detail 
in  Section  9.  The  stepped-FM  system  was  tested  on  a  pipe  test 
field  that  GAR  built  several  years  ago  near  its  Marietta 
facility. 

The  pipe  test  field  contains  two  plastic  pipes  and  four 
metal  pipes,  of  several  diajneters  and  buried  at  various  depths. 
Figure  7.4  is  a  top  view  of  the  pipe  test  field  showing  the 
pipe  diameters,  their  fabrication  material,  and  their  relative 
spacing.  Figure  7.5  is  a  side  view  of  the  pipe  test  field 
showing  the  burial  depths  of  the  pipes.  As  the  pipe  field  is 
several  years  old,  the  pipe  burial  depths  are  not  precisely 
known  due  to  settling  and  other  environmental  effects.  It  is 
estimated  that  the  burial  depths  are  known  to  within  +  4 
inches  . 

The  full-scale  field  measurements  were  qualitative  in 
nature.  They  were  meant  to  demonstrate  the  detection 
capabilities  of  a  stepped-FM  GPR  system.  It  was  felt  that 
quantitative  measurements  would  not  be  accurate  due  to  the 
unknown  soil  conditions  of  the  pipe  field.  The  homogeneity  of 
the  soil  could  not  be  determined,  and  the  exact  burial  depths 
of  the  pipes  were  not  known. 

The  full-scale  field  stepped-FM  GPR  measurements  utilized 
a  time-domain  synthesized  pulse  with  a  bandwidth  of  45  MHz  to 
700  MHz.  A  "normal"  bandwidth  weighting  function  was  used, 
thus  a  pulse  of  three  nanoseconds  in  duration  (as  measured  at 
the  5<2%  voltage  amplitude  points)  was  synthesized.  The 
bistatio  antenna  system  used  for  the  tests  consisted  of  two 
tr i angul ar -sheet  dipoles  with  dimensions  approximately  2.5 
times  those  of  the  antennas  used  for  the  laboratory  scale-model 
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gure 


drawing  of  pipe  test,  field  showing  pipe 
locatl ons . 


measurements.  The  antenna  system  was  a  unit  designed  to  be 
used  with  a  GSSI  short-pulse  radar  system  having  a  three  to 
four  nanosecond  duration  transmitter  output  pulse.  The  antenna 
system  was  adapted  for  use  with  the  HP  8510.  All  measurements 
of  the  time-domain  pipe  signal  returns  were  made  with  the 
antennas  polarized  such  that  the  electric  field  vector  was 
parallel  to  the  pipe  axis. 
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SECTION  8 


MEASUREMENT  RESULTS 


This  section  summarizes  the  measurement  results  obtained 
while  testing  the  candidate  G?R  source  signal  types.  The 
Gcntrolled  dielectric  environment  for  the  scale-model  tests 
provided  a  basis  by  which  the  performances  of  the  GPP.  sysoe.ms 
could  be  compared.  The  resulos  from  the  scale-model  tests  wer 
svaluc^ted  and  the  scale-model  G?R  system  wish  the  beso 
performance  was  identified. 

A  full-scale  version  of  the  scale-model  GPR  with  the  best 
osrfcrmance  was  then  fabricated  and  tested  i.n  a  realistic 
environme.nt  .  GAR  built  a  test  field  several  years  ago 
containing  buried  utility  pipes  to  be  used  for  testing  GPR 
system.s.  The  performance  of  the  full-scale  GPR  was  tested  on 
this  field. 


8.1  PHASE  1  SCALE-MODEL  GPR  MEASUREMENTS  (CLAY  SHALE) 

The  Phase  1  measurements  with  the  sc ale -mo del  GPR  systems 
used  clay  shale  as  the  dielectric  material  in  the  soil-model 
tank.  The  measured  dielectric  properties  of  a  sample  of  the 
clay  shale  are  given  in  Appendix  B.  The  diameters  of  the  pipe 
targets,  and  their  burial  depths  in  the  soil  tank,  were 
consistent  with  a  model  scale  factor  of  five.  Recall  that 
Figures  5.1  and  5.2  show  the  target  locations  withi.n  the  soil- 
mcdel  tank. 


Care  must  be  taken  when  evaluating  results  obtai.ned 
the  scale-model  systems.  When  model-soil  dielectric  p-" 
are  bei.ng  extrapolated  to  a  full-scale  envi r onrae.nt  ,  the 
conduct  ivi  cy  ,  and  hence  the  attenuatic.n  char  acter  ist  i  cs  , 
s.oale  are  scaled  downward  by  the  scale  factor.  Target 
dime.nsions,  burial  depths,  and  the  wavelengt.hs  of  the  .RF 
are  scaled  upward. 


;  n  e  r 


A  summary  of  the  Phase 
sections  which  f  o 1 1 o w .  The 
evaluations,  a  relative  die 


1  GPR  measurements  is  given  in  the 
summary  includes  s ignal - 1 c-c  1  ut t e r 
ectric  ccnsta.nt  '/erificaticn 


measurement  for  the  clay  shale  in  the  test  tank,  and  target 
detection  measurements. 


8.1.1  SIGNAL-TO-CLUTTER  MEASUREMENTS 

The  signal-to-clutter  measurement  procedure  vas  described 
in  Section  7.  The  results  are  given  in  Table  8.1  below. 


TABLE  S.l:  RESULTS  OF  SIGNAL-TO-CLUTTER  RATIO  MEASUREMENTS  FOR 
THE  CANDIDATE  G?R  SYSTEMS. 


"Stem 


mal -to-Clutter  Ratio  f dB ) 


GAR  Short  Pulse  51  .5 
GSSI  Short  Pulse  52.8 
KP  8513  FM-CW  57. '5 
HP  8510  Stepped  FM  57.5 


The  measurements  for  the  HP  8513-based  systems  may  have 
been  subject  to  dynamic  range  limitations.  The  reflection  from 
the  metal  plate  used  as  the  reference  signal  amplitude  for  che 
measurements  was  34-  d3  below  the  calibrated  network  analyzer 
reference  level  for  the  S21  parameter.  Thus,  signal-to-clutter 
ratios  of  58  d3  approach  the  dynamic  range  specification  of  93 
d3  for  the  HP  8513.  Consequently,  the  signal-to-clutter  ratios 
obtainable  with  the  HP  3513-based  syste.ms  may  be  larger  t.han 
those  presented  in  Table  8.1  above. 

The  sidelobe  levels  of  the  signal  returns  from  the  surface 
may  have  been  a  key  factor ,  if  not  the  determining  factor ,  in 
the  GPR  signal-to-clutter  measure.ments .  The  surface  return 
signal  had  sidelobes  associated  with  it  that  extended  into  the 
remainder  of  the  receiver  window.  The  sidelobe  levels  of  the 
surface  return  signal  had  the  same  relative  amplitude 
r e  1  at  1  ons.nips  as  those  measu''sd  for  the  transmitted  pulses. 

The  sidelobe  level  of  the  HP  8513  synthesized  pulses  was  down 
63  d3  at  five  nanoseconds  from  the  main  lobe  peak:  t.he  sidelobe 
level  of  t.he  GA.R  moticcycle  pulse  was  53  dB  below  the  pulse. 

The  sidelobe  level  of  t.he  GSSI  transmitter  pulse  could  net  be 
measured  because  the  transmitter  was  contained  in  a  sealed 
enclosure  . 
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The  surface  return  sidelobe  level  represents  a  practical 
-imit  for  the  minimum  clutter  level  to  be  expected  from  a  G??. 
The  other  sources  of  clutter,  such  as  connector  reflections, 
antenna  reflections,  and  transmitter  noise,  increase  the 
clutter  level  above  that  caused  by  the  surface  return  sideloc 
alone.  The  olutter  level  varies  with  location  in  the  receive 
window  and  it  is  difficult  to  distinguish  contributions  from 
various  clutter  sources. 


8.1.2  GPS  MEASUSEMENT  OF  THE  RELATIVE  DIELECTRIC  CONSTANT  OF 

THE  PHASE  1  SOIL-MODEL  TEST  TANK 

When  the  soil-model  tanx  was  filled,  the  possibility 
existed  that  the  dielectric  properties  of  the  clay  shale  in  t 
tank  would  not  match  those  obtained  in  the  laboratory 
measurements  of  an  individual  sample.  These  laboratory 
measurements  were  always  performed  on  a  clay  sample  that  was 
tightly  packed  into  a  test  cell.  Fortunately,  the  shale  was 
quite  moist  and  as  the  soil  tank  was  filled,  the  clay  was 
tamped  periodically  to  maintain  a  consistent  density. 
Nevertheless,  the  attenuation  property  of  the  shale  tank  coul 
not  be  measured  directly  as  was  done  with  the  soil  sample. 
Therefore,  an  exact  value  for  the  attenuation  of  the  shale  in 
the  tank  was  not  known. 

The  dielectric  constant  and  attenuation  properties  of  a 
given  soil  vary  with  density  just  as  they  do  with  moisture, 
the  density  increases,  so  does  the  dielectric  constant  and  th 
attenuation.  GA.R  engineers  felt  that  if  the  relative 
dielectric  constant  of  the  clay  in  the  tank  was  the  same  as 
that  of  the  laboratory  sample,  then  the  attenuation  of  the  cl 
in  the  tank  would  be  approximately  the  same  as  the  attenuatio 
measured  for  the  sample.  Therefore,  it  was  necessary  to 
measure  the  dielectric  constant  of  the  soil  in  the  test  tank. 

The  relative  dielectric  constant  of  the  clay  in  the  tank 
was  computed  using  the  signal  returns  from  the  targets.  The 
only  measured  radar  parameter  required  is  the  real-time  delay 
between  the  surface  return  and  a  target.  If  the  measurement 
performed  with  a  monostatic  antenna  system,  the  relative 
dielectric  constant  is  then  given  by: 
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e .  =  (ct/2d)^  , 


(S.l) 


where  t  is  the  real-time  separation  between  the  surface  return 
and  the  target  return,  d  is  the  known  distance  from  the  surface 
to  the  target,  and  c  is  the  speed  of  light  in  a  vacuum.  If 
Equation  (8.1 )  is  modified  to  include  the  spacing  between  the 
phase  centers  of  bistatic  antennas,  the  relative  dielectric 
constant  equation  becomes: 

£,  =  (  ct/(  2((s/2)2  +  d2)'^2  )  )2  ^  (g_2) 

where  s  is  the  spacing  between  the  antenna  phase  centers.  In 
the  case  of  the  bistatic  bow-tie  antennas,  the  phase  centers  of 
the  dipoles  were  assumed  to  be  the  feed  points. 

The  accuracy  of  this  method  of  determining  the  dielectric 
constant  of  the  clay  in  the  tank  is  affected  by  several 
factors:  the  accuracy  and  linearity  of  the  GPR  equivalent-ti.me 

receiver  range  display,  the  exact  spacing  between  the  phase 
centers  of  the  bistatic  antennas,  the  centering  of  the  antenna 
system  above  the  target,  and  the  difference  between  the 
expected  target  depth  and  its  actual  buried  depth.  The  H?  8513 
had  the  most  accurately  calibrated  range  display:  therefore, 
the  dielectric  constant  was  calculated  using  the  network 
analyzer  receiver  output  signals.  Furthermore,  the  effects  of 
errors  in  determining  s  and  d  are  minimized  by  selecting 
signals  from  deep  targets  for  computing  t.he  dielectric  ccnstani 
of  the  soil  in  the  tank.  However,  the  shale  had  high 
attenuation  properties  and  the  deepest  targets  were  not 
detectable  at  the  receiver,  so  a  relatively  shallow  target  was 
used  for  the  dielectri?  constant  measurement. 

Figure  3.1  is  a  plot  of  a  received  waveform  from  the  H? 

35  13  using  the  stepped-F?'l  t i.me-dcmain  mode.  .Marker  1  is  placed 
on  the  surface  return,  Marker  2  is  on  a  i .25  inch  diameter 
plastic  pipe  7.8  incnes  deep,  and  Marker  5  is  on  the  return 
from  a  metal  plate  2' .6  inches  below  the  surface.  Using  the 
plate  return  to  compute  t.he  dielectric  constant,  the  param.eters 
for  Equation  (8.2)  become:  d=3.5>i3  m,  s  =  3.l65  m,  t=  1  2 . 6  ns.  and 
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c=3  X  10®  m/ s .  The  relative  dielectric  constant  is  computed  to 
be  11.63. 

The  laboratory  measurement  of  the  dielectric  constant  of 
the  shale  sample  is  presented  in  Appendix  B.  For  the  frequenoy 
band  of  200  MHz  to  2.0  GHz,  the  dielectric  constant  varies 
between  14  and  16.  If  the  measurement  of  the  relative 
dielectric  constant  using  the  GPR  system  was  accurate,  the 
shale  in  the  tank  above  the  plate  target  may  not  have  been 
packed  quite  as  densely  as  the  sample  was  in  the  test  cell. 

This  means  that  the  attenuation,  in  dB/m,  for  the  clay  in  the 
oank  was  probably  less  than  that  of  the  sample. 

The  two  measurement  techniques  yielded  similar  results  for 
the  relative  dielectric  constant  of  the  shale  in  the  test  cell 
and  the  shale  in  the  test  tank.  If  the  potential  sources  of 
error  for  the  test  tank  measurement  are  taken  into  account,  the 
results  of  the  test-cell  measurement  are  an  acceptably  accurate 
representation  of  the  electromagnetic  properties  of  the  shale 
in  the  test  tank.  Therefore,  the  test-cell  parameters  were 
used  for  subsequent  calculation^  involving  the  electromagnetic 
parameters  of  the  shale  in  the  test  tank. 


8.1.3  PHASE  1  SOIL-MODEL  TANX  TARGET  MEASUREMENTS 

This  section  presents  the  target  measurement  data  for  the 
soil-model  tank  containing  the  clay  shale.  Table  8.2  below 
summarizes  the  measured  target  signal  return  amplitudes  from 
each  of  the  GPR  systems.  The  table  readily  allows  comparisons 
of  the  target-detection  performances  of  the  test  radars.  Also 
included  at  the  end  of  this  section  are  photographs  and  plots 
of  the  candidate  GPR  receiver  output  waveforms  showing  some  of 
the  target  signal  returns. 

The  left-most  five  columns  in  Table  3.2  describe  the  model 
pipe  target  environment  and  an  extrapolated  full-scale 
environment  based  on  the  model.  The  model  scale  factor  of  five 
is  chosen  for  the  computation  of  the  full-scale  environmental 
parameters,  which  is  consistent  with  the  scale  factor  used  for 
selection  of  the  model  tank  parameters.  The  left-most  column 
is  the  actual  outside  diameter,  in  inches,  of  the  pipe  targets 
buried  in  the  soil-model  tank.  The  second  column  lists  tbe 
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TABLE  8.2.  SUMMARY  OF  MEASUREMENT  RESULTS  FROM  GPR  TESTS  WITH  SOIL  MODEL  TANK  FILLED  WITH  SHALE  CLAY. 
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depths,  in  inches,  from  the  surface  to  the  tops  of  the  pipes. 
The  third  column  lists  the  equivalent  outside  diameter  of  the 
pipes  in  a  full-scale  environment.  The  fourth  column  contains 
the  extrapolated  depths  of  the  pipes,  measured  in  feet,  for  a 
full-scale  environment.  The  fifth  column  describes  t.he 
material  of  the  pipe  construction:  metal  or  plastic. 

The  remaining  four  columns  present  the  measured  amplitudes 
of  the  signals  from  the  pipe  targets  for  each  of  the  candidate 
G?R  systems.  The  cases  for  vhich  the  pipes  were  undetectable 
are  noted  in  the  table.  The  short-pulse  radar  amplitudes  are 
expressed  in  millivolts  and  the  K?  8513  target  signal 
amplitudes  are  expressed  in  micro-units.  One  micro-unit  is 
^  ,  S>(Z)i2  ,  SQiQth  of  the  ajnplitude  of  the  normalized  signal 
obtained  from  the  HP  8510  S21  a.t  calibration.  The  S2i 
calibration  is  performed  by  removing  the  antennas  from  the 
antenna  cables  and  connecting  the  cables  directly  together. 
Therefore,  the  S21  calibration  amplitude  is  much  greater  tha.n 
the  signals  from  the  antennas. 

Also  included  near  the  bottom  of  Table  8.2  is  the  shale 
clay  attenuation  expressed  in  dB/m  at  1  GHz  and  an 
extrapolation  of  this  attenuation  to  a  full-scale  environment. 
The  value  of  the  clay  attenuation  at  1  GHz  comes  from  the  grapn 
included  in  Appendix  B. 

Care  must  be  taken  when  comparing  the  target-detection 
performances  of  the  GPR  systems  using  the  amplitude  results  in 
Table  8.2.  The  aosolute  amplitudes  cannot  be  compared  between 
radars.  This  is  because  each  of  the  radar  systems  is 
constructed  differently  and  produces  different  output  levels. 
Instead,  the  comparison  should  be  based  on  the  ratios  of  the 
amplitudes  of  the  targets  as  measured  with  each  radar. 

For  example,  the  maximum  target  signal  amplitude  for  each 
GPR  is  from  the  3.88  inch  diameter  metal  pipe.  These 
.measurements  (along  with  the  corresponding  measurements  for  the 
3.38  i.nch  diameter  plastic  pipe)  ca.n  be  used  as  the  reference 
target  signal  amplitudes  for  each  radar.  The  other  target 
signal  returns  ca.n  be  expressed  as  a  fraction  of  the  amplitude 
of  t.he  cor r espondi.ng  reference  level.  These  ratios  are 
summarized  in  Table  8.3. 
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TABLE  8.3.  SUKMARY  OF  NORMALIZED  TARGET  SIGNAL  AMPLITUDES  FROM  GPR  TESTS  WITH  SOIL  MODEL  TANK  FILLED  WITH  SHALE  CLAY 
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*Denotes  undetectable 

**Denotes  detectable  only  with  strip  chart  recording.  Target  waveform  not  detectable  in  receiver  output. 


For  the  GAR  short-pulse  radar,  the  ratio  of  the  2.38  inch 
metal  pipe  return  amplitude  to  the  amplitude  of  the  0.88  inch 
metal  pipe  is  43/1120=0.038.  The  ratios  of  the  same  two  target 
returns  for  the  GSSI  system,  the  HP  8510  swept-FM  system,  and 
the  HP  8510  stepped-FM  system  are  20/700=0.029,  350/8970=0.039, 
and  314/7760=0.040,  respectively.  These  ratios  are  all  within 
3  dB  of  each  other,  as  would  be  expected.  This  is  because  all 
of  the  environmental  factors  affecting  the  signal  returns  are 
the  same  with  each  of  the  radars. 

Other  target  return  ratios  show  greater  variation  between 
the  GPR  systems.  The  ratio  of  the  0.88  inch  plastic  pipe 
return  to  the  0.88  metal  pipe  return  shows  a  1 0  dB  variation 
between  the  radars.  The  target  ratios  for  the  GSSI  system,  the 
GAR  system,  the  HP  8510  swept-FM  system,  and  the  HP  8510 
stepped  system  are  0.143,  0.277,  0.370,  and  0.440, 
respectively.  Several  reasons  may  account  for  this  variation. 

First,  the  pipe  target  was  very  close  to  the  surface, 
causing  the  target  return  and  the  surface  return  to  interfere 
somewhat.  This  may  have  led  to  an  error  in  measuring  the 
return  amplitudes  of  the  0.88  inch  pipes.  Second,  the  GSSI 
radar  used  a  different  antenna  than  the  other  three  systems, 
which  may  have  resulted  in  a  different  near-range  response. 
Third,  the  received  signal  level  is  sensitive  to  the  antenna 
coupling  to  the  ground.  Optimally,  the  antennas  should  be 
flush  with  the  soil  surface.  However,  the  soil  surface  could 
not  be  made  perfectly  flat,  so  that  some  unevenness  existed. 

Another  very  important  result  can  be  determined  from  the 
ratios  of  the  plastic  and  metal  0.88  inch  pipes.  On  the 
average,  the  plastic  pipe  return  is  about  9  dB  lower  in 
amplitude  than  the  metal  pipe  return.  This  implies  that  the 
plastic  pipe  has  a  radar  cross  section  that  is  approximately 
one-third  that  of  the  metal  pipe  in  voltage.  Similar 
variations  of  the  signal  amplitudes  are  observed  when 
comparing  the  returns  from  the  plastic  and  metal  pipes  buried 
at  other  depths.  The  implication  is  that  a  plastic  pipe  of  a 
given  diameter  must  be  closer  to  the  surface  than  a  metal  pipe 
of  the  same  diameter  in  order  to  be  detected. 

The  radar  cross  section  of  a  plastic  pipe  filled  with  a 
liquid,  such  as  water,  will  typically  be  greater  than  that  of  a 
plastic  pipe  filled  with  a  gas.  This  is  due  to  the  higher 
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relative  dielectric  constant  of  the  liquid.  The  plastic  pipes 
in  the  soil-model  tar.K  were  air-filled,  and  no  measurements 
were  made  with  liquid-filled  plastic  pipes. 

Another  phenomenon  was  noted  when  generating  strip-chart 
recordings  using  the  two  time-domain  pulsed  radars.  The  metal 
pipe  signals  were  detectable  at  much  greater  off-center  angles 
than  the  plastic  pipes.  (An  off-center  orientation  is  obtained 
by  first  placing  the  GPR  antenna  directly  over  the  pipe  and 
then  moving  the  antenna  in  a  direction  perpendicular  to  the 
longitudinal  axis  of  the  pipe.)  This  difference  in  off-center 
responses  was  greater  than  that  which  could  be  explained  simply 
based  on  the  larger  return  obtained  in  general  from  the  metal 
pipe  than  from  the  plastic  pipe.  This  means  that  the  radar 
cross  section  of  an  air-filled  plastic  pipe  varies  differently 
with  respect  to  the  off-center  angle  than  does  that  of  a  metal 
pipe.  Thus,  algorithms  to  extract  pipe  returns  from  radar 
receiver  waveforms  should  account  for  these  differences  in  the 
radar  cross  section  patterns  of  the  two  classes  of  pipes. 

A  thorough  measurement  of  the  radar  cross  section  patterns 
of  metal  and  plastic  pipes  of  similar  diameter  with  respect  to 
depth  would  be  quite  interesting.  Results  from  such  an 
investigation  would  be  valuable  to  the  GPR  designer  and  to  the 
designer  of  pipe-detection  signal  processing  algorithms.  Such 
measurements,  however,  were  beyond  the  scope  of  this  program. 

Figures  8.2  through  8.17  are  photographs  and  plots  of  the 
received  waveforms  from  each  of  the  GPR  systems.  The  antenna 
systems  were  placed  in  four  different  locations  on  the  surface 
of  the  soil-model  tank.  The  antenna  systems  were  always 
centered  directly  above  the  pipe  and  polarized  in  parallel  with 
the  pipe  axis  to  ensure  maximum  return  signal.  The  surface  and 
target  return  signals  are  identified  in  each  figure.  The 
location  of  the  antenna  systems  in  each  figure  is  described 
below. 

In  Figures  8.2  through  8.5,  the  antenna  system  was  placed 
above  the  1.31  inch  metal  pipe  buried  at  a  depth  of  7.8  inches. 
Figures  8,6  through  8.9  show  the  returns  from  the  1.25  inch 
diameter  plastic  pipe  at  a  depth  of  7.8  inches  with  the  metal 
plate  located  beneath  the  pipe  at  a  depth  of  21.6  inches.  In 
this  case,  the  antenna  was  placed  directly  above  the  pipe,  and 
the  plate  was  slightly  off-center.  Figures  8.10  through  8.13 
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Figure  8.2: 


Phase  1  measurement.  GAR  short-pulse  radar 
received  waveform  showing  the  1.31  inch  diameter 
metal  pipe  return  at  a  depth  of  7.8  inches. 


SURFACE  AT  1.2  DIV  FROM  LEFT 
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Phase  1  measurement.  GSSI  short-pulse  radar 
received  waveform  showing  the  1.31  inch  diamete 
metal  pipe  return  at  a  depth  of  7.8  inches. 


Figure  8.3: 
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Phase  1  measurement.  HP8510  swept-FM  received 
waveform  showing  the  1.31  inch  diameter  metal 
pipe  return  at  a  depth  of  7.8  inches. 


F  igure  8  .  -1  : 
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Figure  8.5:  Phase  1  measurement.  HP8510  stepped-FM  received 

waveform  showing  the  1.31  inch  diameter  metal 
pipe  return  at  a  depth  of  7.8  inches. 
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SURFACE  AT  1.5  DIV  FROM  LEFT 
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r'igure  8.6:  Phase  1  measurement.  GAR  short-pulse  radar 

received  waveform  showing  the  returns  from  the 
1.25  inch  plastic  pipe  at  7.8  inches  and  the 
metal  plate  at  21.6  inches. 
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Figure  8.7:  Phase  1  measurement.  GSSI  short-pulse  radar 

received  waveform  showing  the  return  from  the 
1.25  inch  diameter  plastic  pipe  at  7.8  inches  anc 
the  metal  plate  at  21.6  inches. 
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Figure  8.8: 


Phase  1  measurement.  HP8510  swept-FM  received 
waveform  shewing  the  return  from  the  1.25  inch 
diameter  plastic  pipe  at  7.8  inches  and  the  metal 
plate  at  21.6  inches. 
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Phase  1  measurement.  HP8510 
waveform  showing  the  return 
diameter  plastic  pipe  at  7.8 
plate  at  21.6  inches. 
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inches  and  the  meta'. 


igure  8.9: 
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Figure  8.10;  Phase  1  measurement.  GAR  short-pulse  radar 

received  waveform  showing  the  return  from  the  3.3 
inch  metal  pipe  at  18  inches. 
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Phase  1  measurement.  GSSI  short-pulse  radar 
received  waveform  showing  the  return  from  the  3.5 
inch  metal  pipe  at  18  inches. 
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.12:  Phase  1  measurement.  HP8510  swept-FM  received 

waveform  showing  the  return  from  the  3.5  inch 
metal  pipe  at  18  inches. 
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Figure  8.13:  Phase  1  measurement.  HP8510  stepped-FM  received 
waveform  showing  the  return  from  the  3.5  inch 
metal  pipe  at  18  Inches. 
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Phase  1  measurement.  GAR  short-pulse  radar 
received  waveform  showing  the  return  from  the  2 
inch  diameter  plastic  pipe  at  15.6  Inches. 
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Figure  8.15:  Phase  1  measurement.  GSSI  short-pulse  radar 

received  waveform  showing  the  return  from  the  2.5 
inch  diameter  plastic  pipe  at  15.6  inches. 
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Figure  8.16:  Phase  1  measurement.  HP8510  swept-FM  received 
waveform  showing  the  return  from  the  2.5  inch 
diameter  plastic  pipe  at  15.6  inches. 
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Figure  8.17:  Phase  1  measurement.  HP8510  stepped-FM  receive 
waveform  showing  the  return  from  the  2.5  inch 
diameter  plastic  pipe  at  15.6  inches. 


show  the  received  signals  from  the  3.5  inch  diameter  metal  pipe 
at  a  depth  of  18  inches.  Finally,  Figures  8.14-  through  8. 17 
show  the  received  waveform  from  the  2.5  inch  plastic  pipe 
buried  at  a  depth  of  15.6  inches. 


8.2  PHASE  2  SCALE-MODEL  GPR  MEASUREMENTS  (CLAY  SCHIST) 

The  Phase  2  measurements  were  performed  with  clay  schist 
in  the  soil-model  tank.  The  measured  dielectric  properties  of 
a  sample  of  the  clay  schist  are  given  in  .i^ppendix  3.  The  tank 
was  filled  using  the  same  procedure  as  was  used  with  the  clay 
shale  of  the  Phase  1  measurements.  The  targets  and  their 
placements  were  exactly  the  same  as  for  the  Phase  1 
measurements . 

The  Phase  2  measurements  consisted  of  the  verification  of 
the  soil  dielectric  constant  and  the  actual  target 
measurements.  A  s ignal - to-c lutter  measurement  was  not 
performed  because  the  schist  did  not  have  a  high  enough 
attenuation  to  make  it  suitable  as  a  matched  load  for  the 
antenna.  A  summary  of  the  Phase  2  measurements  is  given  in  the 
sections  which  follow. 


8.2.1  GPR  MEASUREMENT  OF  THE  RELATIVE  DIELECTRIC  CONSTANT  OF 

THE  PHASE  2  SOIL-MODEL  TEST  TANK 

The  measurement  procedure  to  determine  the  relative 
dielectric  ccnstaint  of  the  clay  schist  was  the  same  as  that 
used  for  the  Phase  1  measurement.  In  this  case,  however,  the 
return  from  the  metal  pipe  36  inches  below  the  surface  was  used 
as  a  reference.  Figure  8.18  shows  the  pipe  return  signal  from 
the  HP  8510  oper at ing -i-n  the  stepped  mode.  Marker  1  is  located 
on  the  surface  return  and  Marker  2  is  on  the  pipe  return.  The 
pipe  return  is  located  22  nanoseconds  from  the  surface.  Using 
Equation  (3.2)  and  s  =  0.0S2m,  one  computes  the  relative 
dielectric  constant  to  be  12.94.. 


Two  samples  of  the 
circuit  sample  cell  in 
properties  of  the  soil 
from  soil  at  the  brick 


schist  were  tested  using  the  open 
order  to  determine  the  dielectric 
in  the  tanx.  The  first  sample  was  taken 
manufacturer’s  location  and  the  second 
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Figure  8 . 18 
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Phase  2  measurement.  Plot  of  the  HP8510  returns 
of  the  surface  (Marker  1).  and  the  6.63  inch 
diameter  metal  pipe  (Marker  2). 


was  taken  from  soil  that  vas  used  in  the  test  t 
sample  produced,  a  relative  dielectric  constant 
22,  which  is  consideraoly  higher  than  the  value 
the  GPR  system.  Consequently,  the  second  soil 
measured  using  two  packing  densities,  firm  and 
effort  to  produce  a  relative  dielectric  constan 
that  was  closer  to  the  GPR  system  measurement, 
of  measurements  are  presented  in  Appendix  3. 


ank .  The  first 
of  approximate!, 
obtained  using 
sa.mple  was 
dense ,  i.n  an 
t  measure.ment 
All  three  sets 


The  third  set  of  .measurements  (dense  packing)  revealed  a 
relati'/e  dielectric  constant  varying  between  11  and  11.5  over 
the  frequency  range  of  290  MHz  to  2.3  Ghz .  However,  t.ne  curve 
does  exhibit  an  undulation  that  is  not  typical  of  such  a 
measurement.  It  was  felt  that  this  undulation  could  be 
attributed  to  the  difficulty  in  producing  a  uniformly  dense 
packing  of  the  soil  sample  in  the  test  cell.  Nevertheless,  th 
s^.:uilarity  between  the  GPR  dielectric  constant  measurement  and 
t.ne  densely  packed  test  cell  measure.ment  of  the  second  sa.mple 
implied  that  the  test  cell  measurement  of  the  atte.nuaticn  ccul 
characterize  the  soil  in  t.he  test  tank  with  a  reasonable  degre 
of  aocuracy.  Therefore,  the  matching  test  cell  measurements 
could  be  used  to  extrapolate  the  scale-model  system  to  a  full- 
scale  environment. 


8.2.2  PHASE  2  SOIL-MODEL  TANK  TAiJGET  MEASURE.MENTS 

The  Phase  2  target  measurements  were  performed  in  a 
similar  man.ner  as  t.he  Phase  1  target  measurements.  Table  3 .  i 
summarizes  the  results  of  these  target  measurements.  This 
table  is  structured  in  exactly  the  same  way  as  Table  8.2. 

Table  8.5  is  similar  to  Table  8.5  from  Section  8.1.5  in  that  i 
lists  the  normalized  received  signals  for  each  GPR,  where 
normalization  is  done . r e 1  at i ve  to  the  return  from  the 
corresponding  3 . 88  inch  pipe  (metal  or  plastic)  for  the  GPR  cf 
:.nterest.  I.n  addition,  photographs  a.nd  plots  of  the  GPR  cutpu 
waveforms  are  included  at  the  end  of  this  section. 

When  examining  Table  8.4-,  one  should  note  that  all  cf  the 
targets  were  detected  with  all  of  the  GPR  systems.  This  is 
due.  of  course,  to  the  lower  attenuation  of  the  schist  in  the 
tank.  For  the  same  reason,  target  signal  return  amplitudes 
f-om  the  Phase  2  measurements  were  in  general  s  ign  i  f  i  ca.n  t  ly 
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TABLE  8.4.  SUMMARY  OF  MEASUREMENT  RESULTS  FROM  GPR  TESTS  WITH  SOIL  MODEL  TANK  FILLED  WITH  SCHIST  CLAY. 
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TABLE  8.S.  SUMMARY  OF  NORMALIZED  TARGET  SIGNAL  AMPLITUDES  FROM  GPR  TESTS  WITH  SOIL  MODEL  TANK  FILLED  WITH  SCHIST  CLAY. 
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greater  than  comparable  signal  returns  from  the  Phase  1 
measurements.  Table  8.6  presents  the  ratios  of  target  signal 
a.mplitudes  measured  for  the  schist  clay  case  divided  by  the 
corresponding  amplitudes  measured  for  the  shale  clay  case. 

An  exception  to  the  trend  of  larger  signal  returns  with 
the  schist  clay  was  observed  with  the  GAR  short-pulse  radar 
returns  from  the  0.88  inch  diameter  metal  pipe  buried  at  4.1 
inches.  The  metal  pipe  return  was  actually  lower  in  this  cas 
by  13  percent.  Contrast  this  with  the  results  from  the 
f r equency-dcm.ain  syste.ms-  The  increase  in  return  amplitude 
frcm  the  3.88  inch  diameter  pipes  in  schist  clay  using  the 
frequency  domain  systems  was  at  least  49  percent,  and  ranged 
high  as  89  percent.  The  GSSI  short-pulse  radar  showed  an 
increase  in  the  metal  pipe  return  of  14  percent,  and  an 
increase  in  the  plastic  pipe  return  of  160  percent. 

In  general ,  the  percentage  increases  in  Phase  2  target 
return  amplitudes  relative  to  Phase  1  measurements  for  the  tw 
frequency-domain  systems  and  the  GSSI  short-pulse  radar  were 
much  higher  thaui  the  comparable  increases  observed  with  the  G 
short-pulse  radar.  These  differences  can  only  be  attributed 
the  differences  in  the  frequency  content  of  the  RF  energy  in  • 
the  pulses  generated  by  the  GAR  short-pulse  radar  as  compared 
to  those  generated  by  the  GSSI  and  the  HP  8513  systems. 


The  ratios  of  deep-target  returns  (the  pipes  with 
diameters  of  2.5  inches  or  more)  to  returns  from  the  3.38 
pipes  for  the  GAR  short-pulse  radar  were  higher  than  the 
ratios  for  the  frequency-domain  systems  (see  Table  8.5). 
ratios  of  the  deep-target  returns  to  the  0.88  inch  pipe  r 
foi  the  GSSI  radar  were  similar  to  those  of  the  GAR  radar 
However,  the  pipe  returns  measured  with  the  GAR  short  pul 
radar  did  .not  exhibit  increases  in  going  from  P.hase  i  to 
2  measur eme.nts  as  l.arge  as  those  exhibited  by  the  frequen 


sa.me 


se 


-■'has 


dcmai.n  syste.ms. 


Figures  3.19  through  8.33  are  photographs  and  plots  of 
received  waveforms  from  each  GPR  systems  with  the  antennas 
five  different  locations  on  the  surface  of  the  soil-model  c 
The  antennas  were  centered  above  each  pipe  of  interest  and 
polarized  i.n  parallel  with  the  pipe  axis  to  maximize  the  si 
returns.  The  surface  and  target  signal  returns  are  identif 
The  photcgrap.ns  are  described  briefly  in  the  paragraph  belc 
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SURFACE  AT  2.4  DIV  FROM  LEFT 

1.31  DIA  METAL  PIPE  AT  3.2  DIV  FROM  LEFT 

4  NS/DIV  REAL-TIME 


Figure  8.19;  Phase  2  measurement.  GAR  short-pulse  radar 

received  waveform  showing  the  return  from  the 
1.31  inch  metal  pipe  at  7.8  inches. 
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SURFACE  AT  1.0  DIV  FROM  LEFT 

1.31  DIA  METAL  PIPE  AT  2 . 2  DIV  FROM  LEFT 

4  NS/DIV  REAL-TIME 


Figure  8.20:  Phase  2  measurement.  GSSI  short-pulse  radar 
received  waveform  showing  the  return  from  the 
1.31  inch  metal  pipe  at  7.8  inches.  (4ns,  Div, 
real  time ) 
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Figure  8.21:  Phase  2  measurement.  HP8510  swept-FM  received 
waveform  showing  the  return  from  the  1.31  inch 
metal  pipe  at  7.8  inches. 
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Figure  8.22:  Phase  2  measurement.  KP8510  stepped-FM  receiv 
waveform  showing  the  return  from  the  1.31  inch 
metal  pipe  at  7.8  inches. 


SURFACE  AT  2.5  DIV  FROM  LEFT 
1.25  DIA  PLASTIC  PIPE  AT  3.5  DIV  FROM  LEFT 
METAL  PLATE  AT  5 . 3  DIV  FROM  LEFT 
4  NS/DIV  REAL-TIME 


Figure  8.23:  Phase  2  measurement.  GAR  short-pulse  radar 

received  waveform  showing  the  returns  from  the 
1.25  inch  diameter  plastic  pipe  at  7.8  inches  and 
the  metal  plate  at  21.6  Inches. 
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SURFACE  AT  1.0  DIV  FROM  LEFT 
1.25  DIA  PLASTIC  PIPE  AT  2 . 3  DIV  FROM  LEFT 
METAL  PLATE  AT  4 . 4  DIV  FROM  LEFT 
4  NS/DIV  REAL-TIME 


Figure  8.24:  Phase  2  measurement.  GSSI  short-pulse  radar 

received  waveform  showing  the  returns  from  the 
1.25  inch  diameter  plastic  pipe  at  7.8  inches  and 
the  metal  plate  at  21.6  inches. 
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igure  8.25;  Phase  2  measurement.  HP8510  swept-FM  received 
waveform  showing  the  returns  from  the  1.25  inch 
diameter  plastic  pipe  at  7.8  inches  and  the  meta 
plate  at  21.6  Inches. 
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Figure  8,26: 


Phase  2  measurement.  HP8510  stepped-FM  received 
waveform  showing  the  returns  from  the  1.25  inch 
diameter  plastic  pipe  at  7.8  inches  and  the  metal 
plate  at  21,6  inches. 
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Figure  8.27;  Phase  2  measurement.  GAR  short-pulse  radar 

received  waveform  showing  the  return  from  the  3.5 
inch  diameter  metal  pipe  at  18  inches. 
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Figure  8. 28:  Phase  2  measurement.  GSSI  short-pulse  radar 

received  waveform  showing  the  return  from  the  3.5 
inch  diameter  metal  pipe  at  18  inches. 
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.29;  Phase  2  measurement.  HP8510  swept-FM  received 
waveform  showing  the  return  from  the  3.5  inch 
diameter  metal  pipe  at  18  inches. 
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Figure  8.30: 


Phase  2  measurement.  HP8510  stepped-FM  received 
wavefoT'm  showing  the  return  from  the  3.5  inch 
diameter  metal  pipe  at  18  inches. 
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.31:  Phase  2  measurement.  GAR  short-pulse  radar 

received  waveform  showing  the  return  from  the  3.5 
inch  diameter  plastic  pipe  at  18  inches. 
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Figure  8.32:  Phase  2  measurement.  GSSI  short-pulse  radar 

received  waveform  showing  the  return  from  the  3.5 
inch  diameter  plastic  pipe  at  18  inches. 
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Figure  8.33:  Phase  2  measurement.  HPS510  swept-FM  received 
waveform  showing  the  return  from  the  3.5  inch 
diameter  plastic  pipe  at  18  inches. 
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Figure  8.34:  Phase  2  measurement.  HP8510  stepped-FM  receive 
waveform  showing  the  return  from  the  3.5  inch 
diameter  plastic  pipe  at  18  inches. 
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Figure  8.35:  Phase  2  measurement.  GAR  short-pulse  radar 

received  waveform  showing  the  return  from  the 
6.63  inch  diameter  metal  pipe  at  36  inches. 
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Figure  8.36:  Phase  2  measurement.  GSSI  short-pulse  radar 

received  waveform  showing  the  return  from  the 
6.63  inch  diameter  metal  pipe  at  36  inches. 
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Phase  2  measurement.  HP8510  swept-FM  received 
waveform  showing  the  return  from  the  6.63  inch 
diameter  metal  pipe  at  36  inches. 
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8.3  FULL-SCALE  GPR  .MEASUREMENTS  WITH  HP  8510 


The  full-scale  GPR  measurements  were  conducted  with  a 
stepped-FM  frequency-domain  system  as  described  in  Section  7 
above.  During  a  measurement  the  bistatic  antenna  syste.m  was 
placed  directly  above  each  pipe  target  and  a  plot  of  the 
resultant  received  waveform  was  generated.  For  the  case  of  ■ 
two  closely  spaced  copper  pipes,  the  antenna  system  was 
centered  between  them. 


Figures  8.39  through  8.43  are  the  plots  of  the  received 
waveforms  from  t.he  various  pipe  targets.  Figure  8.39  is  the 
received  waveform  when  the  antenna  was  centered  between  the  tw 
copper  pipes.  Figures  8.43  and  8.41  are  the  waveforms  from  the 
two-foot  deep  and  four-foot  deep  plastic  pipes  respectively, 
and  Figures  8.42  and  3.43  are  plots  of  the  received  wavefor.ms 
from  the  four-foot  deep  and  six-foot  deep  m.etal  pipes 
respectively.  In  each  of  the  plots,  marker  number  two  is 
placed  on  the  signal  return  from  the  target. 

I.n  Figure  3.39,  marker  .number  two  is  place  on  the  cc.mbine 
signal  retur.n  from  the  two  copper  pipes.  The  two  pipes  were  a 
an  equal  distance  from  the  antenna,  and  therefore  the  returns 
are  comfoined.  .Moving  the  antenna  away  from  the  center  point 
between  the  pipes  revealed  separate  returns  from  each,  because 
the  ranges  to  tne  two  pipes  was  different.  Thus,  t;ie  one-fooo 
spaci.ng  betwee.n  the  copper  pipes  was  greater  tha.n  the 
resolution  cell  of  the  three-nanosecond  pulse,  e.nabling  the 
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8.42:  Full-scale  field  measurement.  HP8510  stepped-FM 

received  waveform  showing  the  return  from  a  6 
inch  diameter  metal  pipe  buried  4  feet  deep. 
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returns  from  the  individual  pipes  to  be  distinguished  vhen  the 
antenna  was  moved  on  the  surface  across  the  pipes. 


Other 
of  target 
amp  1 i tude 
buried  at 
amplitude 
buried  at 


qualitative  observations  demonstrated  the  dependence 
signal  strength  on  target  radar  cross  section.  The 
of  the  return  from  the  two-inch  diameter  plastic  pipe 
four  feet  was  approximately  one-third  of  the 
of  the  return  from  the  six-inch  diameter  metal  pipe 
the  same  depth. 


The  frequency-domain  G?R  used  in  the  full-scale 
measurements  was  able  to  detect  the  six-inch  diameter  meuc::.! 
pipe  at  a  six-foot  depth.  However,  it  is  doubtful  that  she 
same  pipe  in  a  similar  soil  buried  at  ten  feet  could  be 
reliably  detected.  This  is  because  the  equations  modelli.ng  the 
metallic  cylinder  radar  cross  section  (either  Equation  4.2  or 
4. .  4 )  combined  with  the  Radar  Range  Equation,  Equation  4.1, 
predicts  that  the  target  return  signal  power  has  a  1 /R^ 
dependence  as  well  as  a  loss  term  that  is  also  exponentially 
dependent  on  R.  The  1 /R^  term  alone  predicts  a  received  signal 
power  from  a  target  at  ten  feet  to  be  21  .  6fo  of  the  power  from 
the  same  target  at  six  feet.  Since  the  attenuation 
characteristics  of  the  soil  in  the  pipe  test  field  were  not 
known  and  could  not  be  measured  for  the  entire  field,  the  exacc 
effects  of  the  attenuation  term  could  not  be  predicted. 

However,  the  combination  of  the  two  terms  mentioned  above  would 
reduce  the  amplitude  of  the  return  from  a  six-inch  diameter 
pipe  buried  at  ten  feet  to  below  the  level  of  the  clutter 
present  in  the  received  signal. 


The  full-scale  measur e.ments  using  the  stepped-FM 
frequency-domain  G?R  system  demonstrated  its  capability  to 
detect  buried  pipe  targets.  As  expected,  the  qualitative 
results  illustrated  the  dependence  of  the  target  detectability 
on  its  radar  cross  section  and  burial  depth.  The  unknown 
attenuation  characteristics  of  the  soil  in  the  pipe  test  field 
and  the  exact  burial  depths  of  the  targets  prevented  a 
correlation  between  the  scale-model  results  of  the  Phase  i  and 
Phase  2  measurements  and  the  full-scale  measurements. 
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SECTION  9 


SELECTION  OF  THE  OPTIMUM  GPR  SOURCE  SIGNAL  BASED  ON  THE  PHASE  1 

AND  PHASE  2  SOIL-MODEL  TANK  MEASUREMENTS 

This  section  contains  an  evaluation  of  the  results 
obtained  from  the  Phase  1  and  Phase  2  measurements  conducted 
with  the  scale-model  GPR  systems.  The  effectiveness  of  the 
soale-model ing  procedure  is  examined  as  veil.  Finally,  a 
selection  of  the  optimum  GPR  source  signal  ^ype  is  made,  basec 
on  the  results  from  the  scale-model  tests. 


9.1  EVALUATION  OF  THE  EFFECTIVENESS  OF  THE  SCALE-MODEL  TESTS 

The  scale-modeling  technique  used  for  this  program  to 
evaluate  the  GPR  systems  had  advantages  and  disadvantages. 
However,  GAR  engineers  feel  that  it  provided  good  results 
overall.  It  was  effective  enough  to  allow  a  confident  choice 
of  the  best  of  the  four  candidate  GPR  source  signal  types  to  be 
made,  based  on  an  assumed  set  of  design  criteria. 

The  shale  soil-model  tank  of  the  Phase  1  measurements 
provided  the  most  useful  scale-model  results.  The  loss 
characteristic  of  the  shale  was  high  enough  so  that  it  scaled 
to  a  more  realistic  full-scale  environment  than  did  the  loss 
characteristic  of  the  schist  clay.  The  loss  characteristic  of 
the  shale  in  the  soil-model  tank  was  probably  close  to  that  of 
the  clay  shale  sample  measured  with  the  test  cell.  This  can  be 
reasonably  inferred  since  the  experimental  dielectric  constant 
measurements  using  the  GPR  systems  produced  a  similar  value  to 
that  obtained  from  the  test-cell  measurements.  Thus,  a 
realistic  conversion  of  the  Phase  1  model  measurement  results 
to  a  full  scale  scenario  could  be  done. 

Several  laboratory  test-cell  measurements  were  performed 
on  samples  of  the  clay  schist  soil  in  order  to  correlate  the 
atte.nuaticn  properties  of  the  soil-model  tank  to  those  of  the 
laboratory  measurements .  One  set  of  laboratory  sample 
measurements  did  produce  results  that  corresponded  well  with 
t.ne  GPR  dielectric  constant  measurerae.nt  of  the  soil  i.n  the  test 
ta.nk.  Therefore,  it  was  felt  t.hat  the  atte.nuaticn 
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characteristic  of  the  soil  sample  could  be  extrapolated  to  a 
full-scale  environment. 

When  comparing  the  attenuation  properties  of  the  schist 
and  the  shale,  it  was  felt  that  the  shale  provided  a  more 
useful  soil  to  be  used  as  a  scale  model .  Its  higher 
attenuation  properties  allow  the  modelling  of  full-scaled 
systems  with  more  realistic  attenuation  characteristics  that 
does  the  schist. 


engineers 


were  not  able 


to  exactly  determine  the  loss 


because  the  packing  density  of  the  soil  throughout  the  model 
tank  could  not  be  determined.  Furthermore,  we  did  not  know  he 
homogeneous  the  soil  was  in  the  tank,  even  though  considerable 
efforts  were  made  to  maintain  homogeneity.  The  procedure  of 
testing  a  soil  sample  in  a  test  cell  with  the  network  analyzer 
could  not  replicate  the  exact  packing  density  of  the  soil  tank 


In  one  respect,  the  investigators  were  fortunate  that  the 
relative  dielectric  constant  of  the  schist  was  not  as  high  as 
the  value  measured  in  the  test  cell.  The  GPR  measurements  cf 
the  relative  dielectric  constajits  of  the  schist  soil-model  tan 
amd  of  the  shale  soil-model  tank  indicated  that  the  dielectric 
constants  were  very  similar  (12.94-  for  the  schist  and  11.63  fc 
the  shale).  This  meant  that  the  effective  aperture  and  gain  c 
the  triangular-sheet  dipole  antenna  were  essentially  the  same 
for  both  soil-model  tanks.  Thus,  the  radiation  bandwidth  of 
the  antenna  was  about  the  same  for  both  tanks . 


This  observation  is  important  because  it  allows  a  direct 
evaluation  of  which  GPR  system  uses  the  optimum  bandwidth 
shape,  or  equivalently,  the  optimum  pulse  shape.  If  the 
antenna  effective  aperture  and  gain  were  very  different  from 
the  shale  tank  to  the  schist  tank,  the  antenna  would  radiate 
over  a  different  bandwidth  in  one  tank  than  in  the  other.  Th 
would  bias  one  pulse  shape  over  the  other,  because  the 
different  GPR  pulse  shapes  concentrate  R?  energy  i.n  different 
frequency  bands. 


Viewed  from 
the  schist  and  s 
i.nvest  igator  s  to 
intisrss'^.  Tills 


another  perspective,  this  similarity  between 
hale  dielectric  constants  allowed  program 
better  isolate  a  single  soil  variable  of 
variable  of  i.nterest  is  the  soil  attenuation 
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vhich  must  be  scaled  by  the  scale  factor  in  the  scale-model 
environment.  The  reader  should  recall  that  in  contrast,  the 
dielectric  constant  is  not  scaled,  but  remains  the  same  in  the 
full-scale  scenario  as  in  the  scale-model  environment.  Since 
the  two  soils  differed  in  attenuation  much  more  than  in 
dielectric  constant,  one  can  safely  assume  that  many  of  the 
differences  between  the  Phase  I  and  Phase  II  results  are  due 
virtual ly  exclusively  to  this  attenuation  difference. 

3ven  though  there  was  some  uncertainty  associated  with  the 
soil  parameters  in  the  schist  soil -model  tank,  the  Phase  2 
measurements  were  quite  valuable  in  determining  an  optimum  G?R 
source  signal  type.  Those  measurements  revealed  that  the 
effects  on  the  response  of  the  frequency-domain  systems  when 
moving  from  a  high-loss  to  a  lower-loss  environment  were  much 
more  dramatic  than  the  same  effects  on  the  GAR  short-pulse 
radar.  This  was  true,  to  a  lesser  degree,  of  the  GSSI  short- 
pulse  radar  with  respect  to  the  GAR  radar.  The  implication  is 
that  the  pulse  shape  of  the  GAR  short-pulse  radar  had  an  R? 
bandwidth  that  is  less  sensitive  to  the  specific  sell 
environment  ajnd  thus  is  perhaps  better  suited  for  deep  ground 
penetration  in  varying  environments  than  the  pulse  shape 
synthesized  by  the  HP  85i0-based  systems  or  the  pulse  of  the 
GSSI  radar . 

The  ability  to  make  this  comparison  is  precisely  why  the 
pulse  widths  of  the  GAR  radar  and  the  two  frequency-domain 
systems  were  chosen  to  be  the  same.  Recall  that  an  absolute 
evaluation  and  comparison  of  the  GSSI  measurement  results  wich 
those  of  the  other  radars  is  not  possible  because  the  GSSI  unit 
had  a  different  antenna  system  auid  a  pulse  width  that  was  not 
exactly  known. 

The  R?  bandwidth  of  the  monocycle  pulse  of  the  GAR  radar 
provided  a  better  match  to  the  attenuation  characteristics  of 
t.te  soil  than  did  the  bandwidth  of  the  network  analyzer  pulse 
or  chat  of  the  GSSI  snort-pulse.  The  RF  bandwidth  of  the  other 
pulses  had  more  hig.h-f r equency  content  than  the  GAR  monccycle. 
Thus,  a  greater  percentage  of  the  energy  in  the  bandwidth  of 
t.he  HP  3513  pulse  was  attenuated  in  the  soil  than  with  the  GAR 
monocyc 1 e . 
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9.2  DETERMINATION  OF  THE  OPTIMUM  GPR  SOURCE  SIGNAL 


The  selection  of  an  optimum  GPR  system  is  not  a  simple 
task.  Some  of  the  considerations  involved  in  the  choice  vere 
discussed  earlier  in  this  report.  However,  different  RF 
bandwidths  will  be  required  for  different  applications.  A  G? 
designed  for  deep-penetration  must  radiate  low  RF  frequencies 
to  avoid  the  drastic  attenuation  characteristics  of  the  soil 
higher  frequencies.  Conversely,  a  high-resolution  GPR  must  u 
a  wider  bandwidth  than  the  deep-penetrating  GPR  system. 
Therefore,  the  optimum  GPR  must  have  selectable  RF  bandwidths 
to  provide  the  best  performance  in  varying  soil  environments. 
Based  on  the  preceding  analysis  and  the  results  from  this 
measurement  program,  a  selection  of  an  optimum  GPR  system  can 
be  made  as  follows. 


9.2.1  EVALUATION  OF  THE  GPR  SYSTEMS  TESTED 

When  deep  penetration  is  the  goal,  the  best  RF  bandwidth 
shape  for  a  GPR  will  concentrate  most  of  the  RF  energy  i.nto  t 
passband  of  the  antenna  and  the  soil.  Of  the  GPR  systems 
tested,  the  monocycle  pulse  shape,  and  its  associated  RF 
bandwidth,  appeared  to  offer  significant  advantages  in  t.his 
case.  On  the  other  hand,  if  high  resolution  is  the  goal,  the 
a  bandwidth  shape  with  relatively  higher  frequency  content  wi 
be  preferable. 

Ideally,  one  would  want  a  GPR  system  with  a  selectable 
bajidwidth  shape,  a  shape  that  accentuates  the  lower -frequency 
components  for  deeper  penetration  applications,  and  *he  hlghe 
frequency  components  for  higher  resolution  applications.  Zac 
selectable  bandwidth  shape  could  be  different  altOg^cther,  or 
the  overall  bandwidth  shape  (amplitudes  and  phases  of  the 
frequency  components)  could  be  the  same  for  all  of  the 
different  selectable  bandwidths,  with  only  the  center  frequer. 
adjusted  from  case  to  case. 

This  selectable  bandwidth  capability  is  effectively 
provided  by  the  HP  8510  in  terms  of  its  selectable  frequency- 
domai.n  weighting  functions  (see  Section  6. 3.1.  2.1).  In 
contrast,  incorporation  of  such  a  feature  in  time-domain 
systems  such  as  the  GAR  and  GSSI  GPRs  would  be  difficult.  Th 


ease  with  which  this  feature  can  be  incorporated  into  a 
f requency-domain  G?R  is  certainly  a  major  factor  favoring 
selection  of  such  a  unit  as  the  best  GPR  system  for  location  of 
underground  utilities. 

The  clutter  performance  of  the  GPR  system  is  also  an 
extremely  important  factor.  As  illustrated  in  Table  8.1,  the 
frequency-domain  GPR  systems,  based  on  the  HP  8510,  were  found 
to  have  better  clutter  performance  than  the  pulsed  GPR  systems. 
This  was  partially  due  to  the  additional  clutter  associated 
with  the  coupling  of  high-frequency  noise  from  the  short-pulse 
transmitters  nto  the  receiver  through  the  power  supplies  and 
ground  planes.  The  CW  transmitter  of  the  HP  8510  does  not 
generate  these  transients  that  produce  that  type  of  clutter 
problem.  Also,  one  has  difficulty  in  building  short-pulse 
transmitters  that  have  50  ohm  output  impedances.  Such  output 
impedances  are  required  to  terminate  any  undesired  reflections 
that  travel  back  toward  the  transmitter  from  di  scont  inui  t  i-es  or 
mismatches  in  the  RF  transmission  system. 

Higher  average  power  levels  can  be  generated  more  easily 
with  the  FM-CV  GPR  sjotems  than  with  the  short-pulse  systems. 
This  is  due  to  the  fundamental  problems  associated  with 
generating-  narrow,  high  amplitude  pulses.  Also,  low-?RF  pulsed 
radars  must  generate  a  higher  amplitude  transmitter  pulse  than 
a  high-PRF  pulsed  radar  in  order  to  achieve  the  same  average 
power  level.  As  an  example,  a  20-watt  peak  pulse  with  a  five 
percent  duty  cycle  produces  one  watt  of  average  power.  Cne- 
watt  RF  amplifiers  are  easier  to  obtain  and  less  costly  than 
20 -watt  RF  amplifiers. 

The  GPR  dynamic  range  is  another  important  consideration. 
The  frequency-domain  GPR  systems  tested  had  better  dynamic 
range  characteristics  than  the  short-pulse  radars.  Baseband 
short-pulse  radars  i.neyitably  use  broadband  samplers  on  the 
receiver  front  end.  The  samplers  are  noisy  and  have  high  noise 
figures.  T.he  frequency-domain  GPR  syste.ms  can  avoid  the  use  of 
brcadPand  samplers  on  the  receiver  front  end,  usually  resulting 
in  a  lower  receiver  noise  figure. 

The  complexity  of  the  frequency-domain  GPR  system  based  cn 
the  network  analyzer  is  much  greater  than  that  of  the  short- 
pulse  radar  systems.  This  is  true  of  almost  any  frequency- 
domain  radar,  even  if  it  is  not  based  cn  the  HP  8510.  A  trulo 
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high-performance  frequency-domain  radar  must  use  a  synthesize 
frequency  source  and  a  means  of  frequency  control  . 
Furthermore,  the  RF  circuitry  of  the  frequency-domain  G?R  is 
more  complex.  .Also,  significant  signal  processing  is  usually 
required  vith  the  frequency-domain  GPR  in  order  to  obtain 
easily  interpreted  output  information.  In  oontrast,  the  shcr 
pulse  radars  produce  easily  recorded  time-domain  output  data 
that  can  be  immediately  interpreted  to  some  degree  without 
requiring  further  processing.  However,  to  achieve  optimum 
performance,  signal  processing  must  also  be  applied  to  the 
short-pulse  radar  output. 

The  time-domain  GPR  electronic  equipment  is  smaller  and 
lighter  than  its  frequency-domain  counterpart.  Fven  though 
high-power  pulsers  for  short-pulse  radar  systems  require  larg 
power  supplies  than  the  lower-powered  frequency-domain-based 
radars,  the  frequency-domain  systems  are  still  larger  and 
heavier . 

The  short-pulse  radar  systems  are  usually  built  to  be 
quite  rugged.  The  HP  8510,  however,  could  not  be  exposed  to 
the  same  treatment  that  some  GPR  systems  receive  in  field  use 
Of  course,  a  rugged  high-performance  frequency  domain  GPR  can 
be  built.  However,  a  large,  complex,  frequency-domain  radar 
may  still  exhibit  less  field  reliability  than  comparable  shor 
pulse  radar  systems. 

Reliability  is  often  a  function  of  the  system  complexity 
Large  numbers  of  oomponents  within  a  system  usually  implies  a 
low-er  mean  time  before  failure  than  for  a  system  containing 
fewer  components.  To  minimize  reliability  problems  with  a 
f requenoy-domain  GPR  system,  the  more  oomplex  subsystems,  sue 
as  digital  signal  processing  components  and  vibr at  ion-sens i t i 
components,  may  be  permanently  located  in  a  vehicle  in  order 
minimize  handling  by  field  technicians  performing  the  GPR 
measurements.  Additional  reliability  considerations  would 
i.nclude  the  specification  of  environmentally  rugged  component 
For  a  more  complex  GPR  such  as  a  frequency-domain  radar  syste 
selection  of  rugged  ccmpcne.nts  may  be  difficult  and  expensive 

Another  factor  affecting  the  choice  of  a  GPR  technology 
system  cost.  The  frequency-domain  systems  are  usually  far  me 
expensive  than  the  short-pulse  radars.  If  FFT  and  frequency- 
domain  signal  processing  capabilities  are  added  to  short-puls 
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radars  such  as  those  tested,  their  costs  will  increase,  but 
they  will  probably  still  be  significantly  less  expensive  than  a 
high-per fornance  f r equency-domain  GPR  system. 


9.2.2  SELECTION  OF  THE  OPTIMUM  GPR  SYSTEM 

The  following  general  specifications  are  requirements  for 
a  high  performance  GPR  system  designed  to  meet  the  goals  of 
deep  penetration  and  high  resolution.  The  spec i f icat ions  ca.n 
be  applied  to  any  GP.R  technology:  f  r  eque.ncy-domain ,  or  time- 
domain  . 

The  GPR  system  for  utility  location  must  be  designed  to 
operate  with  several  RF  bandwidths.  Three  nominal  RF  bandwidth- 
ranges  are  recommended:  one  for  deep  penetration  and  lowe,- 
resolution,  one  for  medium  penetration  and  medium  resolution, 
and  one  for  shallow  penetration  and  high  resolution.  The 
suggested  bandwidths  given  below  are  meant  as  guidelines. 

We  recommend  that  for  deep-penetrating,  lower-resolution 
app 1 i cations ,  the  RF  bandwidth  should  be  confined  to  a  range 
of  50  MHz  to  200  MHz-  For  medium-resolution,  medium- 
penetration  applications ,  a  bandwidth  of  200  MHz  to  600  MHz 
should  be  used.  For  hig.h  resolution,  a  bandwidth  of  500  MHz  to 
1500  MHz  is  suggested.  These  RF  bandwidths  correspond  to 
approximate  pulse  widths  of  6.7  ns,  2.5  .ns,  and  1  ns, 
respeotively. 

A  multiple-bandwidth  GPR  system  will  probably  require  a 
separate  antenna  system  for  each  bandwidth.  Nevertheless,  for 
maximum  performance,  the  GPR  system  should  use  a  bistatic 
antenna  configuration.  This  may  present  a  problem  at  low  R? 
frequencies  because  each  of  the  antennas  might  be  quite  large. 
I.",  such  a  case,  a  monostatic  GPR  design  can  be  used,  but 
performance  will  probably  not  be  as  good  as  with  a  bistatic 
oonf igurat i on . 


9.3  THE  -MAXIMUM  PERFORMANCE  GPR 

Based  on  this  program,  G.AR  engineers  concluded  that  a 
f  r  eque.no  y-dcmai  n  GPR  based  on  a  s  t  epped- f  r  e  quency  design 
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provides  the  highest  performance  capabilities.  Such  a  syste.m 
can  be  built  to  meet  the  dynamic  range,  sensitivity,  average 
transmitter  power,  and  RF  bandwidth  shape  and  flexibility 
requirements  for  high-perf ormance  GPR  applications  more  readil 
than  its  pulsed,  time-domain  counterpart.  Furthermore,  a 
stepped-FM  GPR  utilizing  a  synthesized  frequency  source 
potentially  offers  superior  performance  over  the  swept-FM 
system  due  to  reduced  phase  noise  effects  and  frequency 
inaccuracies  (see  Section  6.3. 1.3).  However,  a  swept-FM  GPR 
can  offer  significant  improvements  in  measurement  speed  over 
the  stepped-FM  GPR. 

The  RF  bandwidth  shape  of  the  GPR  should  be  consistent 
with  that  of  a  monocycle  pulse,  or  a  similar  shape,  in  order  t 
match  the  antenna  radiation  capabilities  and  the  soil 
attenuation  characteristics.  Such  a  bandwidth  shape  can  be 
achieved  by  modulating  the  individual  transmitted  RF 
frequencies,  or  it  can  be  achieved  by  weighting  the  received  R 
frequency  components.  In  general,  weighting  the  received 
signals  is  the  simplest  technique. 

One  method  of  implementing  a  stepped-FM  GPR  would  be  to 
radiate  each  individual  R?  frequency,  while  keeping  the 
receiver  tuned  to  that  frequency  for  a  time  equal  to  at  least 
the  duration  of  the  transmitted  RF  pulse  plus  the  maximum  rang 
time.  Then  the  transmitter  and  receiver  are  stepped  to  the 
next  frequency,  repeating  the  procedure  until  all  of  the 
frequencies  have  been  radiated  and  received.  After  processing 
the  received  data  from  these  individual  frequencies,  the 
measurement  is  complete.  The  cycle  is  repeated  again  for  one 
next  measurement . 

An  important  performance  goal  in  a  stepped-FM  system  is  t 
raciate  the  individual  RF  frequency  components  with  very  pure 
specrral  contents.  ■-  .A  pure  spectral  conte.nt  is  impossible  to 
achieve,  however,  because  the  individual  transmitted  RF 
frequencies  must  be  turned  on  and  off,  or  modulated,  for  seme 
fixed  time  duration.  The  shape  of  the  pulse  used  to  modulate 
the  RF  frequency  is  typically  rectangular.  Thus,  each  of  the 
RF  frequencies  acts  as  a  "carrier"  frequency  for  the  mcdulatic 
pulse.  To  approximate  the  spectrum  of  a  CW  RF  frequency,  the 
modulation  pulse  .must  be  infinittjly  long  in  duration,  which  is 
clearly  impossible  in  a  realizable  system.  Progressively 
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shorter  modulation  pulses  broaden  the  spectrum  of  the 
individual  RF  frequency  components. 


On  the  other  hand,  the  design  engineer  wishes  to  make  the 
modulation  pulse  narrow  because  this  decreases  the  time  that 
the  GPR  must  dwell  on  each  frequency.  Conversely,  increased 
modulated  pulse  widths  decrease  the  harmonic  content  of  the 
individually  transmitted  RF  components.  A  compromise  must  be 
reached  between  the  modulating  pulse  width,  the  processing 
techniques  to  be  applied  to  the  received  waveforms,  and  the 
desired  measurement  speed. 


Recall  that  the  unambiguous  range  specification  of  the 
stepped-FM  G?R  determines  the  frequency  difference  between  the 
individual  RF  frequencies  that  must  be  used  to  reconstruct  a 
simulated  pulse.  The  specif ication  of  the  overall  RF  bandwidth 
and  the  frequency  component  spacing  determines  the  number  of 
frequency  components  required  to  span  the  bandwidth.  Thus, 
operating  at  the  minimum  possible  unambiguous  range  permits  the 
use  of  fewer  frequency  components,  which  reduces  the  GPR  syste.m 
complexity.  An  added  benefit  of  the  reduced  number  of  radiated 
frequency  components  is  a  potential  increase  in  measurement ,  or 
soan,  speed.  On  the  other  hauid,  a  reduction  in  the  minimum 
unambiguous  range  of  the  GPR  reduces  the  maximum  depth  which 
can  be  "seen"  by  the  radar.  Thus,  a  trade-off  exists  hers  as 
we  1 1  . 


Some  level  of  signal  processing  capability  is  clearly 
required  for  the  stepped-FM  radar  system.  A  minimum 
requirement  is  to  provide  frequency  bandwidth  control,  both  for 
weighting  functions  to  reduce  the  synthesized  pulse  sidelobe 
levels,  and  for  RF  bandwidth  selection  for  different 
penetration  or  resolution  applications.  Additional  signal 
processing  algorithms  should  include  FFT  and  inverse-FFT 
software  to  allow  advanced  processing  techniques  to  be  applied 
to  the  received  GPR  waveforms.  The  GPR  should  have  the 
capability  of  waveform  averaging  to  reduce  noise  effects.  The 
averaging,  of  course,  would  impact  the  GPR  measurement  speed. 
A^so,  it  may  be  desirable  to  provide  a  means  of  calibration  to 
remove  undesired  effects  of  cables  and  connectors,  similar  to 
one  calibration  capabilities  of  the  HP  85')^. 

The  GPR  will  require  measure.ment  data  rates,  or  "scan" 
rates  cf  at  least  '3  Hz,  and  should  probably  include  the 
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capability  of  a  53  Hz  scan  rate.  These  data  scan  rates  will  b 
a  challenge  for  a  synthesized,  stepped-FM  design  to  accornplish 
due  to  the  synthesizer  switching  speed  limitations.  However  a 
swept-FM,  or  a  stepped-FM  system  that  does  not  phase  lock  to 
the  individual  frequency  components  as  it  switches,  can  operat 
at  even  higher  measurement  speeds. 

The  design  and  construction  of  the  GPR  receiver  and 
transmitter  circuits  must  be  done  with  great  care  to  maximize 
the  signal-to-clutter  and  signal-to-noi se  performance.  The 
nu.mbsr  of  RF  connectors  in  critical  transmission  lines  should 
be  minimized  to  avoid  parasitic  reflections.  .\11  RF  oompone.no 
should  provide  good  impedanc-  matches  to  the  transmission 
lines,  or  they  will  also  produce  reflections  at  their  input  an 
output  ports  that  will  contribute  to  the  internal  system 
clutter . 

The  selection  of  the  transmitter  power  level  is  dependent 
on  safety  requirements,  electromagnetic  interference 
regulations,  the  depth  penetration  specification,  and  the  RF 
bandwidt.h  of  the  GPR.  Safety  regulations  and  EMI  requirements 
provide  the  upper  limit  on  transmitter  power.  Given  this 
constraint,  the  GPR  transmitter  should  operate  with  a  power 
level  that  is  high  enough  to  produce  a  recognizable  target 
signal  level  at  the  receiver  from  a  specified  target  geometry, 
at  a  specified  burial  depth,  with  a  specified  soil  attenuation 
Hig.her  transm' tter  powers  will  act  to  increase  the  target 
signal  dynamic  range  above  the  receiver  noise  floor,  assumi.ng 
t.hat  the  receiver  does  not  compress  the  high  level  returns, 
riowever  ,  recall  that  a  GPR  system  is  often  c  1  ut  ter- 1  imi  t  ed  , 
cherefore  increasing  the  transmitter  power  may  be  ineffective 
i.n  helping  to  detect  low-amplitude  targets. 

The  RF  receiver  should  be  designed  to  have  a  low  noise 
figure.  This  will  probably  require  a  low-noise  amplifier  on 
the  receiver  front  end.  Such  an  amplifier  will  reduce  the 
effaces  of  a  noisy  mixer  or  sampler  following  it  in  the  RF 
signal  path.  Iz  may  also  help  to  increase  the  receiver 
se.ns :  t  i  Vi  ty  .  Although  it  is  simple,  and  can  provide  high 
performa.nce  i.n  some  respects,  using  a  sampler  on  the  receiver 
front  e.nd  should  be  avoided  due  to  its  noise  performance  a.nd 
limited  dynamic  ra.nge  capability. 
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The  receiver  must  utilize  I  and  Q  (in-phase  and 
quadrature-phase)  detection  to  provide  a  vector  representation 
of  the  received  waveform  so  that  proper  phase  relationships 
between  the  RF  frequency  components  can  be  maintained.  This  is 
a  stringent  requirement,  because  phase  and  amplitude 
distortions,  especially  if  only  a  few  RF  frequency  components 
are  used,  will  result  in  errors  when  reconstructing  the  pulse. 
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APPENDIX  A 


Measured  Electromagnetic  Properties  of  Various  Soil  Types 


A-1 


GEORGIA  RED  CLAY 


A  large  sample  of  Georgia  red  clay  (Cobb  County,  GA)  was  baked  to  remove 
the  natural  water.  Six  smaller  samples  were  prepared  from  the  dry  clay  and  deionized 
distilled  water  with  approximately  the  following  water  contents  (percent  water  by  dry 
weight)  0,  2,  5,  10,  20  and  30%.  The  measured  values  of  percent  water  by  weight, 
percent  water  by  dry  weight,  and  density  are  in  Table  2. 

The  measured  electrical  parameters  e/,  e/*,  cr.  and  a  are  presented  as  a 
function  of  frequency,  50  MHz  <  f  <  1.25  GHz,  in  the  two  attached  figures.  The  data 
for  the  six  samples  are  plotted  on  the  same  graphs  to  show  the  effects  of  the  water 
content.  These  graphs  clearly  demonstrate  the  increases  in  the  dielectric  constant  e . 
and  effective  conductivity  that  occur  with  increasing  water  content.  The 

attenuations  for  these  samples  are  in  most  cases  much  smaller  than  those  for  the 
Japanese  soils. 


RED  CLAY 
(Cobb  County,  GA) 


A13ER 

PERCENT  WATER 

BY  WEIGHT 

PERCENT  WATER 

BY  DRY  WATER 

DENSIT' 

sir.,  err. 

1 

0 

^  0 

1.-6 

2.0 

2.0 

I. -9 

3 

4.6 

4.8 

1.55 

- 

8.7 

9.6 

1.66 

5 

15.9 

18.9 

1.72 

6 

22.3 

28.7 

i  .  2 
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Attenuation  (dB/m) 
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J  0.05 
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o 
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0.00 


0  L-J 
0.00 


SAND/WATER  MEASUREMENTS 


The  electrical  parameters  of  five  sand/water  mixtures  were  measured,  and  the 
results  are  recorded  in  the  accompanying  graphs.  The  parameters  were  measured  over 
the  frequency  range  50  MHz  <  f  <  3  GHz. 

The  "dry"  sand  was  mixed  with  distilled  water  to  obtain  various  percentages  of 
water  by  weight.  The  highest  water  content  was  obtained  by  completely  saturating  the 
sand.  The  parameters  for  the  mixtures  are  given  in  the  accompanying  table. 


.Note  that  there  are  ripples  on  the  data  for  the  higher  water  content  mixtures. 
The  ripples  are  the  result  of  nonuniform  packing  of  the  sand  in  the  measurement  cell. 


SAND-WATER  MIXTURES 


Sample 

No. 

Percent  Water 
bv  Weight 

Percent  Water 
by  dry 
Weight 

Density 

gm/cm^ 

Nominal 

Er 

1 

- 

- 

1.37 

2.3 

(Dry) 

2 

9.1 

10.0 

1.38 

5 

3 

11.7 

13.2 

1.52 

10 

4 

16.7 

20.0 

1.57 

14 

5 

21.3 

27.1 

1.92 

26 

(Saturated) 


Attenuation  (dD/m)  Conductivity  (S/m) 


JAPANESE  SOILS 


Samples  from  each  of  the  seven  Japanese  soils  were  oven  dried  and  used 
to  determine  the  water  content.  The  results  of  these  measurements  (percent 
water  by  weight,  percent  water  by  dry  weight,  and  density)  are  summarized  in 
Table  1. 

A  brief  physical  description  of  each  sample  follows: 


Sample  1: 
Sample  2: 
Sample  3: 
Sample  4: 
Sample  5: 
Sample  6: 

Sample  7; 


redish  brown,  soil  plus  clay,  fine  particles, 
dark  brown  to  black,  sandy  soil. 

dark  tan  to  gray,  dry  soil,  several  small  rocks  in  sample, 
mostly  large  soft  rocks,  some  dark  tan  soil, 
dark  tan  to  black,  coarse  sandy  soil. 

dark  gray  clay,  very  cohesive,  very  fine  particles,  very 
homogenous. 

black  clay  with  some  sand  and  rocks,  very  cohesive. 


The  complex  relative  permittivity  measured  for  each 

sample,  and  the  effective  conductivity  and  attenuation  constant  a  (d3/m) 

were  computed  from  c^.  The  results  are  presented  as  a  function  of  frequencv. 
50  MHz  ^  f  ^  1.25  GHz,  in  the  attached  graphs. 

The  average  dielectric  constants  (real  part  of  the  relative  permittivitv, 
^')  are  graphed  as  a  function  of  the  water  content  of  the  samples  below. 
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NUMBER 

PERCENT  WATER 

BY  WEIGHT 

PERCENT  WATER 

BY  DRY  WEIGHT 

DENSIT 
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1 

48 

91 

1.32 

2 

16 

19 

1.83 

3 

19 

23 

1.69 

4 

26 

35 

1.68 

5 

11 

13 

1.86 

6 

40 

66 

1.53 

7 

20 

25 

1.36 
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BICKERSTAFF  CLAY  MEASUREMENTS 


The  electrical  parameters  of  two  types  of  clay  obtained  from  the  Bickerstaff 
rick  Co.  were  measured,  and  the  results  are  recorded  in  the  accompanying  graphs, 
i  he  parameters  were  measured  over  the  frequency  range  of  0  -  2.0  GHz.  The  first 
two  sets  of  clay  samples  were  taken  at  Bickerstaff  and  included  a  sample  of  clay 
schist  and  a  sample  of  clay  shale. 


The  first  two  sets  of  samples  had  the  following  properties: 


Density 

Percent  Water  to  Dry  Weight 

Schist 

1.S8  gm/ml 

21 

Shale 

1.92  gm/ml 

14 

Two  additional  sets  of  solid  electric  property  measurements  were  made  using  a 
sample  of  schist  from  the  test  tank.  The  two  sets  of  measurements  were  performer: 
with  different  packing  densities  of  the  samples  in  the  test  cell.  These  measurement 
results  appear  beginning  on  pages  B-9. 
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I  111.'  Ill’  85i0  neiuork  .iiialy/u!  s\stem  is  ,iii  a<.l\ iiiiL'ci.1  jjhI  sophisi icaicd  measur¬ 
ing  instiumoiu  dcsignc(.l  U)  make  mieiowasc  measuiemems  of  many  kinds.  But 
the  basic  pi  incijilcs  ol‘  iis  operation  are  I'aii  ly  simple.  'I’he  information  in  this  par: 
t)!'  the  l  ip  S.^li)  system  mtmual  is  designed  to  help  you  get  the  most  fiom  yotu' 
Ml*  iS5lO  system  by  esphtining  some  of  the  btisie  principles  of  its  operation  and 
the  etiiiipmeni  that  sliould  be  iiseil  ssiih  it.  .•\ciu;il  metisurements  are  described  in 
the  nest  section  of  this  mamuil.  as  tin  Introductory  .Vleasurement  Sequence. 

In  the  present  secticin.  the  IIP  8510  network  antily/.er  system  is  described  and  a 
(\picai  metisuremeni  is  explainetl  in  terms  of  a  system  block  diagrant.  Digital 
micropi oeessing  of  the  dal.i.  sources  eomptiiible  with  the  HP  8510  system,  and 
tile  IIP  8510  sssiem  test  sets  are  also  described  using  block  diagrams. 

Pstiemely  ttcctiraie  and  comities  metisurements  are  possible  with  the  HP  85 10 
system,  and  for  this  reason  ticee.ssories  such  as  cables,  attenuators,  e.xiension  lines, 
adapters.  ;tnd  calibrtition  .and  verification  kits  arc  unusually  important.  .Acces¬ 
sories  u  hich  should  be  used  with  the  HP  8510  system  arc  listed  and  discussed 
aftei  the  system,  its  sources,  and  its  test  sets,  have  been  descia’bed. 
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Operating  and  Programming 


BASIC  PRIiNCIPLES 

Vector  network  analyzers  such  as  tiic  HP  8510  netwoik  analyzer  system  measure 
ilie  magnitude  and  phase  characteristics  of  linear  networks  such  as  filters,  am¬ 
plifiers.  attenuators,  and  atitennas.  As  with  all  network  analyzers,  two  kinds  of 
measurements  are  made:  reflection  measurements  and  transtnissioti  measurements. 
An  incident  signal  generated  by  an  RF  source  is  compared  with  the  signal  trans¬ 
mitted  through  the  de\  ice  or  reflected  from  its  input. 


TRANSMISSION  REFLECTION 


’.NCICENT  ,  TEST 

Transmission  measurements  are  made  by  comparing  the  transmitted  signal  to  the 
incident  signal.  This  results  in  measuretnent  data  on  transmission  characteristics 
of  the  netwoi  k  such  as: 

Insei  tion  Loss  or  Gain. 

Transmission  Coefficient. 

Electrical  Dehiy. 

from  which  Electrical  Length  can  be  obtained. 

Deviation  from  Linear  Pliase, 

Group  Delay. 

[Reflection  measurements  are  made  by  comparing  the  reflected  signal  to  the  inci¬ 
dent  signal.  This  results  in  measurement  data  on  reflection  characteristics  of  tlie 
des  ice  such  as; 

Return  Loss. 

Standing  Wave  Ratio  (SW'R), 

Reflection  Coefficient. 

Impediince. 

Mathematical  anal\sis-of  transmission  and  leflection  data  on  the  swept  response 
of  the  network  also  makes  it  possible  to  deieiinine  the  position  ;md  magnitude  of 
impedance  ch.mges  'sith  respect  to  a  reference  plane.  This  analysis,  called  time 
domain  .imilssis.  is  done  using  Fourier  Transform  principles  and  is  possioie  on  HP 
8510  network  an.il'.  zer  s>siems  ct|uipped  with  lime  domain  Option  OlO, 
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HP  85  10  NPTNVORK  A^AL^  ZER  SVSTE.M 

Tlic  MP  S5I0  network  analyzer  system  lias  four  essential  parts: 

a  source, 
a  tost  sot. 

a  signal  detector  and  analog-to-digit;il  coinerier.  and 
a  digital  microprocessor  and  display. 


The  source  prosides  the  RF  signal.  The  test  set  separates  this  signal  into  an  inci- 
(jont  signtii  sent  to  the  des  ice-under-test  and  a  refei'ence  signal  against  whicii  the 
transmitted  and  reflected  signals  aie  later  compaied.  It  also  receives  ti'ansmitted 
and  reflected  signals  from  the  device-under-test.  The  signal  detector  and  analog- 
to-digital  converter  takes  all  of  these  signals  and  converts  them  to  digital  informa¬ 
tion  for  high-speed  processing.  The  digital  microprocessor  controls  the  system, 
analyzes  the  digitized  signals,  collects  errors,  and  displays  the  results  in  a  variety 
of  formats. 

In  the  HP  8510  netwoik  analyzer  svstem.  these  essential  parts  are  individual  HP 
instruments  configured  together  make  up  the  HP  8510  system: 

HP  8,'4\-series  synthesized  sweeper, 
or 

HP  835\-serics  ssseep  oscillator  with  an  appropriate 
HP  835x\-serics  plug-in. 

HP  851  \-series  test  set; 

HP  85102a  if  detector; 

HP  85  101. A  display  processor. 

Additional  system  components  can  include  hardcopy  output  devices  such  as  a 
printer  and  or  a  plotter,  and  an  HP  series  200  computer  serving  as  an  e.xternal 
controller  for  programmed  operation. 
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OR 


HP835X-SER1ES 
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b\S[  i;M  111. OCR  DIAGRAM 

\  siinpiiricLi  '^ioek  diaai'am  of  [he  MI^  8510  network  analyzer  system  is  shown  in 
higurc  a 

-\s  rigine  a  shous.  ilie  HP  8510  network  analyzer  is  a  high  perfoiniance  vcctoi' 
rccciser  wiili  four  inputs,  two  independent  measurement  clianneis,  and  an  internal 
micioeompu:e;  to  automate  measurement  and  data  processing  operations.  .A 
soecial  -System  Dus  provides  fast  digital  communication  between  the  instruments 
dial  make  laa  die  system,  allowing  the  network  analyzer  to  make  full  use  of  the 
source  and  :esi  sec  capabilities.  This  interface  also  presides  direct  data  transfer 
:o  the  iiarucop;.  des  ice  for  neat,  permanent  records  of  die  measurement  displas'. 

During  a  c  nicai  measurement  with  the  source  operating  in  die  I'anip  sweep 
mode,  the  souice  is  swept  from  the  lower  to  the  higher  measurement  frequency  in 
a  linear  ramp.  .Signal  separation  components  in  the  test  set  apply  a  portion  of'the 
incideiK  signal  and  the  responses  from  the  device  under  test  to  the  first  frequenev 
conversion  stage. 

Digital  comnumiention  between  the  receiver  and  the  test  set  pretunes  the  65  to 
.'^00  MHz  \  oitage-tuned  local  oscillator  (VTO)  so  that  one  of  its  harmonics  mixes 
with  die  stiniuius  to  produce  a  first  IF  frequency  close  to  20  .VIHz.  Fine  tuning  is 
accomplished  b\  comparing  the  IF  frequency  with  tlie  internal  20  MHz  crystal 
reference  and  sweeping  the  local  oscillator  to  track  tiie  stimulus  frequency. 

W  hen  [he  'oc.al  oscillator  reaches  its  upper  frequency  limit,  the  sweep  is  stopped, 
the  local  oscillator  Is  retuned,  phase  lock  is  reestablished,  and  the  sweep  is  con¬ 
tinued.  Since  the  first  local  oscillator  frequency  is  selected  algorithmically  from 
the  known  stimulus  freriuency,  the  measurement  is  free  from  harmonic  skip. 

The  second  fiequency  conversion  produces  an  IF  frequency  of  100  kHz  for 
application  to  die  detection  and  data  processing  elements  of  the  receiser.  Be¬ 
cause  the  frequency  conversions  are  phase  coherent  and  tlie  IF  signal  paths  are 
carefully  matciled.  magnitude  and  phase  relationships  between  the  input  signals 
are  maintained  tlirougiiout  the  frequency  conversion  and  detection  steps.  .Auto¬ 
matic.  fully  caiiiirated  autoranging  IF  gain  steps  iiinintain  the  IF  signal  at  opti¬ 
mum  le\  els  tor  detection  o\er  .i  wide  dynamic  range. 
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Tile  leierence  detector  ciiaiinel  'can  use  either  input  a[  or  a-^  as  the  reference 
signal.  Tile  test  detector  can  use  any  of  the  inputs  as  the  test  "signal.  During  the 
s  Aeep.  tiie  selected  inputs  are  sampled  up  to  401  times,  with  sample  timing  ac¬ 
complished  Oy  setrsdig  the  0  to  10  volt  sweep  output  fi  om  the  source.  W  ith  401 
points  selected,  at  each  positise  0.025  volt  change  in  the  sweep  \oltnge  all  selected 
inputs  are  sampled  and  applied  to  the  leference  and  test  synchronous  detectois. 

The  S'.  I’chronous  sleteciors  dereiop  the  real  '.\i  and  imaginarv  {\)  parts  ot'  the 
>ignar  Tile  \.h  pans  aie  seuuentially  coiueiied  to  digital  sallies  and  lead  P\  tne 
C  end  ai  Processing  Lint  iC'PL  ).  Then  digital  techniques  are  used  that  practical!'; 
chuiiiuiie  drift,  offsets,  and  cicrulariiy  eriois  as  sources  ol'  nieasui eiiient 
uiicei  t.iiiu  . 
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POST-DETECTION  DIGITAL  SIGNAL  PROCESSING 

Post-detection  digital  signal  processing  (Figure  4)  proceeds  under  control  of  die 
CPL .  a  microprocessor  etiuipped  with  256  kBytes  of  RAM.  256  kBytes  of  mag¬ 
netic  bubble  memory,  and  26  kBytes  of  ROM. 

The  CPU  takes  advantage  of  multi-tasking  software  architecture  and  several 
distributed  processors  to  provide  a  very  fast  display  update  rate.  It  accepts  the 
digitized  real  and  imaginary  data  and  corrects  gain  and  quadrature  errors  before 
the  rel'erence  and  test  pairs  are  ratioed  and  stored  in  the  raw  data  array.  If 
.net aging  is  on.  the  incoming  data  is  averaged  with  the  existing  data  as  it  is 
stored. 

While  the  data  acquisition  software  is  continually  filling  the  raw  data  array,  the 
data  processing  soft'^are  is  processing  the  data  for  the  two  independent  display 
channels. 

if  error  correction  is  turned  on.  tne  raw  data  and  error  coefficients  from  the 
selected  calibration  coefficient  set  are  used  in  , appropriate  computations  by  a 
deilicated  vector  math  processor.  Next,  phase  offsets  commanded  by  the  electri¬ 
cal  delay  and  leference  plane  extension  are  added  to  the  data.  If  a  time  domain 
presentation  is  selected,  the  corrected  data  is  consented  from  the  frequency  do¬ 
main  to  the  time  domain  using  the  inserse  Fourier  Chirp  Z  transform  technitiue 
and  stoied  into  the  corrected  data  arrays. 

The  (tiemoi  >'  a/  /a\  s  arc  filled  from  the  corrected  data  array  under  control  of  the 
'.iser  svith  trace  data  for  use  in  vector  computations  with  the  current-  corrected 
data.  If  trace  math  is  selected,  vector  multiplication,  division,  addition,  or  subtrac¬ 
tion  is  performed.  The  resulting  data  are  formatted  according  to  the  FORMAT 
selection,  point-to-point  smoothing  is  applied,  if  selected,  and  stored  into  the  for¬ 
matted  data  arrays.  The  traces  are  now  scaled,  and  output  to  the  display  memo- 
' v  here  the  trace  data  is  combined  with  vaiious  CRT  annotation  data.  A  dedi¬ 
cated  display  processor  asynchionously  converts  the  formatted  data  and  annota¬ 
tions  for  disphy  at  a  flicker-free  rate  on  the  vector-writing  CRT. 

When  the  operating  system  detects  a  front  panel  button  push,  it  executes  the  com¬ 
mand  immediately  (as  when  a  paiameter  change  is  made),  or  it  makes  the  se¬ 
lected  function  the  active  function  and  awaits  input  from  the  knob,  numeric  pad. 
01  .STEP  keys  (as  when  there  is  a  scale,  division  change),  or  it  presents  a  softkey 
menu.  .Selecting  some  functions  aborts  the  data  processing  operation.  For  exam¬ 
ple.  MEASl  RE.MEN  I  RESTART  restarts  all  measurement  related  functions  to 
the  beginning  of  the  dtita  acquisition  group  (a  group  is  that  number  of  sweeps 
needed  to  make  the  metisurement  completely;  how  may  sweeps  are  taken  tlius 
depends  on  the  metisurement);  PRESET  initializes  the  system  to  a  pre-defined 
state. 
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SOI  UC  ES 

The  RE  soiiicc  in  an  111*  S5IU  network  analyzer  system  ctin  be  cither  an  EIP 
,S.1J\-sei  ics  s\  ntliesizeil  sweeper  e'l  111*  835\-seiies  sweep  oscillator  with  an  EIP 
8 '"'w-series  pliig-in  (Eigiires  5  and  0).  I'liese  sources  have  the  correct  rtnaiog 
interface  signals  and  liill  compaiiiiility  with  the  digittil  8510  System  Btis.  If  an 
HP  STox-series  sweep  oscillatoi  is  used,  both  the  sweep  oscillator  and  the  plug-in 
may  need  to  be  retiol'itled  with  certain  later  revisions  of  the  firmware  to  be 
ccmpatible  with  the  HI*  8510  ssoem.  Consult  your  Hcw'lett-Packard  representa¬ 
tive  if  you  need  more  information  on  compatibility  (|uesiions. 

The  8510  system  bus  allows  the  network  analyzei  to  act  as  the  system  coturoller 
by  managing  the  source  using  standaid  llP-il3  protocol.  Capabilities  added  by 
the  sssiem  bus  include  alternate  sweep,  in  which  a  different  freqtiency  range 
may  be  selected  for  each  measuiement  channel,  and  control  of  necessary  sotirce 
functions  using  the  HP  8510  front  ptinel  controls. 

Both  tspes  of  sotirces  can  operate  in  the  Ramp  Sweep  mode,  in  which  the  net¬ 
work  analyzer  ditecis  the  source  to  sweep  in  a  linear  ramp  over  the  selected 
freiiuenc\  range.  I-IP  8-'v4\-series  mstrumeiits  provide  better  performance  in  this 
Ramp  Sweep  mode  thtin  do  IIP  8.55\-serics  instruments,  because  of  the  'Lock- 
and-Roli"  tuning  tcchnic|tie  useil  ri  the  III*  854\  series.  In  this  "Lock  and  Roll" 
technique,  the  first  freciuencv  ol'  'te  sweep  is  set  with  synthesizer  accuracs  tind  a 
linear  analog  sweep  proceeds  to  P  e  stop  frequencs'.  E(.)r  sweep  whdths  less  than  5 
MHz.  fully  locked  synthesizer  pe:  ormance  is  obitiined  over  the  complete  sweep. 
Instruments  in  ilic  Ml*  8.’5\  series  re  v'jveii-loop  't'lCi-itined  sources. 

The  1-lP  8.540-\  can  tilso  opei.itc  in  the  .Stcis  Sweep  mode.  In  this  mode,  syn- 
thesizer-ckiss  fretiuency  .icctirtic;-  and  lepeattibilitv  is  obitiined  by  phase-locking 
(he  source  ;n  each  of  the  up  tv)  -"1  frev|uency  steps  over  the  selected  frcciuency 
rtintio.  This  mode  pros  ivies  the  h  .itest  accuracs  tililivnigh  tit  a  reduced  met’sure- 
Mteni  speevi 
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111’  8510  NehNork  Analyzer 


TLS  r  SETS 

Tiitj  Ml^  S5i\-seiics  test  seii  in  the  i-!l’  8510  network  analyzer  system  litue  three 
main  lanctions.  The}  provide: 

the  input  output  ports  to  connect  die  device-under-test; 
signal  separation  to  separate  the  reference  and  test  signals; 
and  RF  to  20  MHz  conversion. 


'Tile  HP  85!  I. A  frecitiency  con\erter  differs  sliglitly  in  tliat  it  does  not  have  sig¬ 
nal  separation  des  ices,  thus  allowing  custom  configurations.)  The  frec|uency  con- 
\  erter  is  fully  integrated  into  the  signal  separation  path  to  provide  optimum  per- 
foimance.  Taking  the  test-to-reference-signal  ratio  in  S-parameter  test  sets  after 
electronic  switching  eliminates  signal  path  selection  repeatability  errors. 
Partimeter  selection  is  controlled  from  the  network  analyzer  front  panel. 


TLihU’  1  HP  ii5 1 x-Scrics  Te^i  Scis 


Test  Set  .Model 
Nuinber.  Ty  pe 

Test  III  (lilt  Port 
Connector 

Frequency 

Range 

HP  85 1  1  A  Fi'ecitiencv 

3.5mm  (f) 

0.045  -  26.5  GHz 

Con\  erter 

1 

HP  85 12.A  Reflection. 

7mm 

0.500  -  18.0  GHz 

Trtinsmission 

HP  85  i  .'.A  Reflection,, 

3.5mm  (m) 

0.045  -  26.5  GHz 

Transmis.sion 

HP  85I4a  S-Para meter 

7mm 

0.500  -  18.0  GHz 

HP  851  5.A  .S-Para meter 

3.5mm  (m) 

0.045  -  26.5  GHz  | 

1 

NOTE  -  HP  8512  and  HP  85(4  test  sets  are  usable  to  0.045 
GHz.  although  with  degraded  performtince  specifications. 
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Rcrk'dioii/'I  |■;Ill^mis^i()ll  Ics(  Sols.  1  lie  111’  <M5I2A  ;ind  IIP  8513/\  icl'kelioii/ 
1  liiiismissioii  lest  sens  (I'iiiuic  7)  |■)i\)vi<.le  ,iti(i)in;:!ic  ' elec! ion  vd'  S|  |  or  St|.  FuIIv 
ei  I  oi  ■con  ceierl  incasui  c'lnciUs  Tor  oiic-port  devices  can  l;e  made  using  iTie  1-Poi  i 
calibraiioii  piocediire.  Tlie  compreliensis e  One-Path  2-l*oii  caiibi.iiion  pioccdure 
pi  ovii.les  I'lill  erroi'  correction  Tor  tvvo-pori  devices  /)'  ilie  dc'.  ice-imder-tcst  is 
manually  reversei.1.  'I'lie  IIP  <S5I2A  test  set  must  use  ii  20  dB  attenuator  at  the 
device  end  ol‘  the  transmission  return  cable;  the  I  II*  851  3A  lest  set  must  use  a  10 
tIB  aitenualoi'  at  the  de\  ice  end  of  the  transmission  return  cable. 
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S-Parainetor  lest  Sets.  'I'iie  8514A  and  111’  851 5A  S-paiaiiieter  test  sets 
iFigiire  S)  pio^ide  atitoinatic  selection  oT  S|  |.  St].  Sp.  and  .St"’.  i  ii'-  'tinuiius  is 
aiitomatic;il!;.  switched  foi'  foiwai'd  ;ind  rcvei^e  nic^isureinenTs,  allowing  fully 
eiioi’  collected  nieasurements  for  one-poit  devices  and  foi'  two-port  desices 
w  ithout  the  need  manually  to  reverse  the  dc\  ice-iinder-tcst.  Bias  input  and  sense 
connections  are  provided  to  tfllow  testing  aciisc  devices.  Internal  0  to  90  dB.  10 
dB  step  attenuators  are  provided  to  control  the  incident  siiimilus  level  at  the 
de\'ice-undei  -test  input,  without  causing  a  change  in  the  reference  signal  level. 
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Operating  and  Programming 


Custom  l  est  Sets.  'I  o  conl'iguro  signal  separation  of  your  own  design,  use  the 
1 11^  K5I  I  I'letiuency  eoiuciiei  fl'i<Mire  9).  IT  yotii'  test  setup  does  not  follow 
the  eons  emions  of  ilie  rel'Icclion  1 1 .insmission  or  S-pnraineter  test  set,  use  the 
lUiOrriNU  P.ARANU' rUR  seciuenee  t)!'  the  HP  8510  system  to  select  appropri- 
nie  icl'erciice  and  test  inputs  to  be  used  foi'  the  measurement. 


ar  INPUT 

4S  SHl _ 20  XHl 
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ri;si  I’ORi  lu.n  KN  c  aiii.i;.s,  a  i'i  i.m  a  ior^ 

lliuli  MUiiliiy  calilcs.  aiiciutaiois.  ai!a(Mcis.  aiul  oilier  acecssoiics  aie  essomial  if 
one  is  U)  aeltiese  aceuraic.  rc|X-aiablc  mcaMiieiiicnis.  Worn  or  unsiablc  cables 
anti  coiineciors  will  increase  nieasuremeni  errors  due  lo  tlireclivii y.  misnintch,  and 
frei|uene\  lesponse  effeeis.  Cheek  eabies  and  conneeiors  rei’iilarly  and  replace 
:heni  'Ahenewei'  neeessai  y. 

1  esi  poll  reuirn  cables  useti  wiih  .in  (-11’  8510  netuoik  analyzer  sysieni  imisi  be 
tlur.ibic  and  siable.  anti  care  Is  retiuireti  to  atoiti  damaging  lliem.  Cables  can  be 
desuovcti  In  escessite  (less  than  5-incli  ratlins)  beiitls.  liven  with  careful  use. 
caoies  do  weai-  oui  eveniually.  and  for  ihis  reason  all  cables  shoultl  be  treaied  as 
consumable  iienis  lo  be  replaced  as  often  as  necesstiry.  The  most  important 
ch.irticiei  isiie  of  all  cables  is  mitiinuim  iiKignittide  anti  ph.'ise  change  between 
mosements  (flevures)  of  the  etibie.  Replace  a  ctible  when  large  magnitude  and  or 
phtise  chtinges  occur  when  the  cable  is  moved. 

The  cables  recommendetl  below,  in  good  condition,  must  be  used  for  detailed 
pei  foi  intmee  verification  of  the  UP  8510  systeih.  These  ctible  sets  have  low  in¬ 
sertion  loss,  good  elecirictil  mtiicii.  and'high  return  loss,  and  they  are  stable  in  use. 
Foi  other  .li^plicaiions.  .my  high  i|uality  cable  set  can  be  used. 

Recommentletl  eabies  .ind  (when  retiuireti)  20  or  10  tlB  tiiienuators  for  the  test  set 
confi'iur.iiions  thtii  c;m  be  used  in  tin  I  IP  8510  network  analyzer  system  are: 


/T'/’/t'  2.  Cu'/’/t’.'  u7/</  .■i/h'/iiu/lor); 


!  I  esl  Set 

i  Model  Number 

1 

Test  Port 
Connector 

Return  Cables, 
.Allemialors 

1 - 

!  HP8512.A  ! 

t  ,  i 

/  mm  ! 

HP  85i32.A.  8492.A-0:0 

i  HP  851. MX  1 

.15  mm  tm)  i 

HP  85I.M  A.  849.1C-0I0 

1  HP8514A 

7  mm 

HP  85I.72B  (2  in  set) 

i  HP  851 5  A  1 

1 

.15  mm  tm)  j 

HP  85 1  3 1  B  (2  in  set) 
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EXTENSION  LINES 

Exieinal  rcrerence-iignal-ixilh  exiension  lines  on  ihe  lest  set  rear  panels  are  used 
lo  balance  ihe  rererence  and  test  signal  path  lengths  according  to  the  port  1  and 
poi  t  2  connections  to  the  test  device.  These  extension  lines  (and  the  signal  paths 
they  apply  to)  depend  on  the  test  set  and  are  as  follows.  The  standard  lengths 
described  in  the  next  seseral  paragraphs  balance  the  cable  configurations  already 
listed  in  Table  2. 


TEST  SET 

LABEL 

SIGNAL  P.ATH 

Rcflcclion- 

Transmission  Test  Sets 

HP  8512A 

EXTENSION  A 

bl.  bn 

E.XTE.NSION  B 

HP  8513A 

E.XTENSION  A 

‘'^1 

E.XTENSION  B 

bp  b2 

S-Parameter  Test  Sets 

HP  8514A 

E.XTENSION  A 

E.XTENSION  B 

a2 

HP  85I5A 

EXTENSION  A 

•''1 

EXTENSION  B 

a-- 

Reflection/Transmission  Test  Sets:  HP  85  12, A,  HP  85I3A.  When  using  a  stan¬ 
dard  test  setup  (HP  8.1I.'2.A  or  HP  851,^I,A  cable  and  attenuator)  with  a  desice- 
under-test  connected  directly  to  Port  I,  use  the  short  extensions.  HP  part  number 
085 1  2-200 lb.  On  these  test  sets,  one  of  the  lines  is  in  the  test  signal  path,  and  this 
fact  makes  it  possible  to  add  bias  tees,  step  or  fixed  attenuators,  amplifiers,  isola¬ 
tors,  or  other  de\  ices. 

S-Parameier  Test  Sets:  HP  85  14A.  HP  8515.A.  When  using  a  standard  test 
setup  svith  the  de\  ice-under-test  connected  at  the  ends  of  the  HP  851  31B  or  HP 
851.32B  test  port  return  cables,  use  the  long  extensions.  HP  part  number  08414- 
20013.  When  connecting  the  device-under-test  directly  at  Port  1  and  using  a 
sntgle  Hl^  85I.HA  or  HI*  85132,A  cable,  use  the  short  extension  lines.  HP  part 
number  085 1  2-200 1 9.  -  - 

E' tension  Lines  ma>  be  chtinged  to  other  lengths  of  high  qu;ilit>’  cable  '  low  mser- 
;ion  loss,  high  reiuin  loss,  stable  in  use)  m  order  to  balance  electrical  lengths  m 
other  conl'iguiations.  .Signal  path  balance  is  less  important  when  using  the  HP 
,8 ''40  s.nthesized  sweepei.  particularly  in  the  Step  sweep  mode. 
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Al) Al’  l  KKS 

ll'  .uiaptcis  must  tic  usLt.1  to  connect  the  ticviccs  unclcf  test,  n.so  only  high-<.|ualit;/ 
m.luiTiei s  sucli  us  tliose  suppllei,!  in  the  III*  85052.A  (3.5niin)  ttnd  the  Ml*  85054A 
(1  yiK'-N)  caiibi'tition  kits.  Keep  the  mating  siiifaees  eletin.  inspect  all  connectois 
sistuilly  hefoi'e  every  use.  tintl  use  eonneelor  gtiges  to  ecrify  that  the  mating 
loleianees  tire  witliih  speeil'ietiiions.  Always  use  ;i  tort|ue  wrench,  set  to  the  cor¬ 
rect  icsrt|ue.  when  lightening  or  remo\  ing  connections. 

fesi  sets  which  have  ,v5mm  eonnectois  on  the  test  ports  (e.g.  MI*  851 3A,  HP 
X515A)  c;in  he  used  with  lest  port  return  cables  which  have  7nini  connectors  i)y 
using  the  adaiMers  in  ihe  Ml*  851  30A  specitil  3.5mm-to-7mm  tidapter  set.  Tliese 
uuapieis  provide  a  rugged  interface  for  attaching  the  7mm  test  port  return  cauies. 

For  best  results,  these  Ml*  85I30.A  adapters  (not  tlie  ;id;ipters  in  the  HP  85052.A 
3.5mm  ctilibraiion  kit)  should  be  used  if  7mm  calibration  or  verification  devices 
.li  e  Lised  for  c.ilibi  tition  or  perfomiance  verification  of  ;i  3.5min  test  set  (e.g.  HP 
85I3.\  or  HP  8515.A).  The  tidapters  in  the  calibration  kit  are  suitable  only  in  the 
opposite  case,  w  hen  3.5mm  devices  tire  used  with'a  7mm  test  set. 
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CALIBR ATION  KITS 

L  sc  only  ilic  iiighcst  quality  calibration  standards;  devices  which  bas  e  a  known 
response  and  aic  stable  in  use.  Only  if  tlie  calibration  devices  used  ha\e  an  ac¬ 
curacy  eciual  to  or  gieater  than  those  in  the  HP  85050A  (7nini)  and  HP  S5052A 
(aOinni)  calibration  kits  will  they  provide  the  calibration  and  error  correction 
accuracy  needed  to  acliieve  ftill,  specified  measurement  accuracy  with  the  HP 
8510  network  analyzer  system. 

Also  be  awaie  that  calibration  standards,  like  all  devices,  can  become  worn  and 
unsttible  with  use.  When  a  calibration  device  is  no  longer  stable  and  repeatable, 
or  shows  signs  of  connector  damage  or  wear,  it  must  be  replaced.  Detailed  han¬ 
dling  and  storage  instructions  appear  in  the  calibration  kit  operating  and  set",  ice 
manuals. 

Characteristics  for  the  standards  in  the  HP  85050A  (7mm),  HP  85052.A  (3.5mm) 
and  HP  85054.A  (Type-.N)  calibration  kits  are  loaded  from  the  tape  cartridge 
supplied  with  the  calibration  kits.  Characteristics  can  also  be  defined  by  the  user. 
Each  calibration  kit  is  supplied  with  a  data  cartridge  on  which  is  stored  the 
nominal  characteristics  for  each  of  the  calibration  devices  in  the  kit. 

The  HP  85050A  "min  calibration  kit  consists  of  open  and  short  circuit  termina¬ 
tions.  fixed  and  sliding  loads,  a  7nim  connector  gage,  gage  calibi'ation  block  and 
aligning  pin.  extra  precision  6-slot  center  collets,  a  center  collet  extractor,  a  "mm 
connector  torc|tie  wrench,  and  the  device  data  cartridge.  Option  010  adds  a  30 
cm  beadless  airline,  w  hiCh  is  used  for  time  domain  applications. 

The  HP  85052A  3,5mm  calibration  kit  consists  of  male  and  female  open  and 
short  circuit  terminations,  fixed  and  sliding  loads.  7mm-to-3.5mm  adapters, 
matched  3.5mm-to-3.5mm  adapters,  a  .3.5mm  connector  torque  wrench.  3.5mm 
connector  gages  and  gage  calibration  block,  and  the  device  data  cai’tridge.  Op¬ 
tion  OlO  adds  a  15cm  bead'ess  airline,  which  is  used  for  time  domain  applications. 

■file  HP  85054A  Type-\  Cilibration  Kit  consists  of  male  and  female  Type-N 
open  and  short  circuit  terminations,  fixed  and  sliding  loads.  7mm-to-Type-N 
adapters,  and  the  de\ ice  data  caitridge  . 

When  other  calibration  kits  are  used,  nominal  characteristics  of  the  standards  can 
be  defined  by  the  user  from  the  front  panel  of  the  HP  8510,  tising  the  MODlFh 
CAL  KIT  set|uence  desecibed  in  .Vicasuremem  Calibration  part  of  this  manual. 
After  the  calibration  kit  standards  tire  defined,  the  data  can  be  recorded  on  taoe 
tlien  loaded  from  tape  wheiieser  reriuired. 
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\  PRIFICATION  KITS 

Pei  ‘'omi. nice  \  ei  it’ication  btandai  cls  are  used  to  determine  tliat  the  system  can  be 
calibrated  and  pi  odtice  good  measurement  results.  De\  ices  in  tlie  verification  kits 
are  precision  de\  ices  which  slioiild  be  treated  vvitli  care  and  used  only  in  specific 
siiLiations.  not  on  a  day-to-day  basis.  These  devices  ha\e  been  ciiaracterized  on  a 
standards-class  network  analyzer  by  e.xperienced  factory  personnel.  If  you  use 
proper  calibration  and  measurement  techniques,  your  measurement  results  should 
be  comparable  to  the  data  supplied  with  the  devices,  within  the  system 
specifications. 

Onl>  erification  des  ices  which  hav  e  an  accuracy  equal  to  or  greater  than  those 
in  tile  HP  85051. A  (~mm)  and  HP  85053A  (3.5mm)  verification  kits  can  be  used 
to  verify  HP  8510  network  analyzer  system  specifications. 

The  HP  85051. A  ("mm)  and  HP  85953A  (3.5mm)  verification  kits  both  include 
fi.xed  attenuators  (20  dB  and  50  dB  for  7mm.  20  dB  and  40  dB  for  3.5mm)  and 
beadless  and  stepped  two-port  airline  mismatch  standards.  Data  for  the  devices 
includes  a  device  data  sheet  which  lists  fully  error-corrected  data  and  measure¬ 
ment  uncertainty  data  on  all  devices  in  the  kit  at  various  specified  frequencies. 
This  measurement  uncertainty  includes  both  the  uncertainty  of  the  HP  factory 
measurement  system  and  the  specified  uncertainty  of  the  user's  system. 

The  device  data  sheet  with  the  HP  8505 1. A  7mm  verification  kit  lists  data  at  20 
freciuencies.  19  of  them  within  the  specified  range  of  the  HP  8512A  and  HP 
85 14 A  test  sets.  The  device  datti  sheet  with  the  HP  85053,A  3.5mm  verification 
kit  lists  dtita  at  18  frequencies.  The  data  cartridge  contains  formatted  trace  data 
on  the  devices  before  they  were  shipped  from  the  factory,  as  measured  on  a 
standards-class  '(P  8510  network  analyzer  system.  The  formatted  trace  data 
contains  information  for  201  fi'equencies. 

To  veiify  svstem  performance  using  these  standaids.  perform  standard  7mm  or 
3.5mni  two-port  calibiation  and  tiieasurement  pioceduies.  present  corrected 
response  of  standard  device,  then  read  the  marker  at  the  soecified  frequency 
points  .and  compare  our  measuied  data  witli  the  standard  data  supplied  with  the 
devices.  Pel  Cl'  to  the  Performance  Tests  section  of  the  HP  8510  system  manual 
for  detailed  s;. stem  performance  verification  instructions. 
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INTRODlCTlOiN 

This  part  of  the  MP  8510  network  tinnlyzer  system  mtintial  explains  how  to  make 
icl'.ection  and  transmission  measnremenis  in  the  time  domain.  Measurements 
actualh  made  in  the  I'retinency  domain  are  transfoimied  mathematically  into  the 
time  domain  using  the  internal  higli-speed  computer  in  the  HP  8510.  and  this 
icquires  lhat  the  system  be  equipped  with  Time  Domain  Option  010.  either  at  the 
nine  o('  original  shipment  or  by  means  of  the  HP  85012.A  Time  Doimiin  Software 
Package. 

ri'.e  nine  doimiin  band  pass  mode  is  especially  useful  for  measuring  band-iimited 
de'.  ces  and  in  making  fault  location  measurements.  The  rime  domain  low  pass 
mode  simulates  the  traditional  TDR  measiuOiment  and  makes  it  possible  to  deter¬ 
mine  the  type  of  discontinuity  piesent  in  a  device.  Both  modes  are  explained 
here,  as  are  special  time  domain  featuies  such  as  masking,  windowing,  and  gating. 
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GENERAL  THEOR\ 

TIic  '  el.Hionship  hei\'.cen  the  Fi'etjiiency  Domain  response  and  the  Time  Domain 
I  espouse  ot'  a  netssoik  is  described  by  the  Fourier  Transform; 

FREQUENCY  DOMAIN  TIME  DOMAIN 

HlD  - >  h(t) 

It  is  theiefore  possible  to  measure  the  response  of  a  device  under  test  (DUT'i  in 
tlie  Fietiiienc;  Domain  and  then  mathematically  calculate  the  inserse  Fourier 
Tiansform  of  the  data  to  gise  the  Time  Domain  response.  The  internal  high¬ 
speed  computer  in  the  FIP  8510  does  iliis  calculation  using  Chirp-Z  Fast  Fourier 
Transform  computation  techniques.  The  resulting  measurement  is  the  fully  error- 
corrected  Time  Domain  reflection  or  transmission  response  of  the  device  dis¬ 
played  in  near  real  time. 

In  Figure  "0.  the  Frequency  and  Time  Domain  responses  of  the  same  des  ice  aie 
Uisplaved.  The  Frequency  Domain  refledtion  measurement  is  a  composite 
response  of  all  of  tlie  discontinuities  present  in  the  device  under  test. 

Tile  Time  Domam  measurement  sliovvs  the  effect  of  each  individual  discontinuity 
as  a  function  of  time  (or  distance).  The  time  domain  response  shows  that  the 
device  response  consists  of  three  separate  impedance  changes,  with  the  second 
discontinuif V  ha',  mg  a  relfection  coelTicient  magnitude  of  0.013.  This  discontinu- 
itv  IS  located  lo",5  picoseconds  from  the  reference  plane  relative  to  the  speed  of 
light  in  a  \acuum.  iln  the  time  domain  trace  shosvn  in  Figure  "6.  the  display  and 
tlie  marker  show  the  round-trip  time  to  the  reflection  and  back:  335  ps.) 
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I  IML  DOMAIN  MODLS 

I  lie  111*  8510  neluork  analy/er  s\\tcin  lias  two  (Jirrerent  inodes  ol'  operation  Tor 
l  ime  Domain  measm emenis,  Uand  Pass  and  Low  Pass. 

I  he  Hand  Pass  intHlc.  llic  most  general  purpose  mode  of  operation,  gis'es  the  Im¬ 
pulse  lespstnse  sil"  ihe  desiee.  Hand  l*ass  will  work  with  any  device  and  over  any 
rrei|uene\  range  and  is  the  letist  complicated  mode  to  use. 

I  he  I  ou  Pass  mode  is  used  to  simulate  the  traditional  Time  Domain  Reflec- 
lomeier  (1  DU)  measurement.  The  response  gises  the  user  information  to  deter¬ 
mine  the  t\pe  of  discontinuity  present  (U.  L.  or  C).  The  Low  Pass  mode  will  also 
pros  isle  either  the  impulse  or  step  response  of  the  device. 
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lIMt  DOMAIN  BAND  PASS 

1  he  IkificJ  Pass  nioilo  is  so  named  l>ecansc  it  vvill  work  with  band-limited  de\  ices. 
1  his  is  a  distinct  tidsantage  over  traditional  TDR.  wiiich  retinires  that  the  DLT 
he  tible  to  oiiertJie  down  to  dc.  With  Band  Pass  tliere  are  no  restrictions  on  the 
rreqiiency  range  of  the  measurement. 


Reflection  ^leasll^eIncn^s  Using  Band  Pass 

nelore  making  Time  Domain  rericetion  measurements,  it  is  necessary  to  perform 
the  appiopiiaie  metisnremenl  calibration. 

•  Press  PRtSET. 

•  Perform  an  Sj  |  l-PORT  calibration. 

Loiive  the  sliding  load  connected  ;ind  obsei  ve  the  Frequency  Domain 
response  as  the  sliding  element  is  moved. 

•  Press  DOMAIN.  TIME  B.\ND  PASS. 

•  Ihess  VU'IO  to  di'.plav  the  tiace  /tnd  observe  the  Time  Domain 
response  tis  the  sliding  element  is  moved. 

I  he  tvp'L.ii  I  retiuency  Domain  tind  fime  Domain  lesponses  of  a  sliding  load  are 
vhovvn  111  F'lgure  "T. 


I  rClIIICIKV  DnillUIII 


rinic  Diiiihiin  Band  Pn.w 


/■/cn/'f  McaMircnicni  (>l  a  Slii/iiii:  Lnad 
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Move  the  sliding  element  and  observe  the  response  in  both  the  Fret|iiency  Do¬ 
main  .111(1  the  Band  Pass  Time  Domain.  The  Freciuency  Domain  measurement  of 
the  sliding  load  should  change  \ery  little  when  the  slide  is  moved  (unless  the  cal¬ 
ibration  is  bad).  However,  the  Time  Domain  measurement  shows  the  individual 
response  of  the  load  element,  and  it  moves  along  the  horizontal  axis  as  the  slide  is 
moved. 


interpreting  the  Band  Pass  Response  Horizontal  Axis.  In  Band  Pass  reflection 
mcasuicments.  the  horizontal  axis  represents  the  amount  of  time  that  it  takes  for 
an  Miipulse.  launched  at  (lie  test  port,  to  reach  the  discontinuity  and  return.  Thus, 
this  is  the  two-way  travel  time  to  the  discontinuity,  which  in  Figure  77  is  the  load 
clement  of  the  sliding  load. 

The  Marker  reads  out  both  the  time  (x2)  and  the  electrical  length  (\2)  to  the  dis¬ 
continuity.  The  electrical  length  is  obtained  by  multiplying  the  time  by  the 
velocity  of  light  in  a  vacuum  (2.997925E8  m/sec).  To  get  the  physical  length, 
multiply  the  electrical  length  by  the  relative  velocity  of  light  in  the  transmission 
medium. 

In  the  Time  Domain,  the  STI.MliLl'S  keys  (START.  STOP.  CENTER,  and 
SP.A.N)  refer  to  time,  and  tliey  can  be  used  to  change  the  horizontal  (time)  axis  of 
the  display  independent  of  the  frequency  range  chosen.  This  can  be  done  using 
the  knob,  step  keys,  or  the  keypad.  The  keypad  terminators  also  refer  to  time  in 
seconds  (with  the  lowercase  prefixes). 


Interpreting  t.hc  Band  Pa^s  Response  Vertical  .Axis.  The  quantity  displayed  on 
the  \ertical  axis  depends  on  the  format  selected.  Band  Pass  is  PRESET  to  the 
Linear  Magnitude  format  which  displays  the  response  in  reflection  coefficient  (p) 
units.  This  can  be  thought  of  as  an  average  reflection  coefficient  of  (he  disconti¬ 
nuity  o\  er  the  freciuency  lange  of  the  measurement. 

Othei'  useful  formats  are  listed  in  Table  13.  The  Band  Pass  response  gi\es  the 
magnitude  of  the  reflection  only  and  has  no  impedance  information  (R.  L.  or  C). 
This  information  is  axailable.  howevec.  in  the  Low  Pass  response. 


Table  IJ.  Lbclul  Tinic  Dantatn  Band  Fa^''  Formal^ 


FOR. MAT 

PARAMETER 

LINEAR  MAG 
LOG  M  AG 
.SVVR 

REFLECTION  COEFFICIENT  UNITS 
RETURN  LOSS  (dB) 

SWR  units 
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I  anil  Location  Measurements  I  sing  Band  Bass 

1  lie  Uand  Pass  inoile  is  \  ei  y  usel'id  in  making  l  aiili  location  inc.isiii eincnis.  Fig- 
iiie  "S  sliovvs  ilie  Band  l^ass  lime  Domain  mcasiiicmcnt  of  a  length  of  coasial 
cable  basing  itmiliple  discontinuities  and  terminaieil  in  50  ohms.  \otc  the  lespon- 
SC'>  Oi  Cii  eh  discontinuity  atid  of  the  terminating  element. 


S  ”  ;  -  1  “St  Ae 

aE?  ^  •  « 

»  c39  32  _ 


I't'^urc  “b.  C\tNc  Fault  UKatum  M caMirciucni  Band  Pa.'^s 


Also.  i7ecause  the  Band  Pass  motle  will  work  o\er  any  rie()uency  range,  it  can  be 
Lisesl  to  do  fault  location  in  band-limited  transmission  media,  such  its  waveguide. 
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Traii'diiission  MeasuioiiKMils  in  Panel  Pass 

The  Band  Pass  mode  is  also  usefui  in  imikiiig  transmission  measuiements.  Bcfoie 
making  Time  Domain  iinnsmission  niensnrements,  it  is  neccssais  to  perform  tlie 
appropi'iate  measurement  calibration. 

•  Press  PRESET. 

•  Perform  an  S7  i  RESPONSE,  FULL  2-PORT,  or  ONE  PATH  2- 
PORT  calibrafion. 

•  Connect  a  20  dB  coaxial  attenuator  and  observe  the  Frecinency 
Doniain  1  espouse. 

•  Press  DO.MAIN.  Tl.ME  BAND  P.ASS. 

•  Press  .A I  TO  to  display  the  trace. 

The  Fiequency  Domain  and  Time  Domain  responses  of  a  20  dB  attenuator  are 
shown  in  Figure  79. 
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Operating  and  Programming 


Inlcrproting  ll\e  IJand  Pass  Transmission  Response  Horizontal  Axis.  In  Time 
Domain  [lansmission  measiiiemcnts.  tlie  hoi'izontal  axis  is  displayed  in  units  of 
(imo.  The  lesponse  ol'  the  thru  connection  used  in  the  calibration  is  an  impulse  at 
t  =  0  and  with  unit  height,  indicating  that  the  impulse  made  it  through  in  zero 
lime  ;ind  '^vith  no  loss.  When  a  de\  ice  is  inserted,  the  time  axis  indicates  the  prop¬ 
agation  delay  or  electrical  length  of  the  device.  Note  that  in  Time  Domain  trans¬ 
mission  measurements,  the  saliie  displayed  is  the  actual  electrical  length  (not  \2). 
The  Mtirker  reads  out  the  electrical  length  in  both  time  and  distance.  You  must 
multiply  the  dist.mce  number  by  the  relative  velocity  of  the  transmission  medium 
to  get  tile  actual  pnysical  length. 


Inlcriireting  the  Band  Pass  Transmission  Response  Vertical  .Axis.  The  vertical 
axis  displays  the  titinsmission  response  in  transmission  coefficient  units  (t)  in  the 
Linear  .Magnitude  format  and  the  transmission  loss  or  gain  in  dB  in  the  Log 
Magnitude  format.  This  can  be  thought  of  as  an  average  of  the  transmission 
response  o^er  the  fre(|uenc\'  range  of  the  measurement.  For  the  20  dB  attenuator 
example,  the  Band  Pass  lesponse  has  a  magnitude  of  0.10  transmission  coefficient 
units  (-20  dB  insertion  loss). 
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1  IMl.  DOM  \1N  LOW  1'  \SS 

1 'K  I  ovv  I’.iss  niDclc  (0  lime  Domiiiii  is  iiscsl  lo  siimil.nic  (lie  tinijiiion.nl  TDR 
.iie.isiireniem.  1  his  nuHle  ui\cs  the  user  iiil'oi  imiliDii  lo  dcicrminc  (lie  type  of 
Jiseoe:  i;uiii  N  ( R.  ! ..  oi  C)  that  is  present,  l  ow  Pass  provides  ilie  best  resolution 
ii.istest  rise  lime).  ;ind  it  may  be  used  lo  give  either  the  Step  or  Impulse  response 
oi  a  des  lee. 

1  he  Low  Pass  morie  is  less  general  purpose  than  ILiikI  Pass  in  that  it  places  strict 
iimiMiions  on  the  I'reiiaency  range  of  the  measurement.  It  reciuires  that  the 
1  leiiuenc)  Domain  rlata  points  l>e  h.irinonieally  related  from  dc  to  SI  OP 
:  e(|uenc>  (S  I  OP  =  \  \  Sl  AR'l'.  where  =  NLMUL'R  of  POINTS).  The  dc 
lemieiic;.  lesiK'iise  is  esitapoiaied  rrom  the  low  rrei|uency  data.  The  reciuiie- 
:;:e;’.!  to  p.iss  dc  is  the  same  limitation  that  exists  for  trtidiiional  TDR 
mc.isui  emenis. 


Celling  I  reiiueiic,(  iiange  lor  lime  Domain  Lon  Pass 

bcto'c  m. iking  me.ismeme:ns  in  '.ite  1  ow  Pass  mode,  the  rietiuency  range  ol'  the 
■■e.ism  emena  must  be  set  so  that  .S  TOP  =  n  \  S  TART,  ss  here  n  is  the  niimbet'  oi' 
V'"'w  !  h's  can  iie  '.lone  directK  b\  the  user,  or  else  it  wall  be  done  automatic.''i;l;. 

.  :icn  '.n.e  sL'I  TRLQ.  (LO^^  PASI^i  sotiae;  is  pressed.  This  key  is  included  in 
'  c  C  \L.  Menu  and  .loo  al'iei  the  I  IML  LO'^  PA.Sb  sol'tkes.  Because  the  Hi* 
' .1  a.  ;ii  ”01  con.ert  ;o  the  Lo'a  Pass  mode  until  the  SL  T  I  RLQ.  (LOW  PASS) 
■  c  is  Oi  csssd  le.ist  once,  it  is  \  er\  import. mi  that  this  be  done  before  caiib;  ai- 
c.  O;  '.c: '.s  ISC.  going  to  Lv)w  P.iss  will  ch.tnge  the  measurement  freciuencies 
I  'A  ;ii  on  n  01  :  e:  :  Ol  correction. 

i'  c'vi'_  .SLT  I'RLO-  ^L()^^  PA.'^.S)  will  set  the  STOP  Treciuencx  as  close  ns  pos- 
'.iiic  1(1  the  (.due  enieted  b\  the  user,  .ind  it  will  set  the  START  freiiuenc;.  eciuai 
slO!'’  N  M  .11)  cN.impie.  iT  the  user  selects  101  points,  with  START  =  lOO 
Ml!.:,  .ind  MOP  =  'c  oMz.  theiT  pressing  SLT  I'RLQ.  (LOW  PASS)  w;il 
.  , 'I  ;  \R  i  to  nu.d  M!  17.  .=  M  OP  lUI). 

:;c>,:,;sc  ii’c  iowesi  mc.."'U' emetti  ■■e(|uenc>  Tor  the  IIP  8510  is  43  MTlz.  Toi'  e.Tcn 
,;.i'c  o!  \  thcic  ;s  niiiinuim  .laowabic  S’l  Oi*  l  ie(iue:ic\  that  can  be  used,  .ina 

'  _ c '  m  A  Miiz.  I.mie  14  desciaPes  the  minimum  I'renuen.c;-  ini'.ge 

...  ,  us‘.'..  ;  ('■  c.ic'i  '..line  oT  N  when  m. iking  Low  Pass  Time  Dom, 

, ^ .  e m c : ti  s. 
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Fahlc  14.  Miiunniiii  Fraiuciu  v  Fur  Tiii:c  Dniikiiii  Lon-  Pu.'^\ 


MiMBER  of  POINTS  (N) 

MINFMIM  FREQUENCY  RANGE 

51 

45  .MHz  to  2.295  GHz 

101 

45  .MHz  to  4.545  GHz 

201 

45  MHz  to  9.045  GHz 

401 

45  MHz  to  18.045  GHz 

NOTE;  If  the  sotiice  cannot  opeiate  ovei  the  rec|uired 
fietiuency  lange.  the  HP  8510  will  nevertheless  attempt 

the  operation. 

11'  the  STOP  f'leciiiLMicy  eiiieiecl  is  lower  than  the  miniimini  that  is  available  for 
the  value  of  N  selected,  then  pressing  the  SET  FREQ.  (LO^V  P.ASS)  softkey  will 
ciiange  the  STOP  freciuency  to  that  mininutin  value.  For  example,  if  Number  of 
Points  =  201.  start  =  100  \tl-lz.  and  STOP  =  6.00  GHz.  then  pressing  SET 
FREQ.  (LONN  PASS)  will  change  START  to  45  .MHz  and  STOP  to  9.045  GHz  (= 
S  I  . ART  \  201).  Because  of  these  restrictions  on  the  fiec|iiency  range  of  the  mea- 
suiement.  the  Low  Pass  mode  is  most  useful  for  measuring  lowpass  broad  band 
dev  ices. 


Analyzing  Low  Pass  Reflections 

As  mentioned,  the  Low  Pass  mode  gives  the  TDR  response  of  the  device  under 
test.  This  response  contains  information  that  is  useful  in  determining  the  type  of 
discontinuity  piescnt.  Before  making  actual  measurements  in  the  Low  Pass  mode, 
it  is  helpful  to  review  the  Low  Pass  responses  of  known  discontinuities.  Each 
circuit  element  of  Figuie  80  was  simulated  to  show  the  corresponding  Low  Pass 
Time  Domain  .S|  |  response  waveform.  The  Low  Pass  mode  will  give  the 
response  of  the  device  to  either  a  Step  or  an  Impulse  stimulus.  (Mathematically, 
the  Low  Pass  Impulse  stimulus  is  the  derivative  of  the  Step  stimulus.) 

T  hese  Time  Domain  i espouses  were  generated  tising  the  Circuit  Modeling  Pro¬ 
gram  which  is  supplied  with  the  Time  Domain  option  (described  at  the  end  of  the 
1  ime  Domain  section). 
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LONN  PASS  Ki:i  I.IX'I  IONS 
(KLAL  I  ORMA  I  ) 


ELEMENT 

STEP  RESPONSE 

IMPULSE  RESPONSE 

/ 

A 

OPEN 

UNITY  REFLECTION 

UNITY  REFLECTION 

SHORT 

\ 

. V 

UNITY  REFLECTION.  -180- 

UNITY  REFLECTION  -180- 

A 

RESISTOR 
R  >  Zo 

POSITIVE  LEVEL  SHIFT 

POSITIVE  PEAK 

RESISTOR 

\ 

f 

V 

R  <  Zo 

NEGATIVE  LEVEL  SHIFT 

NEGATIVE  PEAK 

l.'iw  l\i'^  Sic;'  end  liiiiuif^c  Mc^i’nn^c  iiycfnrni\  i  I  nl  J" 


2n{}  I  iiiK'  lX)iti.iin  \ k'.isiircimnis 
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■ - - 

POSITIVE  THEN  NEGATIVE  PEAKS 


NEGATIVE  THEN  POSITIVE  PEAKS 

Fi'^urc  Sii.  L(i\y  Sw,”  jiui  fnipul>d  Rc^poii^c  ii'uvcl orni.s  >2  o!  2 
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Reflection  Measurements  in  Time  Domain  Lo'»  Pass 

l  o  make  mcasuiements  the  Low  Pass  mode,  use  the  following  procedure: 
.  Press  PRESET. 

•  Cal,  CAL  I  (  7  mm  )  or  CAL  2  (  3.5  mm  ). 

The  CnI  Type  menu  (Figure  51.  p.  183)  will  be  displayed. 

.  Press  SET  FREQ.  (LOW  PASS). 

•  Perform  an  Sj  j  1-PORT  calibration. 

•  Connect  a  25  .0  airline  and  broadband  load. 

.  Press  DOM  A  I, N.  TIME  LOW  PASS.  SET  FREQ.  (LOW  PASS). 

•  Pi  ess  A I  TO  to  \  iew  the  STEP  response.  Figure  81. 

•  To  view  the  Low  Pass  Impulse  response  of  the  drw  ice.  i^ress  DO¬ 
MAIN.  SPECIF\  Tl.ME.  LMPULSE  (LOW  P  \.SS). 


CO***  •» 
•  -  3  3  3  .  *9  "'w. 


S.;  a. 

®cc  a. a 

laa  a  < «/ 
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Intorproting  the  Low  I’ass  Rosiionse  Horizontal  Axis.  Tite  horizontal  axis  for 
ihc  Low  Pass  nieasorciiient  is  the  2-\vay  ira.xcl  time  to  the  cliscontiiiuity.  the  same 
as  I'oi'  ilic  Band  Pass  mode.  Also,  the  .Marker  rmiction  displays  both  the  time  (\2) 
and  eleetiic.il  length  (x2).  obtained  by  iiuiltiplying  the  time  by  the  veioeity  of 
light  in  a  s  aciium  (2.9b7925LS  nr  sec),  d'o  get  the  actual  physical  length,  multiply 
by  the  relatixe  \  elociiy  of  light  in  the  propagation  medium. 


Interpreting  tlie  Lo'v  Bass  Response  Vortical  .Axis.  The  vertical  axis  depends 
upon  the  t'onv^ii  chosen.  In  the  Low  Pas.s  mode,  the  most  useful  forir.ai  is  RE.AL. 
whicli  displays  the  TDK  response  in  reflection  coefficient  units. 

This  points  out  a  key  difference  between  the  Btind  Pass  and  Low  Pass  modes. 
The  Band  Pass  measurement  is  actually  the  response  of  the  device  to  an  RF  pulse 
with  an  impulse  shaped  envelope.  For  Band  Pass,  the  Inverse  Fourier  Transform 
of  the  (comple.x)  Freciuency  Domnin  data  gives  a  complex  (real  and  imaginary 
parts)  Time  Domain  response,  and  it  is  the  magnitude  of  this  response  that  is 
displa  y  cd. 

In  the  Low  Pass  mode,  because  the  Frequency  Dom.iin  data  is  taken  at  haianoni- 
c.illy  lelated  frctiuencies  dovvn  to  dc.  the  Inverse  Fouiier  Ti'ansform  has  only  a 
retil  part  (the  imtigintiry  part  is  zero).  Therefore,  the  most  useful  format  for  the 
Low  Pass  mode  is  the  REAL  format,  which  displays  the  response  in  reflection 
coefficient  units.  Other  useful  fornuits  .ire  listed  in  Table  15. 


15.  Useful  Tu)k’  Domuiii  Uw  Pass  Formula 


1  FORMAT 

i 

PARAMETER  1 

i  REAL 

!  RLFLECTIOiN  COEFFICIENT  UNITS  i 

!  LOG  MAG 

i  RETl  RN  LOSS  (clB)  i 

SMR 

'  1 

S\A  R  L  N  ITS 

Trace  Bounce.  Depending  on  the  mtigniiude  of  the  response  and  on  the  test  set 
used,  the  Low  Pass  .Step  lesponse  of  ihe  vlevice  may  exhibit  a  phenomenon  called 
display  tr.'ice  bounce.  This  is  normal,  and  it  can  be  impioved  by  turning  on 
averaging  iunder  the  "Response  MENl  ).  This  iiace  bounce  is  caused  by 
loss  of  measurement  dvnamic  range  at  loxv  fiequencies  because  ot'  the  roil  off  of 
tne  coupler-based  test  sets  TIP  85I2A  .md  HP  I(5I4A)  below  500  MHz  'down  -20 
dB  .u  45  Ml-lz).  Tile  tiace  bounce  is  a  factor  ol'  20  times  less  in  the  bridge-based 
lest  sets  (HP  851. 2A  and  HP  <S515A).  whicli  have  flat  magnitude  frequency 
I  espouses  down  to  45  MHz. 
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\b  a  second  example  ol'  L.ow  Pass  reflection  measurements,  consider  tlie  Low 
Pass  Step  response  of  a  30  cm  airline  and  fixed  load,  sliown  in  Figure  82. 


n« 


1- 


F. ptu'c  d2.  Sic;<  Fc'/'n/ixc'  ol  a  JO  cm  Airline  and  Fixed  Load 


The  Lose  Pass  response  at  t  =  0  is  that  of  the  airline  connection.  By  comparing 
this  response  witli  the  ilieoreticnl  Low  Pass  responses,  one  can  determine  whether 
the  :nismatch  pteseiu  is  capticitixe  or  inductive.  The  discontinuity  at  the  first 
connection  of  the  aii'line  is  capacitive.  TIte  upward  slope  of  the  center  section  of 
the  response  is  ctiused  by  the  loss  in  tlie  airline.  The  second  major  response  is  that 
of  the  fixed  load. 
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MASKING 


Marking  is  a  pinsieal  phenomenon  in  which  tiie  Impulse  or  Step  response  of  one 
discominiiii y  affects  ilie  response  of  eacli  subsecinenc  Uiscoiuiniiity  in  the  circuit. 
This  occurs  because  the  enetgy  reflected  from  or  absorbed  in  the  first  discontinu¬ 
ity  ne\er  retiches  the  second.  In  tlie  25  n  airline  e.xaniple  (Figure  81).  the  Low 
Pass  step  resDonse  shows  the  reflection  coefficient  at  the  first  discontinuity  of 
-0. which  is  correct  for  an  imped.ince  of  25  n.  However,  at  the  end  of  the  25 
A.  section  the  response  does  not  return  to  zero  reflection  coefficient,  which  it 
siiould  at  a  50  P.  impedance.  The  reason  is  that  the  step  incident  on  the  second 
response  is  of  less  thtin  unity  amplitude  because  of  the  energy  reflected  in  the 
first  misitiatch. 


As  a  second  example  of  mtisking.  consider  tlie  Time  Domain  response  of  a  3  dB 
attenuator  and  a  short  circuit.  The  Impulse  response  of  the  short  circuit  alone. 
Figine  8.3.  shows  a  return  loss  of  0  dB.  Hosveser.  the  response  of  the  short  circuit 
placed  at  tlie  end  of  the  3  dB  attenuator  displays  a  return  loss  of  -6  dB.  This 
rtilue  actualh'  represents  the  forward  and  leturn  path  loss  through  the  attenuator, 
and  it  illustrates  how  a  lossy  network  can  affect  the  responses  that  follow  it. 
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File  IIP  8510  has  n  feature  called  W  INDOWING  that  is  designed  to  enhance 
Time  Domain  measurements.  I'lie  need  for  Windowing  is  due  to  the  abrupt  tran¬ 
sitions  in  the  Frequency  Domtiin  metisurement  tii  the  START  and  STOP  frequen¬ 
cies.  This  band  limiting  of  the  Frequency  Domain  response  causes  overshoot  and 
ringing  in  the  Time  Domtiin  response.  It  causes  the  (un-Windowed)  Impulse 
stimulus  to  have  a  sintkt),  kt  shtipe  (k  =7r/ frequency  span),  whicli  has  two  effects 
that  iiinit  tin.  useftdness  of  t!te  Time  Domain  measurement: 

(1)  Finite  Impulse  Width.  Tliis  limits  the  ability  to  resolve  between 
two  closely  spaced  responses.  The  effects  of  the  finite  impulse  width 
cannot  be  improsed  without  increasing  the  fi'equency  span  of  the 
measuiement.  .See  Table  16. 

(2)  Sidelobes.  The  Impulse  sidelobes  limit  the  dynamic  range  of  the 
Time  Domain  measurement  by  hiding  low  level  responses  within  the 
sidelobes  of  the  higher  le\  cl  responses.  Tiie  effects  of  sidelobes  can  be 
improsed  by  Windowing.  See  Table  17. 

W  indowing  impioves  the  dynamic  range  of  the  Time  Domain  measurement  by 
modifying  (filtering)  the  Fietiuency  Domnin  datti  prior  to  conversion  to  the  Time 
Domtiin  to  produce  an  impulse  stimulus  with  lower  sidelobes.  This  greatly  enhan¬ 
ces  the  effectis eness  in  viewing  Time  Domain  responses  that  are  very  different  in 
magnitude.  The  sidelobe  reduction  is  achies'ed.  hoss-ever.  as  the  tradeoff  with 
increased  impulse  width.  The  effect  of  Windowing  on  the  STEP  stimulus  (inte¬ 
gral  of  the  impulse  stimulus.  Losv  Ptiss  ihode  only)  is  a  reduction  of  overshoot 
and  ringing  at  the  tradeoff  with  incretised  rise  time. 

Three  Windows  are  asailable:  MINI.MfvM.  NORMAL,  and  MA.XIML'.M.  The 
Window  may  be  selected  by  pressing  DOM.-\IN.  .SPEC1F^  TIME.  The  sidelobe 
Ie\  els  of  the  Time  Domain  stimulus  depend  only  on  the  Window  that  is  selected 
(see  Table  1“).  .MINIMIM  is  esseniially  no  window  and  therefore  gives  the 
highest  sidelobes;  NOR.MAL  (selected  b>  PRESET)  gives  reduced  sidelobes  and  is 
normally  the  most  useful:  M.A.XI.VILAI  gives  the  minimum  sidelobes  and  thus 
pros  ides  the  greatest  dy  namic  range. 
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I 

I 


] 


Impulse  Sideiobe  Level 


WINDOW 

TYPE 


IMPULSE  STEP 

SIDELOBE  SIDELOBE 
LEVEL  LEVEL 


Step  Sideiobe  Level 


MINIMUM 

NORMAL 

MAXIMUM 


-13  dB  -21  dB 

--+4  dB  -|}0  dB 

<  -90  dB  <  -90  dB 


The  sideiobe  reduction  due  to  Windowing  is  achie\ed  nt  n  tradeoff  witit  an  in¬ 
crease  in  the  Step  (10%  -  90%)  Rise  Time  and  tiie  Impulse  (50%)  width.  These 
parameters  also  depend  upon  the  frequency  span  of  the  measurement,  and  they 
can  be  calculated  using  the  approximate  forniulas  given  in  Table  17. 


Tahie  /".  .4  p pro.xiniaie  Formulas  For  Swp  Rise  Time  and  Im/nilsL'  IL.T/t/t 


LOW  PASS 


STEP  RISE  Tl.ME 
( I  O^b  -  90'>b) 


l.M PULSE  WIDTH 
'50%) 


(  1.0  .MINI.VIU.M  WINDOW  I 
{  I.D  NOR.MAL  WINDOW  ; 
{  2.4  .MAXIMUM  W  INDOW  ! 


BAND  PASS 


IMPULSE  W  IDTH  1.20 

'50'^o)  'FREQSP.A.N 


{1.0  MINIMUM  WINDOW 
0.45  V  (  2.2  NORMAL  WINDOW- 

FREQ  SPAN"  {  3.3  MAXIMUM  WINDOW’ 

{  1.0  MINIMUM  WINDOW 
0.60  V  (1.6  NOR.MAL  WIN  DOW-’ 

FREQ  SPAN  (  2.4  .MAXI.MU.VI  W’INDOW’  I 


Alultiply  by  the  velocity  of  light  in  a  vacuum  (2.997925E8  m.sec)  to  get  electrical 
length,  and  then  by  the  relative  velocity  of  light  in  the  the  propagation  medium 
to  get  physical  length. 
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The  purpose  of  windowing  is  to  make  the  Time  Domain  response  more  useful  in 
isolating  and  identifying  individual  responses.  The  window  does  not  affect  the 
displayed  Frequency  Domain  response.  It  is  turned  on  only  when  the  Time  Do¬ 
main  response  is  viewed.  Figure  84  shows  typical  effects  of  windowing  on  the 
Time  Domain  response  of  the  reflection  measurement  of  a  short  circuit. 


WINDOW  MINIMUM  NORMAL  WIDE 


LOW  PASS 
STEP 


LOW  PASS 
IMPULSE 


HAND  PASS 
IMPULSE 


Fi^'urc  iS4.  ElU'c!  ol  li'iihlnu-tni:  mi  Time  Domain  Responses 

ol  a  Shori  Circiiii 
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In  the  Tinic  Domal/i.  (he  RANG!:  is  tlcrincd  as  ihe  loiiglli  in  lime  that  a  measure¬ 
ment  can  be  made  witliout  encouiueiing  a  repeiition  of  the  response  (see  Figure 
83).  1  lie  repetiiion  ol'  (he  Time  Domain  response  occurs  al  regular  intervals  of 
time  and  is  a  conseciuence  of  the  Frequency  Domain  data  being  taken  at  discrete 
freciuency  points  rather  than  being  continuous. 

The  Range  of  ;i  measurement  is  equal  to  l/AF.  the  spacing  between  freciuency 
data  points.  It  is  therefoie  directly  proportional  to  the  number  of  points  and 
iinersely  piopoiiional  to  the  Frequency  Span  (STOP  -  START  freciuency)  and 
can  be  calculated  using  the  follow'ing  formula. 

R.AiNGF  =  l/AF  =  (Number  of  Points  -  1 )/ Frequency  Span 

As  a  sample  calculation,  for  a  201  point  measurement  from  50  MFlz  to  18  GHz 
(S1*AN  =  I  7.d5  GHz),  the  Range  is  (201  -  I)  17.95  GHz  =  II. I  nsec  (3.34  m). 
Thus  the  device  under  test  has  to  lie  3.34  m  or  less  in  electrical  length  for  a  trans¬ 
mission  measurement  (1.67  m  for  a  reflection  measurement)  or  else  an  overlap¬ 
ping  of  the  Time  Domain  responses  (aii.'ising)  \s  ill  occur.  (Remember  to  multiply 
by  the  lelative  selociiy  of  light  in  the  medium  to  get  tictual  physical  length.) 


/•(C’lnc  (V5.  riiiK'  Doniain  McuMin'iiu'iu  Shnwnn;  Rcsiyonw  Rcpcitiidn.^ 


To  increase  the  Time  Doniain  measurement  Range,  it  is  usually  better  to  first 
inciease  the  number  of  points,  because  decreasing  the  freciuency  span  will  leduce 
the  resolution. 
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RLSOLl  TION 


['here  aie  two  different  terms  involving  resolution  in  Time  Domnin;  RESPONSE- 
RESOLLTION  and  RANGE-RESOLUTION.  The  Time  Domnin  Response- 
Resolution  is  defined  ns  the  nbility  to  resolve  two  closely  spaced  responses.  In 
other  woids.  if  two  responses  are  present,  this  is  how  closely  they  can  be  spaced 
and  still  be  distinguished  from  one  another.  For  responses  of  equal  amplitude,  the 
Response-Resolution  is  equal  to  the  50%  (-6  dB)  impulse  width.  It  therefore  is 
inversely  proportional  to  the  frequency  span  of  the  measurement  and  is  also  a 
ftmetion  of  the  window  that  is  tised.  .Approximate  formulas  for  calculating  the 
50'-’u  Impulse  width  are  gi\’en  in  Table  17.  For  responses  that  are  of  different 
.impiitudes.  the  Response-Resolution  will  be  wider. 

Range-Resolution  is  defined  as  the  ability  to  locate  a  single  response  in  time.  In 
other  words,  if  only  one  response  is  there,  this  is  how  closely  you  can  pinpoint  the 
peak  of  that  response.  The  Range-Resolution  is  equal  to  the  digital  resolution  of 
the  CRT  display  which  is  the  time  span  displayed  divided  by  the  number  of 
points.  Maximum  Range-Resolution  is  achieved  by  centering  the  response  on  the 
display  and  then  reducing  the  time  span.  Therefore,  the  Range-Resolution  is  al¬ 
ways  much  finer  titan  the  Response-Resolution. 

lo  illustrate  the  difference  between  these  two  resolution  terms,  consider  a  mea¬ 
surement  with  a  frequency  span  of  18  GHz.  For  Low  Pass,  with  a  Normal 
■'A  indow,  the  Response-Resolution  (Impulse  width)  is  53  psec  (0.6  x  (1  18  GHz)  x 
i.o;  or  16  mm  in  electrical  length  (53  psec  x  2.997925E8  m/sec).  .As  illustrated  in 
Figure  8b.  two  Time  Domain  responses  of  equal  amplitude  separated  by  16  mm 
could  be  resolved  in  (his  Time  Domain  measurement.  (This  indicates  an  actual 
discontinuity  sepaiation  of  8  mm  for  reflection  measurements.)  ’ 

Now  consider  the  case  where  only  one  response  is  present.  Bx’  centering  that 
lesponse  on  the  display  and  adjusting  the  time  SPAN  to  equal  the  50%  Impulse 
wiqth  (53  psec.  16  mm).  Figure  86.  the  Range-Resolution  is  reduced  to  40  um  1 16 
mm  401  points).  The  Range-Resolution  can  be  further  reduced  by  narrowing  the 
time  span. 
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GATING 

The  HP  8510  gating  featuie  gives  the  user  the  flexibility  to  selectively  remove 
reflection  oi'  transmission  Time  Domain  responses.  In  converting  back  to  the 
Fiequency  Domiiin.  the  effects  of  the  responses  outside  the  Gate  are  removed.  In 
a  reflection  measurement,  you  can  remove  the  effects  of  unwanted  mismatches 
or  else  isolate  and  view  the  response  of  an  individual  mismatch.  In  a  transmission 
measurement  you  can  remove  the  responses  of  multiple  transmission  paths. 


Setting  the  Chite.  A  Gate  is  a  temporal  band  pass  filter  used  to  filter  out  unwan¬ 
ted  Time  Doniiiln  responses.  Responses  outside  the  selected  gate  are  not  included 
in  the  trtice.  Tliere  are  three  Gate  indicators:  ST.ART.  CE.NTER.  and  STOP. 
The  Gate  has  a  bandpass  filter  shape,  as  shown  in  Figure  87.  The  GATE  CEN¬ 
TER  indicates  tlie  center  time  (not  frequency)  of  this  filter,  and  the  Gate  ST.ART 
and  STOP  indicate  the  -6  dB  cutoff  times.  Gate  SPAN  =  STOP  -  START. 
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Cousitlci  iisiny  uaiiny  li)  .inaly/c  iltc  response  ol'  a  7inm-lo-3.5mni  adaptei'  con- 
necied  lo  a  .v^miii  airline  and  a  fixed  load.  I  lie  Fretiiiency  Domain  and  ihe 
Hand  Pass  l  ime  Domain  responses  of  such  a  setup  are  shown  in  l"igure  88. 
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W  c  will  now  use  gating  to  analyze  tlie  response  ol'  the  adapter  onl>, 

•  Press  DOMAIN,  riMlI  HAND  PASS.  SPEC! l•^  GATE. 

•  The  three  Gttte  indicators  will  now  appear  on  the  screen.  Press 
GATE  CENTER,  and  use  die  knob  or  ke\pad  to  mo\e  the  center 
indicator  to  t  =  0. 

In  Figure  80.  the  time  domain  display  shows  the  gate  center.  86  ps.  as 
the  Actise  Function. 

•  Press  GATE  SP.A.N  and  use  the  knob  or  keypad  to  adjust  the  Gate 
.Span  to  0.“0  ns. 

•  Pre.ss  G.ATE  ON  to  turn  on  the  Gate. 

The  responses  outside  the  G;tte  will  be  removed.  See  Fiaure  89. 

•  Press  DOMAIN.  FREQL  ENCA. 

View  the  gated  Freciuency  Domain  response  of  the  adapter.  See 
Figure  89. 
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The  dtirker  shaded  trace  in  the  Frequency  Domain  plot  of  Figure  89  shows  the 
Gated  Freciueticy  Domain  response,  which  is  that  of  the  adapter  only.  The  ef¬ 
fects  of  the  fixed  load  ofi  the  measurement  are  reniosed. 
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(..itc  .Sli;i[)e.  lour  clin'ciciK  Ci.iii.-  sli.ipcs  ,iic  ;)\;ii  In  bio;  MINIMIM, 
\()RM\1.  WlDl..  ,iik1  MA.MMl  ^l.  I  noli  ol'  the  C'intes  lia\o  (.lirrcicnt  pass- 
a,;;Kl  li.iiiioss.  ouiol'f  laic.  aiul  bitlolobo  lo\olv  'l  l  iiKlicaios  the  Gate  s|)an  which 
IS  ih.o  lane  iioiuocii  ihc  (aiio  si.m  anil  sU)c>  liulica(ois.  TZ  is  iho  time  botwocii  the 
eske  sjl'  ilic  (j.iio  p.issbantl  and  the  -b  dll  Gale  si(.>p  time.  T3.  eciutil  to  12.  is  ihio 
lime  bei'Aceii  ihe  Gale  slop  time  and  (he  point  wlieie  (he  filter  first  reaches  the 
iesei  of  ihe  iiiehesi  Gate  sielelobe,  1  lie  Gale  eharaeterisiics  for  eticii  Gate  shtipe 
are  hsicsi  in  I'abie  1 8. 
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!  iie  Passoanel  Ri()ple  and  Sielelobe  Lesels  are  deseriptise  of  the  gate  (filter)  shape. 
I  lie  C  uiofr  Time.  '1  2  =  1."'  isee  I'al.ile  18).  indietiies  how  fast  the  gate  filtei'  I’oils 
k'l  r.  1  01  eaeii  gate  shape,  iliere  is  also  a  Mininuini  Gate  Sptin  lTl,jijn  =  -  "^-1 

hid)  gr.  es  a  filter  p.issband  eif  zero,  l  o  enter  a  Gtiie  sptin  smaller  thtin  mini- 
ainm  '.sill  iirodnce  a  disioried  fillet  shtipe  ilitii  will  htise  no  passband.  will  not 
iiiise  a  narrower  slitipe.  mtis  htise  higher  sielelobe  lesels.  tind  svill  gise  an  iiicor- 
leci  inehcaiion  of  gtne  Sl  XR'l  aiiel  S'lOP  times,  'riierefore.  it  is  important  to 
alssass  select  .1  Cj.iie  sistin  ill. it  is  higher  ih.in  the  mimimim  salne.  The  Cutoff 
'line  .me!  the  Miiiinmm  G.ite  Span  tire  iiisersels  proportiontil  to  the  frequency 
sp.m  elf  tile  metisni  emeiil  .is  mehe.iteel  in  I.ilile  18. 

Sesi  icsniis  using  Gating,  it  is  import.int  to  tilsstiss  eeiiier  the  G.ite  around  the 
espeinsc'si  that  son  ssani  to  retain  m  the  me.isui enieiii  .ind  to  make  the  Gate 
'■o.m  ss  leie  eiiongh  to  ineluele  .ill  eif  those  1  espouses.  It  is  .ilso  recomnieiided  to  use 
'lie  ss  lelest  Cj.iie  sli.ipe  peissinle. 
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Operaliiiy  aiul  Prograinmiiii' 


MP  \SL  RP.MPM'  RPCO.M.MPNDATIONS 

W  iicii  making  Time  Domain  measmeiueiUb.  it  is  genei'all>  a  good  practice  to 
mcastire  die  device  within  the  rrequency  range  that  it  is  designed  to  opei'ate. 
There  are  two  leasons  for  this.  First,  tlie  noise  floor  of  the  Time  Domain 
lesponse  is  directly  related  to  the  noise  in  the  Frec|iiency  Domain  data.  There- 
t'ore.  if  mtiny  of  the  Frequency  Donuiin  data  points  are  taken  at  or  below  the 
noise  floor  ol'  the  measurement,  then  the  noise  floor  of  the  Time  Domain  mca- 
suiemeni  will  be  increased.  A  second  reason  to  measure  the  device  with.in  its 
ooeiaiinii  fietniencs'  range  is  because  the  in  band  response  is  normally  of  interest, 
r  le  Time  Domain  metisurement  is  an  average  of  the  response  over  the  frequency 
r.mge  ot'  the  metisurement.  .ind  if  the  Frequency  Domain  data  is  measured  out  oi' 
iMiid.  tiien  tile  Time  Domain  measurement  vviil  tilso  be  the  out  of  band  response. 
However,  since  tlie  Time  Domain  Response-Resolution  is  inversely  proportional  to 
the  fieciuency  span,  it  may  at  times  be  desirable  (with  these  limitations  in  mindj  to 
use  a  fi'etiuency  span  that  is  slightly  wider  than  the  device  bandwidth  to  give 
better  resolution. 


Source  Considerations 

•\  tnough  eithei-  source  vviil  work  well  in  making  Time  Domain  measurements,  the 
biP  8.'A0a  svmliesized  sweeper  has  the  advantage  that  it  provides  greater 
d;  mimic  .iinge  than  the  HP  8350Q  sweeper.  The  main  reason  for  this  is  the 
f  eciueiicy  sttibilitv  of  a  synthesized  source.  Tlie  small  nonlinearities  and  phase 
discontinuities  that  occur  in  the  ramp  sweep  mode  cause  low  level  noise  sidebands 
on  the  Time  Donuiin  Impulse  or  Step  stimulus.  These  interfere  in  measurements 
reviuiring  large  dynamic  range.  Perform  a  TRIM  SWEEP  adjustment  before 
caiibi iuing  to  help  ininimize  these  noise  sidebands.  Adjusting  trim  sweep  is  ex- 
pianned  at  the  end  of  the  section  of  tliis  manual  titled  .Measurement  Calibration. 

la  the  IIP  8,'-iU.A  (synthesized)  step  svveep  mode,  the  improvement  in  source 
vi,,iiiiitv  eimiin.nes  these  noise  sidebands  and  improves  the  Time  Domain  measure- 
ue;u  dynamic  rtinge  iiy  as  much  .is  .'’0  dB.  A  second  improvement  is  that  the  HP 
v'-OA  steppevi  sweep  mode  allows  the  use  of  many  averages  per  point  without 
zw.iiiy  affecting  the  sweep  time,  and  this  lowers  the  noise  floor  of  the  Time 
Domain  measurement.  It  is  recommended  to  perform  a  Trim  Sweep  adjustment 
01 'or  to  calibrating  when  making  measurements  in  the  ramp  svveep  mode  to  mini- 
'lare  pn, ise  (.(isconi muities. . 
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Test  Set  Considerations 

The  bridge-based  test  sets  (HP  851  3A  and  HP  851  5A)  have  two  advantages  o\er 
the  coupler-based  test  sets  (HP  8512A  and  HP  8514A)  when  making  Time  Do¬ 
main  measurements.  First,  the  bridge-based  test  sets  extend  in  frequency  to  26.5 
Gblz.  versus  18  GHz  foi'  the  couplei-based  test  sets.  When  measuring  btoadband 
devices,  this  extia  bandwidth  piovides  better  Time  Domain  Response- Resolution. 

The  second  advantage  is  that  the  bridge-based  test  sets  have  a  flat  response  down 
to  45  MHz.  whereas  the  coupler-based  test  sets  begin  to  roll  off  (but  are  still 
usable)  below  500  MHz.  This  coupler  roll  off  I'educes  the  dynamic  range  avail¬ 
able  at  the  low  frequencies  (-.30  dB  at  45  .MHz)  and  therefore  incieases  tlie  Time 
Domain  noise  flooi’  when  measurements  are  made  at  those  frequencies  (this  caus¬ 
es  the  trace  bounce  in  the  l.ov,  Pass  Step  response). 
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